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PREFACE 

THE concept of the boundary layer, introduced by Ludwig Prandtl in 
1 904, has been a particularly fruitful one. Research on this topic has 
now reached the stage where there is a certain body of fundamental 
definitive information which is unlikely to be superseded to any great 
extent. This relates mainly to the steady incompressible laminar 
boundary layer in two dimensions. Work proceeds, however, and a 
considerable number of papers are being published, on unsteady 
boundary layers, on three-dimensional boundary layers in incompres- 
sible flow, and upon such topics as boundary-layer stability. In com- 
pressible flow, too, where additional important parameters arise, much 
is being done and more remains. 

This monograph is one of a series, each of which is being written by 
an author on the general field in which his own research interests lie. 
It is inevitable, therefore, that there will be a certain amount of bias 
in the choice of material. I have tried to make the book reasonably 
self-sufficient, though lack of space has led to the omission of a number 
of very interesting problems of boundary layers. The topics so axed 
include unsteady boundary layers and boundary-layer stability, bound- 
ary layers on porous walls with suction or blowing, boundary layers in 
three dimensions (including axi-symmetric flow) and boundary layers 
with vorticity in the mainstream. For a discussion of these topics 
reference may be made to more encyclopaedic works, such as Modern 
Developments in Fluid Dynamics (Oxford, edited by S. Goldstein), the 
companion volumes on High Speed Flow (edited by L. Howarth), 
Laminar Boundary Layers (Oxford, edited by L. Rosenhead), and 
Volumes III to V of the series High Speed Aerodynamics and Jet 
Propulsion. 

The purpose of Chapter 1 of this book is to introduce the boundary- 
layer concept, and to show how the equations of viscous flow are 
simplified thereby. The standard boundary -layer parameters, and the 
usual integral forms of the boundary-layer equations, are discussed, 
the incompressible flow forms being introduced as special cases of the 
more general compressible forms. Chapters 2 to 6 deal with various 
aspects of solutions in incompressible flow, commencing with analytic 
solutions for the velocity field, these being solutions which may be 
expressed in terms of functions which satisfy ordinary differential 



vi PREFACE 

equations (Chapter 2). There follow discussions of high-accuracy nu- 
merical solutions for the velocity field (Chapter 3), practical methods 
of calculation (Chapter 4), and an analysis of the factors which might 
govern the choice of a method (Chapter 5). Various types of solution 
of the temperature equation in incompressible flow with small tem- 
perature differences are considered in Chapter 6. Chapters 7 to 9 deal 
with compressible laminar boundary layers, consideration being first 
given to flow with zero pressure gradient (Chapter 7), then to flow with 
zero heat transfer (Chapter 8), and then to flow in which both pressure 
gradient and heat transfer are present (Chapter 9). Finally there is a 
brief discussion (Chapter 10) of some aspects of the problem of the 
interaction between shock waves and laminar boundary layers. 

It is my hope that this book will be of value to a wide variety of 
workers. In the first place I have tried to present the material in a 
sufficiently ordered and logical manner as to make it of value as an 
introduction to boundary -layer theory for young research workers who 
are new to the subject, or to undergraduates who are familiar with the 
elements of classical in viscid fluid dynamics. Secondly, the book should 
be of some value to research workers in this field, since one of the things 
which has governed my choice of material has been the question of 
whether a particular piece of work has been an end in itself or whether 
it has assisted in opening up the way for further advances. Finally, 
I have borne in mind the needs of practising engineers, and have tried 
where possible to indicate the limitations, the likely accuracy, and the 
practical complexity of the methods described for calculating the various 
properties of laminar boundary layers. 

In conclusion, I have great pleasure in expressing my thanks to the 
many people who have helped me, directly or indirectly, in the writing 
of this book. To Dr. M. J. Lighthill, F.R.S., Director of the Royal 
Aircraft Establishment, whose student I was at Manchester University, 
for the wise counsel he gave me then in so many branches of fluid 
dynamics, and whose influence is, I hope, evident in this book. To 
my colleagues at the National Physical Laboratory, for the stimulating 
discussions I have had with them at various times, and most particu- 
larly Dr. J. T. Stuart and Dr. G. E. Gadd. These two colleagues have 
been good enough to offer useful comments on a first draft of this 
book, although the responsibility for its deficiencies remains entirely 
my own. To Mrs. M. E. M. Sayer, for her patient and careful typing 
of the manuscript, and her cheerful approach to the difficult task: of 
reading my writing. To Professor G. Temple, F.R.S., editor of this 
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series, and the staff of the Oxford University Press for the courteous 
way they have dealt with the various problems which have arisen. 
To Sir Gordon Sutherland, F.R.S., Director of the National Physical 
Laboratory, for permission to write this book. The writing has not in 
fact been part of my official duties, and the views expressed are entirely 
my own. Thanks are also due to various bodies for permission to use 
copyright material such as charts and tables. These bodies include 

the Aeronautical Research Council, 

the Clarendon Press, Oxford, 

the Controller, H.M. Stationery Office, Prof. L. Crocco, 

the Institute of the Aero/Space Sciences, 

the Editor, Journal of Fluid Mechanics, 

the Director, National Physical Laboratory, 

the Royal Aeronautical Society, 

the Royal Society, 

the United States Air Force, 

the Editor, Zeitschrift fur angew. Math, und Mech. 

Finally, to my wife and, though they know it not, to my children, for 
so ordering their lives as to make the task of writing this book much 
less difficult. To all these people go my sincere thanks for their much 
appreciated help. 

N. C. 

Hanworth, Middlesex 
June 1961 
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INTRODUCTION 

THIS book is about solutions of the laminar-boundary -layer equations. 
The concept of the boundary layer, one of the corner-stones of modern 
fluid dynamics, was introduced by Praiidtl (1904) in an attempt to 
account for the sometimes considerable discrepancies between the pre- 
dictions of classical inviscid incompressible fluid dynamics and the 
results of experimental observations. As an example, we may remark 
that according to inviscid theory any body moving uniformly through 
an unbounded homogeneous fluid will experience zero drag! 

Now the classical inviscid theories assume that the viscous forces in 
a fluid may be neglected in comparison with the inertia forces. This, 
indeed, would seem a reasonable approximation, since the viscosity of 
many fluids (and of air in particular) is extremely small. However, in 
certain regions of flow, fortunately often limited, the viscous forces can 
still be locally important, as Prandtl observed. The reason for this is 
that* a typical viscous stress is of magnitude p(du/c)'i/), where (JL is the 
viscosity, u is the velocity measured in a direction parallel to that of 
the stress, and y is distance measured normal thereto, so that when 
the velocity gradient (or shear) c)u/dy is large the viscous stress can 
become important even though p, itself is small. It was Praiidtl who 
remarked that in flow past a streamlined body, the region in which 
viscous forces are important is often confined to a thin, layer adjacent 
to the body, and to a thin wake behind it. This thin layer is referred 
to as the boundary layer. When this condition holds the equations 
governing the motion of the fluid within the boundary layer take a 
form considerably simpler than the full viscous-flow equations, though 
less simple than the inviscid equations, and it is the solution of these 
equations with which we shall be presently concerned. 

An alternative method of looking at this concept is as follows. The 
inviscid-fiow equations are of lower order than the viscous-flow equa- 
tions, so fewer boundary conditions can be satisfied in a mathematical 
solution of a given problem. Thus an inviscid -flow solution allows a 
finite velocity of slip at a solid boundary, whereas the solution of the 
viscous-flow equations does not allow such slip. In other words, the 
inviscid-flow solution assumes the existence of an appropriate vortex 
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2 INTRODUCTION 1 

sheet at the solid boundaries. Now in reality this vorticity will diffuse 
outwards from the boundary (in much the same way that heat would 
diffuse from a heated body) and will be convected with the stream. 
Thus, considering now the flow past a flat plate, the time t in which 
fluid travels a distance x parallel to the plate will be of order x/U, where 
U is a typical velocity, and in this time the vorticity will have diffused 
outwards through a distance of order (vt)* = (vx/U)*, where v ^ p,/p is 
the kinematic viscosity. This is an indication of the boundary-layer 
thickness. 

Before turning to a more quantitative discussion, mention must be 
made of the phenomenon of boundary-layer separation. TiWheu the fluid 
is proceeding into a region of rising pressure, it is slowed down by this 
retarding force. In the outer part of the boundary layer, where the 
kinetic energy is large, this results only in a relatively small slowing 
down of the fluid, but the effect on the slower-moving fluid nearer to 
the wall can be considerable, and if the pressure rise is sufficient it can 
be brought to rest, and, farther downstream, a slow back-flow be set 
up. In such circumstances the forward flow must leave the surface to 
by-pass this region, and boundary -layer separation is said to have taken 
placeJ If the region of separated flow is extensive, the separation can 
have a back-reaction on the external flow, which is then quite different 
from what it would have been in the absence of the boundary layer. 
If the separated region is limited, on the other hand, the external flow 
may not be significantly affected, and the flow field may be calculated 
by calculating firstly the external flow (on the assumption of no bound- 
ary layer) and then calculating the boundary layer appropriate to this 
external flow. 

The above qualitative analysis is restricted to low-speed flows, but 
serves to indicate the nature of the boundary layer in a simple way. 
At supersonic speeds, for example, interactions of the boundary layer 
with the external stream become more important, and lead to con- 
siderable theoretical difficulties which will not be discussed at this stage. 
In what follows a quantitative analysis will be given of how the 
boundary-layer equations may be deduced from the exact equations 
of viscous flow, and detailed discussions will be given of some of the 
points touched upon briefly above. 

1. The equations of viscous flow 

We take cartesian coordinates (x,y), with associated components of 
fluid velocity (u, v). The fluid is assumed to have pressure p, density /o, 
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and absolute temperature T, arid these are functions of x and y only, 
in view of the approximation (made throughout this book) that the 
flow is two-dimensional and steady. The equations of motion then 
express the basic physical ideas that for a given element of fluid there 
is conservation of mass, momentum (excepting in so far as the element 
is acted upon by various forces), and energy (excepting in so far as 
work is done by these same forces). For a general derivation of these 
equations, reference may be made to volume I of Modern Developments 
in Fluid Dynamics, High Speed Flow (Howarth, 1953). In steady two- 
dimensional flow the equation of conservation of mass takes the form 

|00+^) = 0. (1) 

This equation is usually referred to simply as the equation of continuity. 
The equations of conservation of momentum in the x and y directions 
(the momentum equations) become 

du . du dp ! #A Y \ fd' 2 u d 2 u 

P U ~~ +pV - = -f- + Jft - "PA +/X -- 

dx dy dx dx \dx- 



" O ' ~~ 'j ' ' O I O I ' ~\ I ' V *"' / 

dx dx dx dy\dy dx] 
and 

f>?> , dv dp , 1 a A . , r . /a 2 ?; d 2 ?;\ 

pU - }- p?' r - = +5^7 +pl +M 7-S + V -T, 

dx rty/ d?/ diy \dx* dy] 

5 A--- +2~- - +7-17-' + 7r-) (3) 
^/ r 1 ?/ c 1 yy dx\dy dx/ 

dn dv- 

where A -+--, (4) 

r'x a?/ 

and JT, l r are the external forces per unit mass of fluid. Finally the 
equation of conservation of energy (the thermal energy equation) takes 
the form 



by J) dx dy dx\ dx dy\ dy 

ij y \/7\// 

(6) 

, . JlduV . /dv\*~] . /dv . duV 2/du 

where 



There are thus four equations, (1), (2), (3), and (5), for the five unknowns 
u, v, p, p, T, and these are soluble, in principle at any rate, when an 
equation of state is defined, relating p, p, and T. For a perfect gas 
this takes the form n =- ( MT (7) 

where $ is the gas constant. 
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In the above equations c. p is the specific heat at constant pressure, 
usually taken to be constant, and k is the thermal conductivity, related 
to the thermometric conductivity K by the relationship 

K = k/pc p . 

For a physical interpretation of the quantities K, k, \L, and v = \i\p, the 
reader is referred to the book Modern Developments in Fluid Dynamics 
(Goldstein, 1938). It will suffice for the present to remark that v is a 
parameter determining the rate at which vorticity is diffused, whilst K 
determines the rate at which heat is diffused. The ratio 

o- -= V\K = fic. 2) /k 

accordingly determines the relative rates of these two types of diffusion, 
and is called the Prandtl number. It is usual to assume that the Prandtl 
number is constant, and this holds trite for air over quite a wide range 
of conditions, the value being about 0-72. It will be seen later that 
considerable simplifications are often possible if it is assumed that the 
Prandtl number is unity, an approximation that is not without value 
for air. 

2. Boundary layers 

It will be noted that the viscosity ^ appears in equations (2) and (3) 
only as a multiplicative factor of velocity gradients, or of powers or 
products of velocity gradients. Accordingly, if the viscosity is small, 
classical fluid dynamics theory, which neglects viscosity, will be valid 
except in regions where velocity gradients are large. Similarly, pro- 
vided the Prandtl number a is not too small, small viscosity implies 
small thermal conductivity, so that the terms involving k in (5) are 
important only where the temperature gradients are large. 

Now it is often found in practice that the regions of high velocity 
and temperature gradients are confined to a narrow region near to solid 
walls, known as the boundary layer, and to a thin wake behind stream- 
lined bodies. In such a domain, considerable simplifications of the 
equations of motion are possible, even though all the terms involving 
jit and k cannot be neglected, as was first shown by Prandtl (1904) in 
a paper of fundamental importance. In low-speed flow it is usually 
possible to consider the development of the boundary layer as a separate 
problem from that of the substantially inviscid external flow, the ex- 
ception being when there is a substantial region of separated flow. 
When the external flow is supersonic, however, there is an interaction 
between the boundary layer and the external stream which must be 
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taken into account. Crudely we may say that though an artificially 
introduced disturbance cannot be propagated upstream in a wholly 
supersonic flow, the presence of a boundary layer, in which the flow 
sufficiently near to the walls will be subsonic, provides a mechanism 
for such upstream influence. Accordingly the external stream does not 
approximate to that obtained in the absence of the boundary layer, 
and cannot be independently prescribed. 

In spite of this difficulty it is still useful to begin by considering 
separately the boundary- layer approximation and the inviscid approxi- 
mation, as use can be made of these results even in certain problems 
in supersonic flow where the interaction between boundary layer and 
external stream is particularly important, as, for example, when a shock 
wave interacts with a boundary layer. 

3. The laminar boundary layer on a plane wall 

We begin by deriving the boundary-layer equations for flow over a 
plane wall. The ;r-axis is taken along the wall and the ?/~axis perpen- 
dicular to it. We assume that the thickness 8 r of the layer in which 
the velocity gradient chi/Py is large, where the velocity u rises rapidly 
from at the wall to a value u ly is much smaller than a typical length I 
in the flow field as a whole. Equally well we assume that 8^, the thick- 
ness of the layer in which the temperature gradient 8T/Py is large, and 
in which the temperature of the fluid changes from the temperature 
of the wall T w to a value 7\, is also much less than I. We shall assume 
at this stage that 8. r and S f are the same order of magnitude, each being 
of order 8 /. This restriction can later be removed, provided 8.,, and 
8, are both /. We let u , p () , 7J, be typical values of velocity u, density 
p, temperature T, and may then deduce the order of v from (1). This 
equation shows that. 

~ C (pr) = -L(pu), 

c)y dx 

which is of order p {} u () jl. Hence, upon integrating across the layer from 
y o to y = 8, we find that pv is of order (p ^ /^)8, so that 



Thus v is small compared with u, within and at the edge of the boundary 
layer. 

We can now examine the magnitudes of the various terms in (2). 
In doing so we remember that the derivative dF/dx of a function F 
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will be much smaller than the derivative dFjdy. In fact 

f 

Thus du du dp 

9U ~d~x pV 8y' fx are 

d 2/ M d/z du 

"**' iylTy are 
and 

, 0A d*u 9A d/x ^Budjji da dv ~ /m /1AX 

IM u |A-^-, 2 --i-, -^- are Oi^uJl 2 ). (10) 

3r ^or r ^ 2 3 c?x ao:^ a?/aa: vro 0/ y 

Ideal fluid-dynamics theory rejects all the terms involving viscosity, 
namely those in (9) and (10), but boundary-layer theory retains the 
former, since S 2 is small as well as /x . We see that these terms are 
the same order of magnitude as (8) when 

), (11) 



p Q u 
where R is the Reynolds number 



Mo ^o 

We have chosen 8 V as the relevant value of 8 since we are considering 
a momentum equation. Upon rejecting the terms in (10), then (2) 
simplifies to 8u 8 8 , u bu 



(12) 

We note that the rejected terms are of order (S/Z) 2 times those 
retained. 

We now deal in a similar manner with equation (3). We can see that 



P u 8 " P v^_ are 0(^F) = Opp) (13) 



$ 2 V dlA $1fc T ^A o A ^Lt rt SV ^LL /-k/M'n *0 w \ /i ^\ 

_c__, lu , *A , 2-7P are ^ ? (14) 

dy z dx oy cy cy oy oy \ o" 5 Lj 

, 3 2 v du, dv /-k/Mn^n 8\ /ie-\ 

and /Ltr, ^--- are Oc^-5-. (15) 

C^iC C73? C/3? \ ' '/ 

Accordingly, the terms in (15) are of order 8/Z times those rejected in 
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(2), and those in (13) and (14) are of order 8/Z times those retained. 
It follows that 

"" 



Now Y will usually be zero. For example body forces can be neglected 
in problems of forced convection, and in problems of free convection 
the body force (gravity) will act in the x direction. Assuming, then, 
that Y = 0, it follows from (16) that the pressure gradient dpjdy is 
small, and the pressure change across the boundary layer is very small, 
being 0(/3 ?^S 2 /Z 2 ), which is neglected. Thus (3) reduces simply to 

p(x,y)=p(x). (17) 

We turn now to the thermal energy equation (5). By identical 
reasoning to that given above it follows that the two terms on the 
left-hand side are of equal order of magnitude. The term v(dpl r dy) 
vanishes by (17), and the term d(kdTjdx)ldx is of order (S/l) 2 times 
d(kdT/dy)ldy. Finally, in <D, the term (du/8y)* is 0(wg/8 2 ), which is at 
least (Z/S) 2 times the other terms; these are therefore neglected. It 
follows that (5) becomes 



When body forces are negligible, an alternative form of (18) is obtained 
by adding u times (12) to (18). This yields 



d L dT , du 

~~-\ 

dy 



We note in passing that the terms pud(c p T)ldx and d(kdTjdy)jdy in 
(18) will be the same order of magnitude if 



8f = 0/-AL) = oiA 

\U p C p / \U ] 



or 8^ = 0(Pe-*), (20) 

where (Pe) is the Peclet number 

Pe = ^. 
" 



8 INTRODUCTION 1.3 

We deduce from (11) and (20) that 



so that the assumption of similar velocity and thermal boundary-layer 
thicknesses is equivalent to the assumption that the Prandtl number 
is not too different from unity. 

The equations (1), (7), (12), (17), and (19) are the boundary-layer 
equations for the five unknowns u, v, p, p, T, which must be solved 
subject to the following boundary conditions. At the wall y = we 
have u v = 0, since there is no slip or normal velocity at a fixed 
solid wall. Further there is a boundary condition on T at y = 0, which 
is usually that either T or dT/dy is prescribed there (or determinable 
from a subsidiary equation). At the edge of the boundary layer we 
may assume that u, T, p are functions only of x since variations in the 
quantities with y outside the boundary layer are small compared with 
those within the boundary layer. Further, the boundary layer tends 
into the mainstream asymptotically, so that the edge of the boundary 
layer is y -> oo, and the boundary conditions at the edge are u -> u^(x), 
T -> T^x), p -> PI(X}. From the practical point of view we may regard 
the edge of the boundary layer as the position where u\u v 0-99 say, 
though mathematically there is no edge as such. 

The solution of (17) is trivial, for if we let y -> oo in the momentum 
equation (12), then provided there are no body forces we find that 



so that the pressure is given in terms of the mainstream velocity and 
density. 

4. The laminar boundary layer on a curved wall 

The equations as derived above are valid only for flow past a plane 
wall. Physically one is often interested in flow past a curved surface, 
such as a cylindrical obstacle say. To consider this flow it is necessary 
to take a more general form of the Navier-Stokes equations than (1) 
to (6), by introducing general orthogonal coordinates x, y, with x 
measured parallel to and y normal to the wall. Since these coordinates 
are non-planar a considerable number of additional terms (involving 
the curvature K I of the wall) appear in the equations. The details are 
here omitted, and the reader referred to Goldstein's book ( 45), Modern 
Developments in Fluid Dynamics, vol. I, where the detailed analysis is 
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given for incompressible flow; the analysis for compressible flow follows 
the same lines. It is sufficient for us to remark that when the boundary- 
layer approximations are made exactly as in section 3, subject to the 
assumption that K I I is of order unity, then the equations reduce to the 
forms (1), (12), (18), and 

1 dp 

- _JL. = K U 2 . 

pey 

Accordingly the only difference between the equations for plane and 
curved walls lies here, since a normal pressure gradient is required to 
balance the centrifugal force. Nevertheless, the total change of pressure 
across the boundary layer is still small, since the boundary layer is 
thin, so the pressure may still be taken to be a function of x only. Thus 
the equations for plane and curved boundaries are identical. 

The practical significance of this fact is that in calculating laminar 
boundary-layer development, all curvature effects may be neglected, 
excepting in so far as these affect the flow outside and at the edge of 
the boundary layer. 

5. Conditions in the mainstream 

Certain relationships between quantities in the external stream will 
be constantly referred to in what follows. It is useful, therefore, to 
derive these relationships here briefly, so that they may be used when 
required, and further details will be found in any standard textbook 
on Gas Dynamics, such as Modern Developments in Fluid Dynamics, 
High Speed Flow, vol. I (Howarth, 1953). 

It is continually necessary to introduce the local velocity of sound, 
for which the usual symbol is a, and in a moving stream this is the 
velocity with which sound waves are propagated relative to the stream. 
It is given by 

(22) 



where the derivative is taken at constant entropy 8. Now in the ex- 
ternal stream, where viscosity and conductivity are negligible, the 
energy equation may be used to show that the entropy does not vary 
along a streamline, except in passing through a shock. Hence, outside 
the boundary layer, the derivative in (22) is taken along the streamlines. 
It also follows that if the entropy is uniform at infinity it will be the 
same everywhere (upstream of any shocks which may be present). 
Such a flow was formerly called isentropic, though Howarth (1953) uses 
the term homentropic and this term will be used here. The condition 



10 INTRODUCTION 1.5 

of homentropic flow may be used, along with the equation of state for 
a perfect gas, (7), to show that 

p oc P v (23) 

in the mainstream, where y is the ratio of the specific heats of the gas 
c p jc v) and is equal to 1-4 for a perfect gas. It is worthy of note that 
air behaves like a perfect gas over a considerable range of conditions. 
With the notation used in developing the boundary -layer theory earlier, 
(23) is written as 



^ (24) 

and it then follows from (22) that a l9 the local velocity of sound outside 
the boundary layer, is given by 

a\ = Yp/p, = Y MT^ (25) 

the second equality being deduced from (7). By making use of the fact 
that, for a perfect gas $ __ 

J( C p C,,, 

we may deduce from (25) that 

(y-l)c^. (26) 



Pi 
It may also be shown from (7) and (24) that 

Pi/Po = (TJT^r-V (27) 

and p/ Po = (^/TyyAy-D. (28) 

From the Navier-Stokes equations, (1) to (6), it is not difficult to 

show, when viscosity and conductivity are neglected and the external 

forces are assumed to be conservative and derived from a potential Q, 

that 



= constant 

along a streamline. This result may be obtained by adding together 
(5), u times (2) and v times (3). Upon neglecting external forces, it 
follows that at the edge of the boundary layer (where ?; 2 <: u 2 ) 

c pTi+z u l = constant, (29) 

which result agrees with that obtained by letting y -> oo in (19). Upon 
making use of (26) it follows that (29) may be written as 

T Hl = TI + ~ = 2l( 1 + ^ M\ } = constant, (30) 

Z,C p ( 4 ) 

where T H1 is the total temperature, or mainstream stagnation tem- 
perature, and M is the local Mach number u^di at the edge of the 
boundary layer. When M 1 is prescribed as a function of x, Tj follows 
at once from (30). If, however, the local velocity % at the edge of 
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the boundary layer is prescribed, the following relationships will be 
useful, namely 



6. Some standard boundary-layer characteristics 

In practice we are often not interested in the full details of the velocity 
and temperature profiles, but only in some overall measure of boundary- 
layer thickness or in some quantity as calculated at the wall, such as 
the heat transferred from the fluid to the wall or the friction al drag 
exerted by the fluid on the wall. We will now discuss some of the more 
interesting (and important) of these characteristics. 

Imagine inviscid flow past a body, the velocity of slip at the surface 
being u^x). In a real fluid the effect of viscosity is to slow down the 
fluid near to the surface, bringing the fluid at the surface to rest and 
causing the streamlines to be displaced outwards by an amount which 
may be calculated as follows. Consider a streamline which is a distance 
h(x) from the wall and, but for viscosity, would have been only 
h(x) ^(x) from the wall. Then the total mass of fluid flowing in unit 
time between y = and y = h is that which would be flowing between 
y = and y = h^ if the velocity were everywhere u = u^x). Thus 



or 

o 



p l u 8 t = J 



Thus streamlines far away from the wall have been displaced by an 
amount , 

81= f(l--^W (33) 

J \ Pi u J 
o 

which distance is known as the displacement thickness. 

In a similar way the difference between the actual rate at which 

h 

momentum is convected between y = and y = h, that is J pu* dy, 

o 
and the rate at which the actual momentum pu would be convected 

h 

by a flow with zero viscosity, that is J puu dy, may be compared with 
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inviscid flow momentum being convected in a region between y 
and y = S 2 . Thus h 





h 



(* I \ 

or 8 2 = -_.|i J dy. 

o 
As h -> oo this yields ^ 

8 I =fj5L(,_U, (34) 

J p 1 u l \ uj 





which is known as the momentum thickness. 

Two other measures of boundary -layer thickness which are of some 
interest are the kinetic -energy thickness, which is a measure of the 
defect of kinetic -energy flow within the boundary layer, and the en- 
thalpy thickness which is a measure of the defect of heat flow within 
the boundary layer. The kinetic-energy thickness is defined as 



, (35) 

PlMl u 



and the enthalpy thickness as 

(36) 





The skin friction at the wall, r wt is defined as the viscous force per 
unit area acting at the surface, and is given by 

(37) 

A position of zero skin friction is one for which dul'dy changes sign near 
the wall, and hence, since u = at y = 0, the local velocity must be 
of opposite sign on either side of such a point. It follows that when 
the boundary layer is of finite thickness the forward flow must 'separate ' 
from the wall in order to by-pass the backward flow. For this reason 
a position of zero skin friction is referred to as a position of boundary- 
layer separation. The serious effects of boundary -layer separation have 
already been referred to. 

The rate of heat transfer to unit area of wall, y = 0, is 

(38) 
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Often either Q w or T w is given, and it is desired to calculate the other. 
However, at high speeds a relationship between Q w and T w may be 
used, expressing the balance between heat transfer and radiation. 

7. The momentum integral equation 

Equation (12) expresses local conservation of momentum in the x 
direction at points within the boundary layer. When integrated across 
the boundary layer from y = to y = oo the resulting equation ex- 
presses conservation of momentum in the boundary layer as a whole. 
Upon neglecting the body force X, substituting for cJp/dx from (21), 
and integrating from y to y = S we obtain 
388 

r du 7 . r du , r <z% , 

J pU ^ dy+ J *%* - J *^*- 





8 

/ 

o 
8 



d'!/-r u ,. (39) 



Now J pv ^ dy == [puv]l - J u -j- ( P v) dy, 

' 





and by (1) this equals 

8 8 



\u.'-(pu)dyUi \---(pu}dy= (uu } )-^(pu)dy 

J CX J CX J CX 

00 

8 8 

^dx P u ( u ~- u ^ d y 



Hence (39) yields g S 

TW = I (pi u ip u )-d,jc dy ~dx J p u ( u - u i> 



upon formally letting 8 -> oo. This equation is called the momentum 
integral equation, and was first derived by Karman (1921) by more 
physical arguments. Afterwards, K. Pohlhausen (1921) derived the in- 
compressible form by the procedure adopted above. 

8. The kinetic- energy integral equation 

If equation (12) is first multiplied by the local velocity u, and then 
integrated from y = to y = 8, an alternative integral equation is 
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obtained. Upon doing this we find 

cj c\ 5> <j 

j u dy+ J^g) dy. 
oo 



The second terms on both sides of (41) are integrated by parts. Thus 

8 8 

f d( 8u\ 7 f ldu\ 2 , ,. nx 

utp,\dy= H dy, (42) 





and 



8 8 



f ^ 2 r ^ r i ' 

j dy * l J dx J 2 c 




8 



o 
Also 8 3 g 

r 3 /i 2^ ! r a / 2 2x^, ^i r ^ /^N 

I pu (j>u ) dy : = I pu~ f (u HI) dy~\-U'i - I pu dy. v"*^/ 

J 8x 2 J dx dx J 

00 

Upon substitution into (41) from (42) to (44) we find that 
8 88 

2 dx J dx J J \8y/ 



2 



or 

00 



f P^ /-, ^ 2 \ J I of /^\ 2 7 , ^1 f / 

J ^ii'-^^r J %) ^ +2Mi ^ J (p 



(45) 
Since, for a perfect gas, the equation of state is (7), it follows that 

p = p&T = Pl &T l9 
so that (45) becomes 

d{ 3 f pu L u*\j} 

T-\PI U I - L \ l i] d y\ 

dx ri X J p lWl \ ull j 

; 

-* <> 



o 
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This may be rewritten in terms of the energy thickness 3 3 and the 
enthalpy thickness 8 4 as 



3, () 



The incompressible form of this equation was first derived by Liebenson 
(1935). 

9. The thermal-energy integral equation 

This equation, as its name implies, is obtained by integrating the 
thermal energy equation (19) across the boundary layer. Now for any 
function F, the integral across the boundary layer of the expression 

{pu(dF/dx)+pv(dFldy)} becomes 

333 3 

[pu^-dy+ (pvdF = f P u^dy+[pvF}l- 

J &% J J & 

000 

Upon making use of ( 1 ) this becomes 

3 33 

J ("S^D dy = J l*""^-*' J 



dx 
o 

Choosing F = c p T lf = c^T+l^ 2 , (48) 

where T H is the total temperature, (19) becomes 



or ^ J p(c^ H -c p T m ) dy = ~*(^) ( 4 ) 



This equation expresses the fact that the rate of change of flux of total 
temperature must be balanced by the heat being transferred across the 
boundary. 

10. Incompressible flow 

For flow of liquids, or of gases at sufficiently low speeds, the density 
p will be sensibly constant unless the temperature differences are great. 
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When this approximation is made the boundary-layer equations take 
a much simpler form. The continuity equation (1) becomes 

^ + -^=0. (50) 

dx tiy 

Similarly the momentum equation (12) becomes 

du du 1 dp . d 2 u du, . d~u lf ... 

u -- \-v " ---- *L -\- v -. ^= u , 1-4-v -. (51) 

dx by p dx dy* dx dy 2 

Here we have neglected the body forces, and have taken the viscosity 
as constant, since JJL is a function of the temperature, assumed not to 
vary much. Finally, the thermal energy equation, (19), remains. If 
typical temperature changes are small like the square of typical 
velocities (in non-dimensional form, of course), then (19) is trivial and 
vanishes. When, however, we consider the case of 'warm walls', with 

I\T 

f - 0(V /a), 

where U/a is the relevant non-dimensional velocity, then (19) reduces to 



We notice that equations (50) and (51) are two equations for the velocity 
field (u, v), which are uncoupled from the equation (52) for T. The 
relevant boundary conditions are 

u v = at y = 0, 

u -> u^x) as y -> oo , 

T -> TI = constant as y -> oo, 

T = T w (x), prescribed, at y = 0. 

In a similar way the integral equations simplify. We find that 
(42) becomes j -, 



or = (8 1+2 8 a ) + ?. (53) 

u dx 



In the same way (47) becomes 

00 

P |-KS 3 ) = 2j^) 2 %) (54) 



and (49) becomes , 

(55) 
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11. Crocco's transformation 

Various attempts have been made to transform the laminar-boundary- 
layer equations into a form more amenable to calculation. These have 
proved most useful in restricted fields, though none has been found of 
universal value. We begin by considering a transformation due to 
Crocco (1939), in which x and u are taken as new independent variables, 
with the viscous stress T - fji(du/dij) and the enthalpy / c p T as 
dependent variables. 

We may write the continuity, momentum, and energy equations, (1), 
(12), and (18), as * 

(pu) + --(pv) = 9 

dx dy 

du , du dp , fir 
pU - - + pV = --\ , 

r dx r dy dx dy 

-, dl , dl dp , 1 dl dl\ , r 2 

and pu \-pv w A la - M . 

dx dy dx a dy\ dy) p 

These equations are then transformed to new independent variables, 
x, u, and after some algebra (for which the reader is referred to Crocco's 
original paper) they become 

dy LJL dp dr ,_ 



(57) 
dxc)r* f ~' { } 

and (l-~a)^T^+T^^+a\~viJipu C ^+iJiv^(-~+u\ = 0. (58) 
du du \du 2 ] dx dx\du ] 

Equations (57) and (58) are simultaneous equations for / and r in 
terms of u and x. After these have been integrated, it- is a straight- 
forward matter to calculate u 

. , Cfi(I)du 

11 (x u} = I - 

J T(X, u) ' 



after which v follows from (56). 

This form of the equations will be used in discussing the compressible 
laminar boundary layer on a flat plate in Chapter 7. 

12. von Mises's transformation 

A related transformation, in which x and = J pu dy are taken as 
independent variables, was given by von Mises (1927) for incom- 

853502 C 
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pressible flow. Upon transforming, the equations become 

du p l u l du l d/ du\ /enx 

u -- ^L_J _^ = u -{wu , (59) 

dx p ax difj\ di/j] 



, , ..i , /flnx 

and ^ te+ ^ = l V J + ww ' (60) 

By adding these equations together we obtain, after some algebra, 



ex 

The equations in this form have proved useful in considering the 
heat transfer through a laminar boundary layer in incompressible flow 
(Chapter 6), also the laminar boundary layer on a flat plate in com- 
pressible flow (Chapter 7). 



ANALYTIC SOLUTIONS FOR INCOMPRESSIBLE 

FLOW 

1 T has already been pointed out that in incompressible flow the momen- 
tum and energy equations become uncoupled, so that the momentum 
equation is independent of the temperature field. When this equation 
has been solved to yield the velocity distribution, the thermal energy 
equation, which is linear in temperature T, must be solved. In Chapters 

2 to 5 various aspects of the first half of this problem will be considered, 
and solution of the thermal energy equation discussed in Chapter 6. 

In this chapter we consider solutions of the incompressible momen- 
tum equation which reduce to the integration of one or more ordinary 
differential equations which can, in principle, be solved to any prescribed 
accuracy. Naturally such a situation only arises when the physical 
picture is extremely simple, such as in flow parallel to a flat plate, or 
when only a restricted domain is under consideration, such as the flow 
near to the forward stagnation point of a bluff body. Mathematically 
the solution will depend upon ordinary differential equations only if 
the velocity depends not upon x and y separately, but upon some 
suitable non-dimensional combination of these and of the basic fluid 
properties /z, p, v. The obvious choice is ?//S, where 8 is the boundary- 
layer thickness, and since this is typically of order (vx/U)* the non- 
dimensional parameter is y( U/vx)^ 9 where U is a representative velocity. 

In spite of the simplicity and obvious limitations of these solutions 
they are none the less useful as exact solutions of the boundary -layer 
equations and, further, may serve as a valuable starting-point for a 
numerical investigation of more complex situations. 

1. Flow parallel to a semi-infinite flat plate 

We consider firstly the case of flow parallel to a semi -infinite flat 
plate, the flow field being otherwise unbounded and the plate assumed 
to have zero thickness. Then the everywhere uniform flow will be 
disturbed only within the thin boundary layer near to the plate, and 
the velocity outside this layer will accordingly be uniform, say equal 
to u Q . We take axes x, parallel to the plate, and y normal to it, with 
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the origin at the leading edge. Then the boundary-layer continuity 
equation (50) and the momentum equation (51) become 

8X + ^ = }, (61) 

du . du 8 2 u 

u \-v = v~ 

dx dy di 

with boundary conditions u = v = at y = 0, u -> U Q as y ~> oo. These 
equations may be reduced to one ordinary differential equation by 
introducing a stream function 0, such that 

#0 ^A ,-. 

u = _r. v r (62) 

r dy dx 

and writing t// = (UQVX)*/^), f] - 

Upon substitution into (62) this yields 



(63) 



where dashes denote derivatives with respect to 77, and when these 
values are substituted into (61) one obtains, after some algebra, 

r+ff" = 0. (64) 

The boundary conditions become/ = f = at 77 = O,/' -> 2 as 77 -> oo. 
The first accurate solution of this equation was given by Blasius (1908), 
by a method which required a series expansion for small 77 and an 
asymptotic expansion for large 17, sufficient terms being derived to 
yield a region of overlap. A better method of integration was later 
given by Topfer (1912), and this has been used by Goldstein (1980) and 
Howarth (1938) to yield extremely accurate solutions. By expanding 
a solution of (64) about the origin it is found that 

oo 

/=* a" +V+, (65) 

where 

1 1 11 375 

S = 2t> a i = -5J' == 8T' ffi 3 = -n7' 

and the constant a is to be determined from the condition that /' -> 2 
as 77 -> oo. Now (65) may be expressed in the form 

/ = * a(a*7,)+ = **F(oL*i)) 9 (66) 
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where F(rj) is the solution when a 1. By numerically integrating 
(64), using starting values derived from the series (66) with a = 1, the 
function F(r]) may be obtained to any desired accuracy. Then the 
boundary condition at 77 = oo yields 

lim/'(7?) = of lim F'(a*r)) = u* lim F'(rj), 

TJ-+OJ 17-^00 7/->OC) 

so that at = {P"(ao)}-*. 

The value of a has by this means been determined as 

a = 1-32824, 

whereupon numerical integration of (64) yields the function f(rj). The 
velocity distribution through the boundary layer, which by (63) is 
proportional to /'(?}), is tabulated below. 

TABLE 1 



The Blasius velocity profile 



o 

0-1 

0-2 
0-3 
0-4 
0-5 
0-6 
0-7 
0-8 
0-9 
1-0 



uju 



0-0664 

0-1328 

0-1989 

0-2647 

0-3298 

0-3938 

0-4563 

0-5168 

0-5748 

0-6298 



M 
1-2 
1-3 
1-4 
1-5 
1-6 
1-7 
1-8 
1-9 
2-0 



u]u 

0-6813 
0-7290 
0-7725 
0-8115 
0-8460 
0-8761 
0-9018 
0-9233 
0-9411 
0-9555 



2-1 

2-2 

3 

4 

5 

6 

7 

2-8 

2-9 

3-0 



0-9670 
0-9759 
0-9827 
0-9878 
0-9915 
0-9942 
0-9962 
0-9975 
0-9984 
0-9990 



The skin friction at the wall is given by 



The displacement thickness is 

OO 00 

c> r /i u \ j i vx \i r /a 

8, = \l \dy=z\ Y (2 
1 J I uj J \uj J V 



o 
1-72077( 



u 



r/ 
' J 



and the momentum thickness is 
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and upon making use of (64) this may be easily shown to yield 



This result could have been very readily deduced from the momentum 
integral equation (53). When the external velocity u-^ is constant, so 
that % = 0, we see that S 2 is directly proportional to the integral of 
the skin friction r W9 given by (67). 

2. Flow near the stagnation-point of a cylinder 

Sufficiently near to the forward stagnation-point in two-dimensional 
flow past a cylinder the external velocity u may be written as 

U 1 = j8 1 a;/c, (68) 

where x is measured along the cylinder and y perpendicular to it. Then 
the continuity and momentum equations become 

du dv \ 

~dx^"dy = 

, (69) 

du , du 02 , - , B 2 u 
u -- l-v = f$x /c 2 -f-v 
&r by ri I -r d y 2 j 

with boundary conditions u v = at y = 0, u -> % as y - oo. We 
introduce a stream function i//, such that u, v are given by (62), and 
VK rn/e 



Then 



and upon substitution into (69) we find 

/'"+jT-/' 2 +i = o, 

the boundary conditions being 



This equation was derived by Blasius (1908), and improved numerical 
solutions have been obtained by various workers, including Heimenz 
(1911), Howarth (1934), and Tifford (1954). The first derivative/' is 
tabulated as /i in Table 2. The skin friction and the momentum and 
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displacement thicknesses are easily shown to be 



3, = 0-64790(J)*, 
\Pi/ 



S 2 - 0-29234^* 

\Pi/ 

It can be seen from (70) that 8 2 u/dx 2 is identically equal to zero, so that 
the solution is an exact solution of the complete Navier-Stokes equa- 
tions at points sufficiently close to the forward stagnation -point, where 
(68) holds. As an indication of the region over which this approxima- 
tion holds, we may remark that for flow past a circular cylinder the 
development of the boundary layer follows this solution over roughly 
the front quadrant. 

3. The Falkner-Skan similarity solutions 

It was shown by Falkner and Skan (1930) that the two solutions 
given above are special cases of a general class of solutions characterized 
by the mainstream velocity 

,,/ _ O f^l^m 
u l ~- Pm( x l c ) - 

Upon introducing a stream-function 



then the momentum equation (51) may be shown to reduce to 

/'"+i(m+l)//"+m(l-/' 2 ) = 0. (71) 

Solutions were given by Falkner and Skan for a range of values of m. 
Later Hartree (1937) studied (71) in greater detail, first making the 
linear transformation 



so that it becomes 

* 



with boundary conditions 

F = ~ = at 
a JL 

For positive values of m, (72) has a unique solution, but more than one 



= ~ = at 7 = 0, 4J 
a JL ctj. 
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solution exists for negative values of m. Hartree gave a family of 
solutions for positive values of m, and for negative values of m from 
zero decreasing to 0-0904. For this last value d^F/dY 2 becomes zero 
at Y = 0, so that 

= when y = for all x. 
dy J 

In other words the external velocity is such that the boundary layer 
is everywhere on the point of separation. For the numerical values of 
F(Y) for various values of m see Hartree (1937). 

A second family of solutions of (72), the so-called 'lower-branch' 
solutions, was given later by Stewartson (1954), again for negative 
values of m in the range 0-0904 < m < 0. These solutions all have 
the property that (d 2 FjdY 2 ) w < 0, so that 

< when 11 = for all x. 
chf ^ 

Accordingly these solutions can only correspond physically to flow in 
a laminar boundary layer beyond separation. These solutions, involving 
the possible singularity at separation (Goldstein, 1948; Stewartson, 
1958), will not be further discussed here. 

4. Series solutions from a stagnation-point 

We now generalize the solution given in section 2, which was valid 
at points sufficiently close to the stagnation-point in flow past a round- 
nosed cylinder. Consider the case of flow past a symmetrically placed 
cylinder, for which the mainstream velocity is assumed to take the 
form 

1 Z, 2w-+l > c 

It was first noted by Blasius (1908) that a solution of the momentum 
equation (49) could be obtained by assuming an expansion for the 
stream-function 

IQ \ I 

I, (74) 

where the functions F 2n+1 are seen to depend only upon the distance y, 
normal to the wall. It follows that the local velocity is 

(75) 



(76) 





and the skin friction at the wall is similarly derived from 
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Accordingly the skin friction is obtained as a simple power series in . 
Upon substituting from (75) into the momentum equation (51) and 
equating powers of f , a series of equations for F 2v +1 are obtained, and 
these may be solved in turn. Heimenz (1911) developed the work of 
Blasius and put F l and F 3 in universal forms, independent of the 
constants /? 2 +i> hut did not succeed in dealing similarly with F 5 . 

The major progress to date was made by Howarth (1 934) who showed 
that all the F 2n hl could be reduced to sums of universal functions. 
This is done as follows : 



'-MM 



O3 ^ J 04 
Pi Pi 



. (77) 



The equations and boundary conditions for the universal functions will 
be found in Howarth 's paper, together with tabulated values of the 
functions / 1? / 3 , # 5 , h 5 , and & 7 . The functions / 3 depends upon / 1} and 
h 5 in turn depends on both / t and / 3 . It follows that the accuracy 
obtainable in successive functions decreases rapidly, so that although 
Howarth 's values of / x are given to four figures, those of / 3 are given 
only to three figures, g 5 and h 5 to two figures, and & 7 to one figure. 
Howarth 's values were later improved, and further functions calculated, 
by Frossling (1940) and Ulrich (1949). More recently TifTord (1954), 
using high-speed computing apparatus, obtained values of all the uni- 
versal functions shown in (77) to higher accuracy than obtained hitherto. 
The universal functions being known, the first six terms in the series 
(75) for the velocity follow at once, as do the first six terms in the 
series (76) for the skin friction. Other boundary-layer characteristics, 
such as displacement and momentum thickness, are given as follows. 



0-0 
0-2 
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0-000000 
0-226612 


0-00000 
0-12511 
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TABLE 2 
Universal Junctions 

05 *5 07 M 

0-00000 0-00000 0-0000 0-0000 
0-10722 0-01413 0-0962 0-0173 


0-0000 
0-0016 


2.4 

09 ^9 
0-0000 0-0000 
0-0884 +0-0112 


0-4 


0-414456 


0-21288 


0-17781 +0-01170 


1563 +0-0030 


0-0044 


0-1413 -0-0079 


0-6 


0-566281 


0-26881 


0-21838 -0-00106 


1879 - 


0-0286 


0-0096 


0-1669 


0-0417 


0-8 


0-685937 


0-29972 


0-23665 0-01767 


1994 


0-0637 


0-0174 


0-1740 


0-0760 


1-0 


0-777865 


0-31250 


0-23990 0-03306 


1980 


0-0925 


0-0271 


0-1700 


0-1019 


1-2 


0-846671 


0-31327 


0-23415 0-04419 


1896 


0-1102 


0-0369 


0-1604 


0-1157 


1-4 


0-896809 


0-30702 


0-22390 0-04989 


1782 


0-1159 


0-0452 


0-1489 


0-1176 


1-6 


0-932348 


0-29746 


0-21226 0-05040 


1665 


0-1114 


0-0506 


0-1375 


0-1101 


1-8 


0-956834 


0-28712 


0-20115 0-04682 


1558 


0-0997 


0-0525 


0-1276 


0-0965 


2-0 


0-973217 


0-27755 


0-19158 0-04061 


1469 


0-0839 


0-0510 


0-1195 


0-0798 
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0-26950 
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0-0466 


0-1132 


0-0628 
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0-0507 


0-0402 


0-1087 


0-0470 
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0-25859 


0-17401 0-01891 
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0-0367 


0-0330 


0-1055 
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2-8 


0-997046 
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0-25003 
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250 


0-0003 
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0-1000 


0-0003 


4-6 


0-999997 


0-25001 


0-16668 0-00008 
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0-0001 


0-0003 


0-1000 


0-0001 


4-8 


0-999999 


0-25001 


0-16667 0-00003 0-1250 


0-0001 


0-0002 


0-1000 


0-0001 
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0-16667 0-00001 
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0-0000 
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0-1000 
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0-0000 


0-1000 


0-0000 


, 


*$ 


a 


99 011 Ml 


*j, 


'ill 


gii 


m n 


n' n 


0-0 


0-0000 


0-0000 


0-0000 0-0000 0-0000 


0-0000 


0-0000 


0-0000 


0-0000 


o-oooo 


0-2 


+ 0-0019 


0-0125 - 


0-0062 0-0824 +0-0074 


0-0041 


0-0165 


0-0237 


-0-0359 


0-0103 


04 


-0-0112 


0-0288 


0-0124 0-1299 -0-0145 


0-0368 


0-0371 


0-0514 


0-0721 


0-0206 


0-6 


0-0311 


0-0525 


0-0190 0-1511 0-0489 


0-0799 


0-0651 


0-0863 


0-1095 


0-0307 
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0-1330 
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The displacement thickness is 



so that 



l IP2TH- 

~ 



The constant coefficients in the bracket are easily obtained from a 
knowledge of the asymptotic forms of F 2n+1 (r^) as r\ tends to infinity. 
No similar simple expression for 8 2 has been obtained. 

Some of the more important functions are tabulated here. In Table 2 
are shown the first derivatives of all the universal functions in (77), 
for use in computing velocity profiles. Table 3 gives the values of the 
second derivatives at the wall, for computing the skin -friction. Table 4 
gives the asymptotic values at infinity, for computing the displacement 
thickness. 

TABLE 3 

Second derivatives of universal functions 

/i(0) = 1-23259 gfj(0) = 0-5399 g'^((}) = 0-5.1 00 

/J(0) = 0-72445 AJ(0) = 0-1520 ^(0) ^ 0-1323 

0(0) - 0-63470 A;J(0) = 0-0572 ^(0) =-- 0-0742 

7^(0) - 0-11918 yj(0) - 0-0607 ^(0) = 0-0806 

= 0-57920 gj(0) = 0-0308 ^(0) = 0-1164 

- 0-1829 Wi,(0) = -0-1796 

= 0-0076 wJ^O) = 0-0516 

TABLE 4 
Limiting forms of universal functions for large rj 



77-0-64790 
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IT? + 0-02847 


7i 9 = 0-1876 


^11 


= 0-1856 


^ ,7 _j_ 0-07002 


& 9 = -0-1153 


^11 


= -0-2479 


0-07390 


j 9 = 0-3567 


3 11 


= 0-3731 


17 + 0-0824 


g 9 = -0-1397 


0ii 


= 0-4424 


-0-1809 




tn\i 


= -0-6750 


0-0995 




Wn 


= 0-2075 



The functions already computed are sufficient to give six terms of 
the series in any problem. In general, however, the convergence will 



28 ANALYTIC SOLUTIONS FOR 2.4 

be such that more terms are required as separation is approached, and 
accordingly, since these could only be obtained at the cost of excessive 
labour, alternative methods must be used to continue the method to 
separation. These will be discussed in Chapter 3. 

A similar procedure can be adopted for flow past a cylinder placed 
non-symmetrically, when the mainstream velocity takes the form 



1 

and the stream-function is expanded in a series 



For details of this case the reader is referred to Howarth's paper 
(1934). 

5. Series solutions from a sharp leading edge 

We consider now the case of flow past a body which has a sharp 
leading edge, such as a flat plate placed parallel to the stream. The 
mainstream velocity is written in the form 



:, f - -, (78) 

o c 

and a solution sought by expanding the stream -function in a series 

o " ^ I 
__ i/AAl ' (79) 

the local velocity thus being 



o 
and the skin friction being derived from 



As for the case of expansion about a stagnation-point, substitution 
from (79) into the momentum equation (51) leads to a series of equations 
for F n which may be solved in turn. It appears to have been noted 
first by Tifford (1954) that the solution for the general case may be 
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expressed in terms of universal functions by writing 
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(82) 



Not all the universal functions given above have been calculated, but 
/o> /i? ^2> ^3> an d #4 have been tabulated by Howarth (1938) and g%, g^ 
and & 4 have been calculated but not tabulated by Tani (1949). Accord- 
ingly in the general case, only the first three terms can be obtained 
from the functions at present calculated. However, if the four functions 
</ 3 , h 3 , h 4 , and^ 4 could be calculated, the first five terms would be known 
for arbitrary mainstream. 

6. Gortler's modified series expansions 

Gortler (1957 a) has given an alternative series method which may 
be applied to cases with any shape of leading edge. He introduces new 
dimensionless independent variables 



X 

1 f 
= - u^x) dx, 



yu^x) 



and writes the stream -function as 



Substitution into the momentum equation shows that F(g, rj) must 
satisfy the equation 

8F 



where 



2u' 1 (x)ju 1 (x)dx 



It follows that the only explicit entry of the mainstream velocity is 
via the function 
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Gortler indicates that when the mainstream velocity takes the form 
(78), then ]8(f ) takes the form 



i 

where these /3 n are different from those of section 5, but can easily be 
deduced from them. A solution for F in terms of universal functions 
may be obtained by writing 



where 



= jB a / 2 +j8?/ u 



and similar expressions for F and F & . 

Gortler (19576) has calculated these universal functions, and pointed 
out that the series expansion appears to converge rather better than 
the usual Blasius series, discussed in section 5. For the case of flow 
past a round-nosed cylinder, the Gortler series does not converge more 
quickly than the Blasius series, discussed in section 4. For details see 
Gortler (1957 a). 

7. Meksyn's technique 

It is perhaps valuable at this point to make brief mention of an 
interesting technique developed by Meksyn (1950, 1956, 1960, 1961), 
and applied by him to the calculation of laminar boundary layers. 
The method was later extended by Merk (1959) and applied with some 
success to the calculation of heat transfer through a laminar boundary 
layer. A discussion of the method, with particular reference to the 
heat transfer problem, will be presented in Chapter 6, section 9. 

It may be remarked that the method is capable, in principle, of any 
required accuracy; in practice fair accuracy is attained with a fairly 
crude approximation, but considerable ingenuity is needed to improve 
on this. For full details reference may be made to the book by Meksyn 
(1961), New Methods in Laminar Boundary Layer Theory. 



NUMERICAL SOLUTIONS FOR 
INCOMPRESSIBLE FLOW 

APART from those .special cases, considered in Chapter 2, where the 
boundary-layer equations may be reduced to one or a number of 
ordinary differential equations, one is faced in general with the problem 
of solving non-linear partial differential equations, and this is an ex- 
ceedingly tricky numerical procedure. For this reason very few precise 
numerical solutions have been obtained. The importance of these few 
may perhaps best be illustrated by reference to the approximate method 
of Pohlhausen, to be considered in Chapter 4, section 1. For many 
years this method was regarded as being perfectly satisfactory for 
calculating the development of a laminar boundary layer in the presence 
of an arbitrary pressure gradient. Then, following the classical experi- 
mental work of Schubauer (1 935), it was found that when Pohlhausen 's 
method was applied to Schubauer 's experimentally observed pressure 
distribution, it failed to predict separation by a considerable margin, 
whereas separation clearly took place in the experiments. It was then 
apparent that the method was inadequate in regions of rising pressure. 
In consequence comparison with at least one accurate solution has since 
been regarded as essential before reliance could be placed upon any 
approximate method. 

It was to answer this need that Howarth (1938) considered a precise 
numerical solution for the flow in which the velocity at the edge of the 
boundary layer is given as 

/ ~A 

= j8 (lf). (83) 

This solution will be considered in section 1 of this chapter. 

Later Tani (1949) considered the three cases for external velocities 



with n 2, 4, and 8. The method used was precisely that of Howarth, 
and, although the accuracy is probably not quite so great (particularly 
for n = 8), there is no doubt that it is far more accurate than one could 
hope for by approximate methods. 
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In the meantime Hartree (1939 a, 6) had re-solved Howarth's prob- 
lem, using a differential analyser. Having thus confirmed the accuracy 
both of Howarth's calculations and of his own procedure, he then 
applied it to a numerically prescribed velocity distribution, derived 
from Schubauer's experimental results. This work gives a precise 
numerical solution for a flow containing a stagnation-point and separa- 
tion, whereas the solutions of Howarth and Tani did not have a stagna- 
tion-point. The only drawback to the solution was that arising from 
the uncertainties of twice numerically differentiating a tabulated func- 
tion, namely the experimental velocity distribution. 

More recently, additional solutions have been obtained for flows with 
stagnation -point and separation for analytically prescribed velocity 
distributions. Curie (1958 a) has considered the case 



by a semi-analytic method similar to but a little simpler than that used 
by Howarth above, and Terrill (1960) has used an electronic computer 
to obtain a solution for the case 

U-L = /3j sin f . 
We shall now consider some of the details of these various solutions. 

1 . The solution of Howarth 

As remarked earlier, Howarth (1938) considered the case of an ex- 
ternal velocity ,(*)=&(!-). 

By reference to Chapter 2, section 5, we see that this is the special 
case for which 

ft = -&, & = ft, = ]8 4 = ... = 0. 
Accordingly the stream function, 



= udy, 



may be expanded in a series 
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where 



and J^ 4 = 8 4 g 4 . 

Howarth derived the equations satisfied by these five functions (and 
by further functions in the expansion) and integrated several of them. 
The first three equations take the form 

fo+fof" = 0, 



and *?+/o*5- 

with boundary conditions 

/o = fi =-- h 2 = 0, fo=fi = h' 2 = 0, when 77 = 0, 
/o -> 2, /i ~> J, ^ -> 0, as 7^ -> oo. 

By integrating these equations in turn, together with the next four, 
and by roughly integrating a further two, the detailed velocity profile 
was obtained for a considerable distance from the leading edge. In 
particular the skin friction at the wall, which is given by (81) as 



was found to equal 

)...}. (84) 



The values of the second derivatives at the wall are perhaps best indi- 
cated by remarking that when 8 = 1 the series in brackets in (84) took 
the form 

1-328 1-021 0-069 0-056 0-037 0-027 0-021 0-017 0-015..., 
so that the series is converging rather slowly in this neighbourhood. 

It is clear, however, that subsequent terms form only a relatively 
small correction to this series. For this reason Howarth suggested repre- 
senting the terms after say the seventh in an approximate manner. 
Thus he approximated 



by = tf pFM+AtfWr,). (85) 
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In other words, it was assumed that the shapes of the subsequent terms 
F'niy) (n > 6), could be adequately represented by a universal form 
B(rj), the various multiples of f n being summed to give A(). By 
examining the functions F' 5 (rj), FQ(T]), and, as far as he had calculated 
them, -^7(7?) and F' 8 (rj), Howarth found that these were all multiples 
of the same functions, to a considerable accuracy. It remained, then, 
to determine A(), and this was done in each of the following two ways. 
Firstly, Howarth substituted from (85) into the momentum integral 
equation, and upon assuming the shape B(r)), this led to a first order 
non -linear ordinary differential equation for A(g), which took the form 

dA P+QA + RA* 

W = S+TA ' (86) 

where P, Q, R, S, T are polynomials which are determined by the 
known coefficients in (85). This equation was integrated by standard 
numerical techniques, and it turned out that 8+TA became zero and 
dA/dg infinite at a position where the skin friction had decreased almost, 
but not quite, to zero. This happening is clearly bound up with the 
question of the singularity at separation, discussed earlier. Howarth 
then completed this calculation of the skin friction by an alternative 
procedure, and concluded that separation occurred when = 0-120. 

The function A(g) was then calculated by the following alternative 
method. By considering the various boundary conditions at the wall 
(implied by the boundary layer equations), it is easily seen that the 
simplest one which is not automatically satisfied by (85) is that obtained 
by twice differentiating the boundary-layer momentum equation with 
respect to y, then putting y = 0. This yields 

'dtf 



By substituting into this relationship, an alternative equation for 
is obtained. By essentially this process, Howarth obtained almost 
identical values for A(), the separation position again being = 0-120. 
It should be noted that the two methods used to determine the func- 
tion A() are quite independent, and the close agreement between the 
two sets of results is valuable confirmation of their accuracy. This has 
been further confirmed, beyond all doubt, by the subsequent work of 
Hartree (1939 a) and Leigh (1955), who solved the same problem on 
a differential analyser and a modern high-speed computing machine, 
respectively. 
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2. The solutions of Tani 

In an attempt to provide further accurate solutions of the laminar 
boundary-layer equations, for cases in which the second streamwise 
derivative of the external velocity, d^ujdx 2 , is non-zero, Tani (1949) 
later considered the cases 

Wi(s)=j8 (l-f) fa = 2,4,8). 

As in Howarth's example, these are special cases of the family con- 
sidered in Chapter 2, section 5. When n = 2, so that 

*!(*)= 0o(l- a ), 

we have fa = & = ]8 4 = ... - 0, & - -]8 . 

Thus, by (82), we may write 



with /. 

^o = / o? 

* = o, 



and ^ 4 = 8 4 A: 4 . 

Tani calculated the first seven non-zero functions F , F&..., F 12) and 
then continued the solution by assuming a universal shape for the 
subsequent terms and determining the amplitude function A(g) from 
the momentum integral equation. Similarly, for the cases n = 4, 8, he 
calculated six and five functions respectively before the approximate 
continuation. In each case he obtained an equation of the form (86), 
finding that S-\-TA became zero and dAjdt; infinite just before separa- 
tion was reached. It was easy, however, to extrapolate the extremely 
small distance to separation, and the separation-points in these three 
cases were determined as 

| = 0-271 (n = 2), f = 0-462 (n = 4), = 0-640 (n = 8). 

3. The solutions of Curie 

As has already been mentioned briefly, Curie (1958 a) has developed 
a method for calculating solutions of the laminar boundary-layer equa- 
tions when the external velocity is that appropriate to flow past a 
cylinder symmetrically placed in a uniform stream. The external 
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velocity then takes the form of equation (73), with $, u, and (dujdy) u , 
given by (74), (75), and (76) respectively. In any chosen case the first 
six terms in these series expansions are known, since the relevant 
universal functions have been tabulated by Tifford, and in all the cases 
for which calculations have been made the series seem to converge well 
until just upstream of separation. For example, in calculating the 
solution for the case 



Curie estimated by reference to approximate solutions that separation 
would occur in the vicinity of 0-66, and showed that at this position 
the series (76) for the skin friction takes the form 

^} oc 0-8135 0-8331 + 0-0896 0-0033 0-0073 0-0064.... 



!3u\ 
W 



We note that the fourth to sixth terms are of order 1 per cent, of the 
first two, so that even if the subsequent (unknown) terms continue to 
decrease rather slowly they will presumably be only a relatively small 
correction to the first six. 

Following the idea used by Howarth (1938), Curie assumed that the 
shape of the subsequent functions F 2n+1 (n ^ 6) could be approxi- 
mated by F u . The accuracy of this approximation can only be tested 
by comparing say F' 7 (r))/F 7 (ao), F' 9 (-r])IF 9 (ao), F^rf/F^ao), to see whether 
there is any indication of a tendency towards universality in shape as 
n increases. In the examples he calculated, this tendency appeared to 
be present to a certain degree. On the basis of this approximation we 
write the velocity as 



(88) 
Pi 

and the skin friction as 

0), (89) 



where A() is to be determined. To do this, Curie suggested using the 
second of the two methods used by Howarth, that is by satisfying the 
first boundary condition at the wall which is not identically satisfied 
by a form (88). This, as in Howarth's case, is 

U\ 

(90) 
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From (88), by four-fold differentiation with respect to y = (vc//y*7?, 
we find A 4 

)- (91) 



Noting that, in (89) and (91), A(g) is multiplied only by constants, we 
can easily eliminate A(g) to give (d*u/ty*) w in terms of (dujdy) w , and 
this yields 



where 6 2)l+1 - F& +1 (0) - - gU FJ n+1 (0). 

Upon substituting from (92) into (90), an equation for (du/dy) w results. 
This equation may be put into non-dimensional form by writing 

T = (^\^(\ 

4 

and defining P() 2 ^zn+i P ?i ' +1 j 

o 

whence it becomes T - = P-4- =~-~ 2 1 . 



Upon formally integrating with respect to , this yields 

I 



f 
J 



f 4 , 

where Q = 2 f P ^ = 



The method of solution is therefore as follows. Given the external 
velocity u^x), the coefficients /3 2/( +i are known, so that the first six 
terms in the series (89) for T, the function Q(), and the constant 
2-Fj 1 (0)/jF T i 1 (0) in (93) are easily obtained. The series may be used to 
obtain the values of T when is small, and the solution is continued 
by numerical integration. A simple, and convenient, procedure for 

I 
doing this is to replace the integral [ T d by its Simpson's rule 

6 

equivalent, whereupon if T($) and T(-\-h) are known, T(-\-2h) is 
given by solution of a quadratic equation. This method appears to 
work particularly well, provided the constant 2JF T J 1 (0)/^ 7 J 1 (0) is not so 
great that the error in Simpson's rule is multiplied by too large a factor. 
When T(g) is thus determined, the details of the velocity profile easily 
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follow. The function A() is given from (89), and then u from (88). 
An a posteriori check on the accuracy of the solution may be made, if 
desired, by using the momentum integral equation. 

Curie used this method to calculate an accurate solution for the case 
(87). The values of T(g) at f = 0-250, 0-275, 0-300 were given by the 
known terms of the series to six significant figures. With these starting 
values, the equation (93) was integrated with steps of A = 0-025 and 
Af = 0-050, Richardson's ^-extrapolation being used to obtain im- 
proved values. Using smaller intervals as separation was approached, 
the integration proceeded without difficulty, and separation was pre- 
dicted at = 0-655 r 

Curie also gave solutions for two cases of the family of velocity 
distributions 



for which the value of the constant ^ii(O) was zero. He showed that 

JFJJO) = 1-02966 4-58707a 107-422a 2 
is zero when a = 0-07885 and a = 0-12156, 
and that in each case (87) takes the form 



Solutions were obtained, therefore, without numerical integration, and 
separation found at 

= 0-6647 when a = 0-07885, 
and = 0-6245 when a = 0-12156. 

4. The solutions of Hartree, Leigh, and Terrill 

The solutions discussed above were all obtained by semi-analytic 
means. We now turn to a number of precise numerical solutions ob- 
tained by predominantly non-analytic means, either by numerical 
analysis or by means of an automatic computing machine. We begin 
by considering a solution by Hartree (1939 a) for the case of a linearly 
decreasing velocity, which had earlier been considered by Howarth. 
The main idea of Hartree 's approach was that of approximating the 
partial differential equation (51) by an equivalent ordinary differential 
equation, derivatives with respect to x being replaced by finite differ- 
ences and the integration being carried out with respect to y, that is 
across the boundary layer. In other words the distribution of velocity 
at one station is used to calculate the distribution at a downstream 
station. 



3.4 INCOMPRESSIBLE FLOW 39 

The integration was started at 8| 0-40, Howarth's series having 
been used to give the velocity profile at that position. The profile at 
8f = 0-80 was derived in one and two steps, after which Richardson's 
A 2 -extrapolation formula was used and the results smoothed. Similar 
corrections for the finite size of the steps of integration were applied 
at 8 = 0-88, 0-94, 0-956, and 0-958. Now according to Howarth (1938), 
separation occurs when 8f = 0-960, and Hartree encountered con- 
siderable difficulties as this position was approached. These difficulties 
pointed very strongly to the existence of a singularity at the position 
of separation, and indeed it was because of them that Goldstein's 
investigations into flow near to a position of boundary-layer separation 
were carried out. Hartree determined the position of separation by 
two independent methods, and concluded that it is close to 8 = 0-959, 
in excellent agreement with Howarth's result. 

This same problem was later considered with even greater care by 
Leigh (1955), using the E.D.S.A.C. machine of the University of Cam- 
bridge. Various special tricks were successfully used to get over some 
of the difficulties which Hartree had met, and it was established that 
separation took place where 8 = 0-95854. The main object of Leigh's 
work was to improve on the accuracy of the earlier solutions, by making 
use of the greater storage capacity of the machine, so as to enable more 
information to be obtained regarding flow in the vicinity of separation. 
For details of the techniques see Leigh (1955). 

These precise numerical solutions are extremely valuable in that 
they confirm (by independent methods) the considerable accuracy of 
Howarth's procedure, as well as throwing some light on the singularity 
at separation. 

We now comment on a further solution obtained by Hartree (19396), 
on a differential analyser, in which the external velocity (or pressure) 
distribution was that found experimentally by Schubauer (1935) in 
some experiments with an elliptic cylinder. The procedure used was 
essentially that developed for the case of a linearly retarded velocity 
save that in this case the integration was done mechanically rather than 
numerically. One practical difficulty arose, however, in that the as- 
sumed external velocity, being given from experimental data, had to 
be smoothed and later differentiated numerically. There is a certain 
latitude of interpretation in doing this, which may not normally be 
critical but appeared to be so in this particular case. For example, 
Hartree smoothed and differentiated the values of %(#) given by 
Schubauer, and then calculated the development of the boundary layer 
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by numerical integration. With the minor axis of the ellipse c as the 
representative length, Schubauer's experiments indicated that the 
boundary layer separated at roughly x/c = g = 2-0, whereas in Har- 
tree's calculations no separation was predicted. It was found that an 
extremely small change in the assumed values of u-^x) in the vicinity 
of = 2-0 sufficed to cause a theoretical prediction of separation at 
the correct position. Accordingly we may regard Hartree's work as 
providing an accurate numerical solution for a case with a numerically 
prescribed external velocity, but it should be treated with some caution, 
in view of the uncertainties discussed above. 

A final, and important, solution on an automatic computing machine 
is that due to Terrill (1960), who considered the case of an external 



The method used involved the application of Gortler's transformation 
(1957 a) followed by the technique so successfully used by Hartree 
above, in which derivatives with respect to x are replaced by finite 
differences. The size of step used was rapidly decreased as separation 
was approached, and this was predicted at 

f = 1-823, 

Terrill also carried out some calculations which throw some light on 
the nature of the singularity at the separation position. For details 
see Terrill (1960). 



PRACTICAL METHODS OF CALCULATION FOR 
INCOMPRESSIBLE FLOW 

IN practice one may not be interested in a solution of the boundary- 
layer equations as accurate as would be obtained by using the methods 
discussed in the preceding chapters, and certainly one would not 
normally wish to spend the considerable time required to obtain such 
solutions. Accordingly it may be sufficient to use methods which are 
quicker to apply, but yield less accurate results. Some of these will 
now be considered. 

1. Pohlhausen's method 

The method developed by K. Pohlhausen (1921) is not now used much 
in its original form, having been improved upon and superseded by 
other methods. It is included here, however, because it was historically 
the first general method to be developed, and because the idea used by 
Pohlhausen has found so many applications in various branches of fluid 
dynamics. 

In this method the boundary-layer equations are not solved every- 
where, but are satisfied at the wall, at the edge of the boundary layer, 
and on an average by satisfying the momentum integral equation. The 
velocity u is assumed to take the form 



where 8 may usefully be regarded at this stage as the effective boundary- 
layer thickness, but in the final formulae is no more than a convenient 
parameter. Then the skin-friction r w , the displacement thickness 3 19 
and the momentum thickness 8 2 are equal to 



Ul S J 
I 



(94) 
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Now the boundary conditions to be satisfied by the true local velocity 
u are 



and uu v , , - -, ... ~> as y ~> oo. (96) 

The idea used by Pohlhausen is to choose /(T?) to be a polynomial satis- 
fying some of the boundary conditions, (95) and (96), leaving 8 to be 
determined so that the resulting approximate u satisfies the momentum 
integral equation. If we write 

A = ~?4, (97) 

V 

it can easily be shown that the above conditions lead to the equation 

"', (98) 



where g(A.) and h(A) are universal functions which depend only upon 
which of the boundary conditions (94) and (95) are satisfied, and are 
independent of the data of the problem. 

There are two criticisms which can immediately be levelled at Pohl- 
hausen 's method as presented above. The first is the arbitrary choice 
as to which boundary conditions shall be satisfied and the fact that 
the outer conditions are satisfied at y 8. Pohlhausen chose a quartic 
for u which satisfied the first two conditions of (95) at the wall, and 
the first three conditions of (96) at the edge of the boundary layer. 
This leads to a method which, when applied to the cases for which 
accurate solutions are known, yields accurate results in regions where 
the pressure is decreasing but becomes less accurate as separation is 
approached, so that the predicted distances to separation are typically 
30 per cent, too high. Curie (1957) suggested that better results would 
be obtained in the region near to separation if an additional boundary 
condition at the wall were satisfied by choosing a quintic profile. Such 
a procedure leads to improved results throughout the region down- 
stream of the pressure minimum, so that a typical predicted separation 
distance is say 6 per cent, too high only. However, the method then 
breaks down entirely near to a stagnation-point. Earlier Schlichting 
and Ulrich (1940), in connexion with some work on boundary -layer 
stability, had chosen to represent u by a sextic, satisfying still one 
condition more at the edge of the boundary layer. This procedure also 
leads to better results near to separation, say 15 per cent, error in the 
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separation distance, but again leads to a breakdown of the method near 
to a stagnation -point. 

A second criticism of the Pohlhausen approach is that the second 
derivative of u v that is u'[, appears explicitly in the formulation. Now 
if % has been obtained from experimental data it will in many cases 
be extremely difficult to obtain even rough estimates of u'{. Fortunately 
the method can be recast, by a procedure due to Holstein and Bohlen 
(1940), so that u{ does not appear, and it is necessary only to derive u[ 
by numerical differentiation. Holstein and Bohlen introduce a para- 
meter A, analogous to A, denned by 

S2 /S \2 

A = X = ? A. (99) 



V 



The momentum integral equation is written in the form 

r w S 2 u^ S 2 dS 2 Sf , L 8 t \ 

. - . . -- _ j H, ^,_j_ 

JU/M! v dx v \ o 2 / 



i A 



which we write as I = i (8|)+A(T+2), (100) 



where I = -H!_? and H = ^. 

Thus 



= L. (102) 

Now from (94) and (99) we know r w , 8 1? S 2 , and A as functions of A. 
Accordingly by (101) and (102) we can deduce I, H, and L as functions 
of A, and (102) reduces to a first-order equation for A, as can be seen 
by writing it in the form 

(103) 



The functions gr(A), ^(A), for use with the method in its original form, 
may be found in Pohlhausen 's original paper or in either of the books 
by Goldstein (1938) and Schlichting (1955), whilst the functions H(X), 
Z(A), L(X), can be found in Schlichting 's book. We may mention here 
that Tani (1941) indicated that the function L( A) is approximately linear, 
so that (103) integrates analytically. This considerable simplification 
has been used by several other authors, notably Thwaites (1949). 
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2. Timman's method 

Timman (1949) attempted to obtain an improved procedure for 
calculating laminar boundary -layer development by assuming a velocity 
profile of the form 

00 
U j,, . _ r ./ , 9 v 7 ^2/r , -, 9 v / -, n A v 

~ J(T]) I e~ 7 i (ct, -\-Crj*... ) ar) e ' (O-j-aiy ...), (1U4) 

Ui J 

1 17 

where ?? = y/8. This form of profile is an improvement over the poly- 
nomials assumed in Section 1 in that uju^ tends to unity asymptotically. 
Accordingly only the boundary conditions at the wall need to be further 
considered, and these serve to determine a, b, c, d. 
From (95), substituting from (104), these yield in turn 

1 b |V7r(a+fc...) = 



(105) 



2(b~d) = -!^L= -A 
v 

2(c-a) = 

11 S 2 (fin In' 1 Aft 

12(2d-b) = ^a[f + a^-lf 
v [dx \% o dx i 

Initially Timman considers only the first three conditions of (105), and 
these determine a, b, c as functions of A when d is set equal to zero. 
The method then proceeds much as in Pohlhausen's method. From 
(104) we can obtain Sj/8, S 2 /S, and r^8/jLc% as functions of A. and also A 
as a function of A. Thus we obtain I, H, and L, defined in (101) and 
(102), as functions of A, with the momentum integral equation reducing 
to the form (103). 

Timman found that this solution yielded accurate representations of 
the flow in a Blasius boundary layer and of the flow near to a stagna- 
tion point. In regions of retarded flow, however, the method was much 
less satisfactory, as application to the case % = u {l (%/c)} showed. 
Accordingly he suggested that in regions of adverse pressure gradient 
the condition d should be replaced by 

2d-b = 0. 

This condition was chosen to ensure that the complicated fourth condi- 
tion in (105) should be satisfied at the separation position, and inci- 
dentally led to values of a, 6, c, d, which were continuous at the join 
A = 0. This procedure leads to results which agree quite well with the 
accurate calculated results discussed in Chapter 3. The values of the 
relevant functions, required for integrating (103), are to be found in 
Timman's paper. 
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3. The methods of Howarth and Walz 

In the methods of Pohlhausen and Timman the velocity profile at 
any station is chosen as one of a singly-infinite family of velocity 
profiles whose shape is characterized in either case by the parameter A 
and whose scale is characterized by 8 = (vAfu'^. The choice as to 
which profile is to apply at each station is made by satisfying the 
momentum integral equation. Similarly the accurate solution of the 
boundary-layer equations for the case % = i%(l (x/c)}, derived by 
Howarth (1938) and discussed in Chapter 3, section 1, yields a singly- 
infinite family of velocity profiles, whose shape may be characterized 
by the parameter A, and with scale represented by say the momentum 
thickness 8 2 ( v ^l u i)^> A method on these lines was given by Howarth 
(1938), by regarding a general velocity distribution as being replaced 
by a series of linear portions. His general result can easily be reduced 
to the form (103), the functions //(A), Z(A), L(X) being obtained from 
Howarth 's tabulated results. These values, of course, will all be for 
negative A, since the velocity is everywhere retarded, but by considering 
the associated accelerated flow u = %{l + (#/c)} it would be possible, 
as Walz (1941) pointed out, to extend the functions to cover a range 
extending from the stagnation point to separation. 

Walz (1941) also suggested a method in which the velocity profiles 
are chosen to be those of the Falkner-Skan similarity solutions. As 
distinct from Howarth 's method, where each value of x corresponds to 
unique values of H, I, L, and A, each complete similarity solution for 
a given value of the velocity exponent yields single values of H, I, 
L, A for all x. However, the complete family of similarity solutions 
yields values for H, I, L, which are continuous functions of A. 

4. The method of Thwaites 

It was pointed out by Thwaites (1949) that if one wishes to calculate 
only the boundary -layer thicknesses and the skin-friction distribution 
it is not necessary to introduce explicit assumptions concerning detailed 
velocity profiles, as the methods previously discussed all do. Provided 
suitable correlations are defined between the overall boundary-layer 
characteristics H, I, L, and the shape parameter A, one can easily obtain 
A as a function of x by numerical integration of (103), after which S 2 
follows from (99) and then 8 X and r w from (101). 

Each accurate solution of the laminar boundary -layer equations will 
correspond to a different set of functions H (A), /(A), L(X). Similarly any 
approximate method using a single-parameter family of profiles will 
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correspond to a definite set of functions. It follows from this that an 
approximate method will yield an accurate answer for a given problem 
if and only if the values of H, I, and L for the accurate and approximate 
solutions agree closely enough. Now the idea used by Thwaites was to 
examine and compare the functions // , I, and L, corresponding to all 
known exact and approximate solutions of the laminar boundary-layer 
equations. By considering the accuracy to which these various solutions 
are known, and how typical each is of a practical boundary layer, 
Thwaites defined what are virtually the optimum values of the functions 
for general use. 

Thwaites found that as regards H(X) and l(X) the solutions lay very 
closely together for positive A, that is for regions of increasing velocity 
upstream of the pressure minimum in other words. Downstream of the 
pressure minimum, for negative A, the curves deviated considerably, 
so that values of A corresponding to boundary-layer separation, that 
is values for which Z(A) = 0, were obtained ranging from A = 0-068 
to A = 0-157. Fortunately the two most important solutions of those 
available to him, namely Howarth's solution and Hartree's analysis of 
Schubauer's experimental results, lay quite close to each other, and 
Thwaites 's choices of H and I lay close to both of these as well as being 
a rough average of the approximate solutions. 

The values of L(X), derived from (102), for the various solutions were 
much closer than those of either H(X) or l(X), both for positive and 
negative A, and in addition the function L(X) appeared to be linear to 
a high degree of accuracy. Thwaites accordingly suggested defining 

L ( A ) to be L(X) - 0-45-6A, (106) 

the coefficients being chosen to give the best overall agreement with 
the available solutions. With L(\) given by (106) it is possible to 
integrate (103) in terms of a single quadrature, to yield 

(107 ) 

(108) 



v 
Having thus obtained 8 2 , (101) then yields 



and 
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The idea of expressing L(X) as a linear function had been suggested 
earlier by both Walz and Tani, but Thwaites's coefficients are, if any- 
thing, more acceptable, since they are not based upon a single solution, 
and yield the computational simplification of requiring only integral 
powers of %. 

Some improvements to Thwaites's functions were suggested by Curie 
and Skan (1957). Their idea may be illustrated thus. Thwaites deter- 
mined his H (A) by considering the ratios of exact Sj to exact 8 2 . Since 
the predicted 8 X is obtained by multiplication of H and the approximate 
S 2 , Curie and Skan suggested that more acceptable values for H(X) 
would be derived from the ratios of exact 8 X to approximate S 2 . This 
argument would be complicated further if one accounted for the fact 
that an approximate value of A, derived from (107) and (108), is used 
in practice. In this way, and by similar arguments applied to l(X), Curie 
and Skan suggested improvements to Thwaites's functions in the region 
near separation where the functions for individual solutions differ most. 
These modified functions are given in Table 5. 

TABLE 5 
Universal functions for Thwaites's method 

m l(m) H(m} m l(m) H(m) 

0-25 0-500 2-00 0-040 0-153 2-81 
-0-20 0-463 2-07 0-048 0-138 2-87 

0-14 0-404 2-18 0-056 0-122 2-94 

0-12 0-382 2-23 0-060 0-113 2-99 

0-10 0-359 2-28 0-064 0-104 3-04 
-0-080 0-333 2-34 0-068 0-095 3-09 

0-064 0-313 2-39 0-072 0-085 3-15 
-0-048 0-291 2-44 0-076 0-072 3-22 

0-032 0-268 2-49 0-080 0-056 3-30 

0-016 0-244 2-55 0-084 0-038 3-39 
0-220 2-61 0-086 0-027 3-44 

+ 0-016 0-195 2-67 0-088 0-015 3-49 

0-032 0-168 2-75 0-090 3-55 

The modifications of Curie and Skan also made use of the accurate 
solutions given by Tani (1949), which had not been available to 
Thwaites. Subsequently Curie (1958 a) has given the additional accurate 
solution discussed in Chapter 3, section 3. The greater number of 
accurate solutions now available makes it much easier to avoid the 
possibility of being unduly biased by the few accurate solutions avail- 
able to Thwaites. Curie and Skan (unpublished) have examined the 
functions H and I which would lead to exact prediction of 8 X and r w 
for each of the accurate solutions discussed in Chapter 3, and have 
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indicated that with the exception of the modifications near separation, 
the original values given by Thwaites probabJy reflect these as accurately 
as any single set of functions can do. 

5. The method of Stratford 

The method developed by Stratford (1954) is based upon the idea of 
dividing the boundary layer into outer and inner portions, for each of 
which a solution is obtained which joins smoothly onto the other. This 
idea was first used by Karman and Millikan (1934). 

Stratford begins by considering the outer part of the boundary layer, 
where the flow is nearly inviscid and so, by Bernoulli's equation, the 
total head is almost constant along streamlines. Thus if the pressure 
is constant between x = and x = X Q , with a pressure rise downstream 
of x = XQ, then we have 

(*P a W = (iP*'W-(?-3>o)+Aff, (109) 

where the stream-function 1/1 is given by 



= Judy, 



(110) 



and AI/ is the small change in total head due to viscosity. Now the 
shape of the outer part of the boundary-layer velocity profile is not 
greatly changed by the pressure rise, so the viscous forces are given 
approximately as though there were no pressure rise. Hence 

AJ/ = (^)^-(MW, (111) 

where U B is the Blasius solution for a boundary layer with no pressure 
gradient. From (111) and (109), remembering that u = U D when x = # , 

wehave ' 



Stratford also indicates how this result may be deduced algebraically 
by writing the boundary-layer equation (51) in the form 

(p+J/>^ 2 ) = /z 7 -, (113) 

where djds denotes differentiation along a streamline, and considering 
the Taylor expansion 

The first term on the right-hand side is independent of the pressure 
rise, since the effects have diffused outwards only downstream of x = x . 
Similarly the coefficient of (x X Q ) is equal, by (113), 
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and so again is independent of the pressure rise. Hence, if terms of 
order (x x Q ) 2 are neglected, the left-hand side of (114) is independent 
of the pressure rise, so that 



which is (112) rewritten. 

By similar arguments, or by differentiation with respect to 
can also deduce from (112) or (115) that 



(116) 

, H ,.__, 

and =[ -- (117) 



in the outer part of the boundary layer. 

Sufficiently near to the wall, in contrast, viscous forces become very 
important, so that the resultant viscous force on an element of fluid 
at the wall must exactly balance the pressure -gradient force. Thus, 
by setting y in the boundary -layer equation (51) we find 

\ dp 

AX y ' 

By further considering the boundary conditions at the wall, Stratford 
indicates that in general /P3 \ 

~ =0, (119) 



but ^0. (120) 



By satisfying (118) and (119), Stratford suggests that at separation the 
velocity distribution near to the wall may be approximately obtained 
by writing , , 

u = ?y*+W, (121) 

2ju- ax 

where h is to be determined. 

With an inner profile given by (121), and an outer profile given by 
(115) to (117), it is postulated that there should be continuity of i/r, u, 
dujdy, d 2 u/dy 2 at the position of the join, ^ = ^ jy u = Uj. These four 
conditions serve to determine the values of Up ^, h^ and x at the 
position of separation. In general the solution of the four equations 
is difficult, since such functions as U B , du B \dy, and d 2 u B /dy* cannot be 
expressed in simple analytic terms. A considerable simplification is 
effected, however, if the assumption is made that the inner layer is 
sufficiently thin that the Blasius velocity profile u B (x,\l>) is sensibly 
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linear throughout /f ^ ^. This condition, which restricts the applica- 
tion of the method in a manner which we shall discuss later, leads to 
the simple result 



(122) 



Ixjji 

By assuming that these values hold at the join, Stratford derives the 
simple result / ,^ . 2 

~2J~^) =^0-0065 at separation, (123) 



where C p is the pressure coefficient, denned as 

n P~Po 



The condition on the inner layer is such that the method should strictly 
only apply when C p is less than about 0-04, with an absolute upper 
limit of about Oil. It has been found, however, that the method yields 
useful results over a larger range of pressure coefficients. 

Stratford also gave an asymptotically exact solution for the case when 
the adverse pressure gradient downstream of x = X Q is constant but 
asymptotically large, and the distance x s x to separation tends asymp- 
totically to zero, as does the thickness of the inner layer. For details 
of this solution reference may be made to Stratford's paper. The 
principal result is that 

x^cJ^Y = 0-0076 at separation. (124) 

\dx / 

The numerical coefficient is not very different from that obtained by 
the more physical approach discussed earlier. Stratford suggested that 
in general (124) should lead to good predictions of separation. 

Stratford also considered a second, very special, type of flow, in which 
the pressure gradient is zero upstream of x = X Q and is such that the 
skin-friction is identically zero downstream of x = X Q . He suggests that 
a more appropriate form than (121) for the inner layer would be 

1 dp 



An analysis along the same lines as that for the earlier case leads to a 
similar result to (123) or (124). The numerical coefficient becomes 
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0-0049, and a rigorous mathematical analysis for an asymptotically 
large adverse pressure gradient of the appropriate form leads to a value 
0-0059. By carefully interpolating between these special cases and 
empirically adding further terms, Stratford also derived a considerably 
lengthier formula than those quoted above, which should yield exceed- 
ingly accurate predictions of separation. Later, however, Curie and 
Skan (1957) indicated that an equally accurate, but considerably simpler, 
empirical improvement of (124) could be obtained by replacing the 
numerical coefficient by 0-0104; separation is then predicted to occur 
where 

= 0-0104. 
dx 

It should be pointed out that the method is based upon a considera- 
tion of the downstream deformation of a known Blasius velocity profile. 
Accordingly the method is not, strictly speaking, applicable to flows 
containing a forward stagnation-point. It was pointed out by Stratford, 
however, that in such a flow the velocity profile at the pressure mini- 
mum will be reasonably similar to a Blasius profile of suitably chosen 
scale. Accordingly the method may be, applied downstream of the 
pressure minimum if we replace x by x X L say, where X L is chosen so 
that the true momentum thickness at the pressure minimum x = x m 
is equal to the momentum thickness of the equivalent Blasius layer 
at x = x m x L . Since the momentum thickness S 2 * s usually given 
accurately by Thwaites's formula 

Sf = 0-45vtef 

o 
this condition is easily seen to yield 

m 

(125) 

where u m is the value of % at the pressure minimum. 

The method as developed by Stratford was designed merely for pre- 
dicting separation. No estimate of the boundary -layer thicknesses or 
of the detailed distribution of skin-friction between the leading edge 
and separation is given. This restriction has been partly removed by 
the method which will now be discussed. 

6. The method of Curie 

This method is a generalization of Stratford's method, and leads to 
a simple method for predicting the detailed distribution of skin- 
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friction. The inner layer is assumed to be of the form 

u = r 4-1^ 2 -4-a n (126) 

2 dx 

where r w is the skin-friction jji(c)ul8y) w , and a(x) and the constant n are 
to be determined. This profile generalizes Stratford's form in two ways, 
namely that the skin-friction is not assumed to be zero and that the 
value of n is not fixed a priori. The outer velocity profile is given, as 
in Stratford's method, by (112), (116), (117), and (122). The four join- 
ing conditions, continuity of /f, u, dufdy, d 2 u/dy 2 , are sufficient to deter- 
mine r w9 a, y^ i/j,j as functions of x, with n as an arbitrary parameter. 
After the relevant algebra has been performed the following results are 
obtained: dp^ _n-l 

and 

^n / . 

L dx 



(128) 

where T#(X) is the skin-friction at station a: of a Blasius boundary layer. 
The position of separation, where r w = 0, is obtained from (128) as 

_ 0-00405( m -l)'(n+3) 

n*(n+l)(n-2) ' ( } 

so that the predicted separation position depends upon the choice of 
the parameter n. When n = 6 the right-hand side becomes 0-0049, and 
when n = 4 it becomes 0-0065, the values obtained by Stratford. In 
view of the fact that, as Curie and Skan have shown, the best overall 
predictions of separation are obtained when the right-hand side of 
(129) is set equal to 0-0104, it would seem most natural to choose n 
accordingly, so that 



With this value of n, the distribution of skin-friction is derived from 
(128), which now becomes 

02 ^V (130) 



It follows that when x^C p (dC p ldx) <2> has been obtained at any station 
x, the value of T W /T B follows by solution of a simple quartic algebraic 
equation. 

For flows with an initial favourable pressure gradient (130) may be 
used to determine the distribution of skin-friction downstream of the 
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pressure minimum x = x mt provided T B is interpreted as the skin- 
friction of a Blasius boundary layer at a station x X L , where X L is 
given by (125). 

The method was applied by Curie (1960) to consider some of the 
cases for which exact solutions are known, and good agreement was 
obtained, both as regards the distribution of skin-friction and, in 
particular, the position of separation. 

7. Use of the energy integral equation. The methods of Tani and 
Truckenbrodt 

The energy integral equation, due to Leibenson (1935), is 

I^ = T' (131) 

, s Fut u*\, 

where S 3 = 1 -- -}dy 

J %\ ujj 
o 

oo 

and D = J~J Ay. 

o 
Upon writing a = 3 3 /S 2 , (132) 

D 

and = |, (133) 



(131) becomes 8 2 ~ (ottf8 2 ) = 2j8vwf. (134) 

dx 

Upon multiplication by 2aul , this becomes 



(135) 

ctx 

which integrates to give 

X 

X 81 = 4v J ocfiul dx, 



plus a constant of integration which is zero. Thus 

X 

31 = ~ ttf f ajSwf dx. ,^i36) 

o ^ /,. i *'" t '\,, 

Now Truckenbrodt (1952) suggested that in many cases it should be 
sufficiently accurate to regard a and /3 as constants.' It is fairly easy 
to see that a should not vary very much, as the following argument 
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indicates. Now a, by (132), is the ratio of the energy thickness to the 
momentum thickness, each of which are integrals with integrands zero 
at the wall. Accordingly little contribution arises from the region very 
close to the wall. But it is just in this region that the greatest changes 
of profile shape take place, and the outer part of the boundary layer 
retains approximately the same general shape regardless of pressure 
gradient. Hence S 3 /S 2 = a does not vary very much. This argument 
is borne out by the fact that for a stagnation-point a 1-63, for a 
Blasius boundary layer a 1-57, and for a typical separation profile 
a = 1-52, so that a varies by a maximum of about 3 per cent, about 
its mean value. Similarly the value of ft does not vary as much as might 
be expected. The reason for this is presumably that the decrease in 
the contribution to the dissipation integral as the skin-friction decreases 
to zero at separation is partially balanced by an increase in momentum 
thickness and a decrease in the local mainstream velocity. The values 
of ft are 0-209 for a stagnation -point, 0-173 for a Blasius layer, and 
0-157 for a typical separation profile, so that ft varies by no more than 
14 per cent, about its mean value. 

Since a and ft do not vary excessively it seems reasonable to replace 
them by their mean values, which will be close to the values appropriate 
to a Blasius layer. At positions downstream of the pressure minimum 
the most important contributions to the integral in (136) will arise 
from the region close to the pressure minimum, and the constant values 
chosen for a and ft will be particularly accurate in this region. With 
OL = 1-57, ft = 0-173, (136) becomes 

X 

81 = 0-441 vMf 6 J u\ dx, (137) 

o 

which is identical to Thwaites's formula (107) save for a 2 per cent, 
change in the numerical factor, corresponding to a change of 1 per cent, 
in the resulting value of 3 2 . 

Tani (1954) has developed a method of calculating laminar boundary 
layers which uses both the momentum and kinetic-energy integral 
equations. With the notation already used we may write the momen- 
tum equation in the form 

||(St) = M2+#)A, (138) 

and the kinetic-energy integral equation in the form (135) may easily 
be rewritten as -, /> 
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Hence, upon eliminating (djdx)($l), we find that 

Z - (#-l)A + ^-A^. (139) 

a a% 

We may regard (136) and (139) as new forms of the two integral 
relationships. Tani integrates (136) in the approximate form (137) to 
yield S 2 (and A) as a function of x. He then assumes that Z, H, a, ft are 
related to each other by the same relationship as holds for Pohlhausen's 
family of profiles, but determines the dependence upon A from (139). 
Thus, as a first approximation he ignores the last term and by replacing 
H, a, j8 with the relevant functions of Z, a relationship between Z and 
A is obtained.*)" Having thus related Z, //, a, /?, and A, these values can 
be used to estimate the neglected term, and a second approximation 
calculated. When this procedure has been repeated until sufficient 
accuracy is obtained, the relationship between a, /?, and A may be 
substituted into (136) to obtain an improved relationship between A 
and x, but Tani remarks this will not in practice be necessary. The 
universal relationships between Z, H, a, jS are shown in Table 6. 

TABLE 6 

Universal relationships for Tarn's method 
H l a p 

3-50 1-53 0-157 

3-33 0-029 -53 0-157 

3-18 0-059 -54 0-158 

3-04 0-089 -54 0-159 

2-92 0-120 -54 0-161 

2-81 0-150 -55 0-164 

2-72 0-180 -56 0-167 

2-63 0-210 1-56 0-170 

2-55 0-235 1-57 0-174 

2-49 0-260 1-58 0-179 

2-43 0-283 1-59 0-184 

2-33 0-322 1-60 0-193 

2-27 0-347 1-61 0-200 

2-25 0-356 1-62 0-203 

Truckenbrodt (1952) introduces an alternative shape parameter to 
H Sj/Sa and a = S 3 /S 2 - He writes (139) in the form 



(140) 



f It is interesting to note that according to this first approximation^^ m^wittfJB[ 
a, /? are functions of A alone, separation "occurs when A = -0-082, \s$l!eh i|^api- ., 
to the value appropriate to the known exact solutions than ia PohUhausen'-a '] 

A =-0-157. ^ <> -;; 
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and introduces a shape parameter $ defined by 

dS ^ d^dx 
dx ~~ a(H-l) 

and 8 = when A = 0. He also defines T(8) as 



so that from (139), putting A 0, T(0) = 0. Then (140) becomes 

fJ'v ti 11 A 

(JUdj M/1 li>1 A 

or ^-L^LT(S)==^. (141) 

dx A% % 

He then assumes that the relationship between T, S, H, Z, a, p is 
exactly that which holds for the Falkner-Skan similarity solutions, 
but determines the dependence upon A and x from the integral relation- 
ships as expressed in equations (137) and (141). The universal relation- 
ships are shown in Table 7. 

TABLE 7 
Universal relationships for Truckenbrodt's method 

S H I a 8 



-0-018 


4-04 





L-52 


0-157 


-0-015 


3-20 


0-095 


53 


0-159 


-0-010 


2-88 


0-149 


55 


0-163 





2-59 


0-220 


57 


0-173 


+ 0-010 


2-43 


0-277 


59 


0-185 


0-020 


2-32 


0-320 


62 


0-197 


0-030 


2-24 


0-357 


63 


0-208 


0-040 


2-16 


0-390 


66 


0-219 



By examination of the similar solutions, Truckenbrodt found that 
the function T(8) could be written as 

T = cS, c = 2-87 when 8 > 0, 
c == 3-53 when 8 < 0, 

so that (141) becomes ;r~"'l~T~^ ~ ~^' 

ax /\U-I u-t 

which may be written 
dx 



( "\ ( r \ 

L, r cu( j \ u\ r cu\ , 

.S'exp \ -^dx = -*exp i dx . 
I r J A% % [ J A^! j 

x 'y, ' v 'r, ' 
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Now it may easily be shown that, with A defined by (137), 



where 



. X . 

f cu l J ( x ) 
ex P ~dx\ = -(, 
(J A% J fto) 

v o?i ' 

f X V -, 

If I 

|(#) = j I u\dx\ 
\ ' 



c 
and d = - = 8-0 when 8 < 

0-441 . (142) 

= 6-5 when 8 > 
It then follows, after some algebra, that 

g^, (143) 



l 

The integration starts from the leading edge, where = 0, and if this 
is taken as the datum position, (143) yields 



If u^) is an increasing function of f, then S will be positive, whereas 
S will be negative if %() is a decreasing function of . The value of 
d is determined accordingly from (142), and the integration proceeds 
easily until a point is reached at which S changes sign. At this point 
changes discontinuously (because of the changed value of d), and the 
integration then proceeds by putting ^ in equation (143) as the down- 
stream value of at this point. 

It should be remarked that the two methods described above, due 
to Tani and Truckenbrodt, arose essentially as simplifications to a rather 
complex two-parameter method due to Weighardt (1946, 1948). A 
further method based upon that of Weighardt, designed in this case to 
improve rather than simplify it, is due to Head (1957 a), and we shall 
now consider this method. 

8. The method of Head 

This method was designed to be considerably more general than any 
of the preceding methods, in that velocity profiles are predicted with 
considerable accuracy even when distributed suction is applied. Head 
prescribed a two-parameter family of velocity profiles dependent upon 
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the parameters I and A. From the assumed profiles, chosen carefully 
to give the correct form in a number of representative cases, the various 
integral properties, a, /?, H, can be deduced as functions of I and A, and 
Head gives charts of these relationships. He also makes use of both 
the momentum and kinetic-energy integral equations, writing these as 



CLJC 

A = -- 

dx % d|/i> 

which are similar to (138) and (140) respectively. It is assumed that 
starting values of I and A are known. These will be I = 0-360, A 0-085 
at a stagnation -point, and I = 0-221, A at a sharp leading edge. 
From the known values of I and A, values of //, a, j3, are read from the 
charts, whence (144) and (145) yield d(Bl/v)ldx and d(at)[dx. It follows 
that, corresponding to a small increment in x, the increments in 8 2 and 
y. are found. Once the new value of 8 2 * s known, A follows. Then, by 
reference to the chart for a(A, Z), the values of a and A enable the new 
Falue of I to be read off. The values of // and [$ follow and the solution 
gtt the new point is complete. By repeating these steps the development 
;>f the boundary layer may be calculated by a step-by-step procedure. 
Head applied his method to several of the cases, considered in Chapter 
3, for which accurate solutions are known, and obtained good agree- 
ment. He also did some calculations (Head, 19576) for flows with 
listributed suction, and estimated the effects of distributed suction 
ipon boundary-layer stability. These points, however, are outside the 
scope of this book, and will not be considered further. 



COMPARISONS OF APPROXIMATE METHODS 
OF SOLUTION 

IN the preceding chapter no less than ten different methods have been 
given for calculating the development of a two-dimensional incom- 
pressible laminar boundary layer. In view of the fact that little more 
remains to be done in this particular field, and there would appear to 
be little room for any improved methods, it seems desirable to make a 
critical comparison of the various methods, so that readers who are 
interested mainly in obtaining numerical results will be able to con- 
sider more readily which of the methods best suits their requirements. 
There would appear to be three things which might determine the 
choice of a method. Firstly, the question arises as to how much informa- 
tion is required. For example, is it desired merely to calculate the 
position at which the boundary layer separates, or the momentum or 
displacement thickness at various stations ? Or, on the other hand, is 
it required to obtain an accurate estimate of the detailed manner in 
which the boundary -layer velocity profile develops? Secondly, the 
question of accuracy arises. Is a rough estimate, say within 15 per cent., 
adequate, or is it essential to calculate some property or other to within 
say 5 per cent, at most? Thirdly, when those methods have been 
rejected which do not match up to the requirements on these two 
counts, the choice between the remainder can be determined by 
questions of speed and simplicity. We will now consider these three 
points in turn, and will first group the methods in the following way. 

A. Methods based on the momentum integral equation: 

Pohlhausen (section 1) 

Timman (section 2) 

Howarth (section 3) 

Walz (section 3) 

Thwaites (section 4) 

B. Methods based on outer and inner solutions: 

Stratford and Curie (sections 5 and 6) 
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C. Methods based on the kinetic -energy integral equation: 
Tani (section 7) 

Truckenbrodt (section 7) 

Head (section 8) 

1. Information provided by various methods 

We begin by remarking that if only the position of boundary-layer 
separation is required, any of the methods is able to provide an estimate. 
If, further, the detailed distribution of skin-friction is sought, again any 
of the methods is satisfactory, with the qualification that for flow which 
starts from a stagnation-point rather than a sharp leading edge, the 
Stratford-Curie method is only used downstream of the pressure 
minimum. 

If any further information is required the Stratford- Curie method 
must be rejected, but any of the remaining methods are able to give 
predictions of the boundary-layer momentum and displacement thick- 
nesses. If, alternatively or additionally, the energy thickness is required 
the methods in group A fall down, strictly speaking. However, it would 
be such a simple process to make the relevant extensions to each of 
these methods that we may perhaps accept them from the practical 
point of view. 

The only remaining possibility (of any great likelihood) is the detailed 
development of the velocity distribution, and the methods of Thwaites, 
Tani, and Truckenbrodt fail to answer this requirement. 

2. Comparison of accuracy 

Before assessing the accuracy of the various methods it is necessary 
to consider carefully just what accuracy is desired in any given problem. 
On the one hand it should be pointed out that the position at which 
boundary-layer separation takes place is somewhat difficult to deter- 
mine experimentally. In this sense it may not be necessary, in certain 
circumstances, to have very great accuracy. On the other hand, if one 
is interested in a theoretical prediction of boundary-layer stability, 
extremely accurate detailed velocity profiles must be obtained. It 
follows that the ensuing remarks regarding accuracy must be interpreted 
in the light of the particular problem under consideration. 

A natural starting-point is to compare the methods of type A above, 
which are all reducible to a form in which three functions, I, H, and L, 
are assumed given as functions of A only. Of course, I, //, and L are 
not universal functions of A. In fact, for any given pressure distribution 
the exact solution of the laminar boundary-layer equations would yield 
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a set of functions. In a sense, therefore, the question of accuracy is one 
of luck, whether the assumed functions happen to be sufficiently close 
to those relevant to the case in question. 

An immediate consequence of this is that no one method of this type 
can be described absolutely as the most accurate or least accurate. 
For example, Pohlhausen's method would probably be considered by 
almost everyone to be the least accurate of the five methods of type A, 
yet there is no doubt that pressure distributions could be devised for 
which Pohlhausen's method would yield more accurate solutions than 
any of the other methods. In practice, however, there will be limits 
to the range of pressure distributions which might arise, and accord- 
ingly in the range of variation of I, H, and L as functions of A. It was 
on this basis that Thwaites calculated these functions for all known 
exact and approximate solutions, and then suggested values for general 
use which are 'average' in the sense of being far removed from the 
extreme limits of variation. 

One of the important results obtained by Thwaites was that the 
functions obtained from the various solutions agree well for favourable 
pressure gradients (A > 0), but differ as separation is approached in 
the presence of unfavourable pressure gradients (A < 0). Thus it 
matters little which of the methods is used from the forward stagna- 
tion point to the pressure minimum, but predictions of boundary-layer 
development downstream of this point will depend vitally upon the 
method used. Since the various curves do not differ much in their 
overall general shape, it is roughly true that the method which yields 
the best prediction of the position of separation will also tend to yield 
the best overall picture of skin-friction distribution. It is therefore 
natural to compare the predicted positions of separation for cases in 
which the position has been accurately determined. This has been done, 
not only for the methods of type A but also for the other methods 
discussed. The results, expressed non-dimensionally by quoting the 
pressure coefficients at separation, are given in Table 8. It appears 
from this Table that the most accurate methods (at least as regards 
predictions of separation) are those of Stratford, Tani, and Head, and 
that there is little to choose between these. It is also apparent that of 
the methods of type A, those due to Timman and Thwaites are the most 
accurate. 

With regard to predictions of momentum thickness it would probably 
matter relatively little which of the methods is used. In methods of 
type A the momentum thickness is determined mainly by assumptions 
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Pressure coefficients at 
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0-089 
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nade concerning the function L(X) in equation (102), and, as Thwaites 
howed, this function does not vary much between various exact and 
ipproximate solutions. Thwaites 's method predicts the momentum 
hickness to within 1 or 2 per cent, at all stations in the case of the 
inearly retarded velocity, but is in error by about 6 per cent, near to 
L stagnation-point. 

Predictions of displacement thickness are given from the momentum 
hickness upon multiplication by the shape factor U. Accuracy in pre- 
Licting displacement thickness is thus dependent upon accuracy in the 
r alue of H. Methods of type A assume that // depends only on A, and 
/re accordingly of uncertain accuracy. The best is likely, a priori, to 
>e that of Thwaites, since it is 'geared' to the known accurate solutions. 
Comparisons between Thwaites 's method and the exact solutions indi- 
ate that displacement thickness is usually predicted fairly accurately 
xcept near to separation where errors of up to about 10 per cent, can 
'ffcen arise. On the other hand, such comparisons as are available 
adicate that Head's method predicts displacement thickness to an 
ccuracy of say 2 per cent, or thereabouts. Although detailed com- 
arisons are not at present available, one would expect that methods 
ke those of Tani and Truckenbrodt would be intermediate in accuracy, 
Lnce some attempt is made to improve upon Thwaites's assumption 
hat H depends upon A alone. 

With regard to predictions of velocity profile, the only method which 
lakes any serious attempt to consider this is that of Head, and we 
ontent ourselves here with the comment that the accuracy achieved 
5 considerable. 
. Ease of computation 

By far the simplest method of those considered is that due to Stratford 
nd Curie. To determine the skin-friction at a given position it is neces- 
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sary only to calculate the value of x 2 C p (dC p ldx) 2 , after which the ratio 
of skin-friction to Blasius skin-friction follows by solution of (130), a 
simple quartic equation. 

The methods of type A, particularly that of Thwaites, are also fairly 
simple to apply. For this method, given the velocity u^x) at a number 
of stations, it is required that u\, u\, and du^/dx be calculated, and a 
simple quadrature (107) then yields 8|. The remaining characteristics 
then follow from prescribed universal functions. Other than Thwaites 's 
method, the methods of type A are rather longer, in that either non- 
integral powers of % are required or the quadrature is replaced by 
numerical solution of a differential equation. 

The methods of Tani and Truckenbrodt are also rather longer than 
that of Thwaites. Thus, although S 2 is calculated as a function of x 
exactly as in Thwaites 's method, the relationships between H and A, 
I and A, etc., must be estimated by approximate solution of the kinetic- 
energy integral equation. 

Head's method is the longest of all, as the simple method above for 
calculating S 2 is rejected, and the method requires the simultaneous 
step by step integration of two first-order ordinary differential equations. 

4. Conclusions 

In summing up what has been discussed above, we begin by remark- 
ing that for calculations of skin-friction alone the obvious choice is the 
method of Stratford and Curie, since this method is both as accurate 
as any other and considerably simpler than any other. By this method 
the skin-friction is given by equation (130). 

If, further, either the momentum or displacement thickness is re- 
quired it is suggested that these be calculated by Thwaites 's method, 
provided that errors of up to about 5 per cent, in 8 2 and 10 per cent, 
in S x are acceptable. If similar errors in skin-friction are accepted the 
whole calculation can be done by Thwaites 's method. 

If the above errors are rather too great, they may be somewhat 
reduced by using the slightly longer methods of either Tani or Trucken- 
brodt. The probability is that Tani's method is slightly more accurate 
and marginally simpler. 

Finally, if greater accuracy than this is required, or if velocity profiles 
are called for, it is suggested that the whole calculation be done by 
Head's method. 

One factor which has not been discussed above, which is fundamentally 
important and very relevant to any possible extensions of a method, 
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is the extent to which the approximations used correspond to physical 
reality. In this respect it should be remarked that the approximations 
of Stratford's method are physically realistic, and that Thwaites's 
quadrature for momentum thickness rests on a firm physical basis by 
virtue of the work of Leibenson and Truckenbrodt. 



SOLUTIONS OF THE TEMPERATURE 
EQUATION AT LOW SPEEDS 

IN low-speed flow, provided the difference between the temperature 
of the stream and that of the wall is not too great (so that the density 
is sensibly constant), the velocity and temperature equations may be 
solved separately, as was previously shown in Chapter 1, section 10. 
The equations governing the flow then become 

du , c)u du-i , ti z u ,-, Aft . 

u --- ^.v = Ui _L-J/_ (146) 

c)x^ r dy l dx^ % 2 v ' 



, c)T , dT v ,-._. 

and u --- \-v- = A -A --- , (148) 

* 



where the last term in (148), representing frictional heating, may 
usually be neglected. Then (146) and (147) may be solved to obtain 
the velocity field (u, v), the methods discussed earlier being available 
for this purpose, and the linear equation (148) for T remains. We now 
consider some solutions of these equations. 

1. Forced convection from a flat plate 

We begin with uniform flow parallel to a heated flat plate where the 
flow has zero pressure gradient, and the velocity % at the edge of the 
boundary layer is constant, U Q . Then the solution of (146) and (147), 
discussed in Chapter 2, section 1, is 



(149) 



where /( ??) is to be regarded as a known function. 

', (150) 



We write TT 
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where suffix w refers to values at the wall and suffix zero to (constant) 
values in the mainstream, and look for a solution in which 6 is a func- 
tion of <rj alone. Then, upon substituting for u, v, T from (149) and 
(150) into (148), we find (when the wall is at constant temperature) 

0"+or/0' = 0, (151) 

after neglecting the frictional heating term. The boundary conditions, 
T = T w when T? = 0, T -> T as rj -> oo, when substituted into (150), 

become 

0(0) = 0, 0->l as 7?->oo. (152) 

The solution of (151), subject to the boundary conditions (152), was 
given by E. Pohlhausen (1921), and is 

i? t 
r l 

= a (a) exp c 
J I 

(153) 

where ao( a ) = I J expl o- fj 

Since /( 17) satisfies (64) we may write 

/=-/"'//", 

whence (153) becomes 

(154) 






Pohlhausen calculated a (cr) for a range of values of a, which are shown 
in Table 9, with the approximate formula 

TABLE 9 
The function a (a) 

0-6 0-552 0-560 

0-7 0-585 0-589 

0-8 0-614 0-616 

0-9 0-640 0-641 

1-0 0-664 0-664 

1-1 0-687 0-685 

7-0 1-29 1-26 

10-0 1-46 1-43 

15-0 1-67 1-64 
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The local heat-transfer rate to the plate, 






2. Temperature of plate thermometer in moving fluid 

The problem here considered is that of determining the temperature 
which a flat-plate thermometer will read when no heat is being trans- 
ferred to it from the fluid. Viscous heating is not neglected, as it is the 
only reason why the wall temperature rises above that of the stream. 
Since T w is not prescribed it is convenient to define an alternative non- 
dimensional temperature, so we write 



(157) 

Then upon substituting into (148), including the frictional heating 
term, we find that ^ , , = __2af ff2 

the boundary conditions being 

-> as 17 ~> oo, 

and 0' = when 77 = 0. 

The solution, first given by E. Pohlhausen (1921), is 



J (D*\ j (/") 

T] *- 



so that the temperature of the wall is 

77 2 



(159) 



The function ^(0) has been calculated by Pohlhausen, an 
in Table 10, with the approximate interpolation formu|| 

0(0) = 4a*, 
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with which approximation (159) becomes 



TABLE 10 
The function 

a 6(0) 4a* 

0-6 3-08 3-09 

0-7 3-34 3-34 

0-8 3-58 3-58 

0-9 3-80 3-79 

1-0 4-00 4-00 

1-1 4-20 4-19 

7-0 10-06 10-6 

10-0 11-86 12-6 

15-0 14-14 15-5 

Before passing on, two points are worthy of mention. Firstly, when 

a = 1, (158) takes a particularly simple form, namely 

9(n) = 4-/' 2 

so that the temperature distribution is given by (157) as 



V 2 

or T + - = constant. (160) 

2c p 

This result, in fact, holds when o- = 1 for flow at arbitrary Mach number. 
Secondly, since the thermal energy equation is linear, solutions can be 
superimposed, so that for flow past a flat plate with heat transfer and 
frictional heating, the solution is a suitable linear combination of the 
solutions given in this section and section 1 of this chapter. 

3. Heat transfer near a stagnation-point 

Near to the forward stagnation-point of a cylinder the velocity 
distribution, as was shown in Chapter 2, section 2, may be written as 



(161) 



' = ta 

and a solution sought by writing 

T-T,., 



( 162 ) 
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where T is the temperature in the mainstream, constant in low-speed 
flow. Upon neglecting frictional heating, substitution from (161) and 
(162) into (148) leads to the equation 

O'+qfB' = 0, 
the boundary conditions being 

= when 7? = 0, 

6 -> 1 as 77 -> oo. 



and 

The solution is 



6(i,) = ^(a) J [expf-a //<*>?)] *1> ( 163 > 

ft L \ r , /J 



where 



This solution, due to Squire, is reported in the book by Goldstein (1938), 
together with values of a 1 (a) for a range of values of a. These results 
are shown in Table 11, with the approximate form 

ai (<j) == 0-570o- 0>4 . 

TABLE 11 
The function a 1 (cr) 
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0-9 

1-0 

1-1 

7-0 

10-0 

15-0 



The local heat-transfer rate to the cylinder in the region of the stagna- 
tion-point is 



a l(<*) 


0-570(7'* 


0-466 


0-465 


0-495 


0-495 


0-521 


0-521 


0-546 


0-546 


0-570 


0-570 


0-592 


0-592 


1-18 


1-24 


1-34 


1-43 


1-54 


1-68 



(164) 
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4. The solutions of Fage and Falkner 

The solutions given in sections 1 and 3 of this chapter are two special 
cases of a more general class of solutions, considered by Fage and 
Falkner (1931), As before we write 

TT 

jf _ w __ 

TT ~~ 
\ w 

where the wall temperature T w may be a function of x. Then upon 
substitution into the temperature equation (148), and again neglecting 
frictional heating, we find 

80 , 6-1 dT w , 80 d*B /lftKX 

<u, -- L-U -- 4~v = K . (165) 

8x h T W -T, dx + 8y Sy* ( > 

Fage and Falkner consider the cases in which the mainstream velocity 
is given by % = ^^ 

so that, as in Chapter 2, section 3, the velocity distribution is given 
in terms of a stream-function, and 



JL v -- r 

8y dx 



where 



is the solution of (71). Fage and Falkner also assume that the 
wall temperature takes the form 

T w = Zl+e(a/c). (166) 

Then upon substituting into (165), and looking for a solution in which 
is a function of rj alone, we find that 

8"+v[%(m+l)f6'-nf'(0-I)] = (167) 

with 0(0) = 0, 0(oo) = 1. Fage and Falkner integrated (167) for a range 
of values of m and n, for a Prandtl number a = 0-77, and making the 
approximation that 



that is, replacing the curve of u against y by its tangent at the wall. 
The value of a is, of course, proportional to the skin-friction. 

For the case of a wall at uniform temperature, n = 0, and (167) is 
integrable for arbitrary m without further approximation, as Squire 
has pointed out, the result being shown in equation (175). 
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5. Lighthill's method 

Lighthill (1950 a) generalized the work of Fage and Falkner, subject 
to their approximation to the velocity profile, to the case of arbitrary 
mainstream velocity and arbitrary wall temperature. Lighthill begins 
by making a transformation of the von Mises form, with new inde- 
pendent variables (%,$), where 

(168) 

is the mass flow between a given point (x, y) and the wall (x, 0). Then 
the energy equation (148), with frictional heating neglected, takes the 
form 

(169) 



dx V i 
Fage and Falkner 's approximation to the velocity, 

u = ?2$.y, (170) 

is written by Lighthill in the form 

U=:(^^ 9 (171) 

so that (169) becomes 

T i % i %nn\ 

(172) 



This equation was integrated by Lighthill, using the Heaviside opera- 
tional technique, and his value for the local heat transfer to the wall is 






where T e is defined as the difference between the wall temperature T w (x) 
and the mainstream temperature T (constant at low Mach number), 

that is T (x) = TV(X)-TI. (174) 

The integral in (173) is a Stieltjes integral, which may be regarded as 
a shorthand notation for 



o " o 

when the only discontinuity in g(t) occurs at t = 0. 
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Lighthill points out that the approximation (171) is essentially a 
high Prandtl number approximation. As the Prandtl number increases 
the thermal boundary layer becomes thinner, relative to the velocity 
boundary layer, so that for the purpose of calculating the thermal 
boundary layer it becomes increasingly accurate to replace the velocity 
curve by its tangent at the wall. Accordingly (173) is expected to be 
asymptotically exact as the Prandtl number tends to infinity. 

Lighthill tests the accuracy of (173) at a Prandtl number a = 0-7, 
by comparison with the similarity solutions. With a mainstream 

velocity ! = &,(*/')" 

and a uniform wall temperature T w , it can easily be deduced from 
(167) that 



(175) 



a function only of m and cr. Exact results for a = 0-7 and various 
values of m are compared in Table 12 with values deduced by Lighthill 
from (173). 

TABLE 12 

The function a m (0-7) 

OL m OL m 

m (accurate) (approx.) 

-0-0904 0-199 

-0-0654 0-253 0-232 

0-292 0-300 
+ | 0-331 0-360 

| 0-384 0-435 

1 0-495 0-587 
4 0-813 0-995 

It is immediately clear that agreement is excellent, about 3 per cent, 
error, when m = (zero pressure gradient), and is reasonable in the 
vicinity of a stagnation point, m = 1, where the error is about 19 per 
cent. On the other hand the agreement is poor near a position of 
separation, where m = 0-0904. Lighthill indicates the reason why 
the error varies with m in this manner. When m = the velocity 
boundary-layer profile is essentially linear over a considerable portion 
ill the boundary layer. For other values of m, where there is a curvature 
in the velocity profile at the wall, there is a greater deviation from the 
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tangent at the wall in the outer parts of the boundary layer. Thus 
when m > the assumed velocity is greater than the true value, so 
the heat transfer is overestimated. For m < the assumed velocity is 
less than the true value, and the heat transfer is underestimated. As 
separation is approached this effect is excessively accentuated, since 
LighthnTs assumption is equivalent to putting the velocity equal to 
zero everywhere. 

The aerodynamically desirable range is < m ^ 1, since negative 
values of m lead to the possibility of boundary -layer separation, and 
further have velocity profiles with a point of inflexion and a resultant 
possibility of instability. Now in this range the error in the predicted 
heat transfer varies between + 3 per cent, and +19 per cent. Lighthill 
points out that if the predicted heat transfer be multiplied by a factor 
0-904, so that the numerical factor 

-'{(*) I}- 1 = 0-538 

in (173) is replaced by 0-487, the error will lie between +7 per cent. 
Accordingly the accuracy is good when this empirical correction is made. 

6. The work of Liepmann and Curie 

Liepmann (1958) has derived an expression for the heat transfer 
which is identical in form to Lighthill 's result (173), and is obtained by 
use of the thermal energy integral equation, (55). This equation may 
be written in the form 

oo 

$- { u(T T 1 )dy= ^51, (176) 

dx J pc p 

o 

where Q w (x) is the heat transfer to the wall 



Liepmann, in place of Lighthill 's approximation (170), writes the 
velocity near to the wall as 






and by writing dy = dT, 

Q 

(176) may be transformed into the form 
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where 6 is the usual non-dimensional temperature as defined in equation 
(162). Liepmann then considers the special case in which the wall 
temperature is equal to its zero heat-transfer value T when x ^ x Qy 
say, and increases discontinuously to a new constant value downstream 
of this position. In such a special case Q/Q W is approximately a universal 
function of 9 alone, so that (177) becomes 

(178) 



where 



1 

a = f ~^ d6 (179) 

J Q/Qw 



is a constant. When r w is known as a function of a;, (178) may easily 
be integrated to yield Q w (x). Liepmann expressed this integral ana- 
lytically, and then derived the value of Q w (x) for an arbitrary wall 
temperature by integrating the contributions from a distribution of 
elementary changes similar to the above, and so obtained a result 
identical to (173), with the constant 9~*{(J)!}~ 1 = 0-538 replaced by 
(fa)*. We note that the value a = 0*215 is obtained by setting 
QIQ w = (l02)i jn (179), and this leads to a value (fa)* = 0-523, which 
is within 3 per cent, of LighthuTs value. Liepmann also shows how 
the same approach may be used to derive a formula which is valid in 
the vicinity of separation. 

Curie (196 la) has shown how, by an approach modelled on that 
of Liepmann, Lighthill's method may be considerably improved in 
accuracy. The thermal-energy integral equation (176) is used, and the 
velocity u approximated by 

p,u = r w y + - - y 2 . (180) 

This leads to the result 

. 00 

^ 

(181) 

In the special case of a region of zero heat transfer with T w = T 19 
followed by a discontinuity in wall temperature to a new constant 
value, the temperature profile is represented by a universal shape, 
characterized by a boundary-layer thickness 8 which varies with 
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position. It follows that 



and 



Thus y(T-TJdy = ~- 9 (182) 



and f y g (r-gi)dy = -<*~ (183) 

g Vu> 

where a and 6 are positive constants. Substitution from (182) and 
(183) into (181) yields 

d 



dx\ Ql 2Q 3 W I pkv * W9 

which is to be compared with Liepmann's equation (178). An alterna- 
tive way of writing this equation is to introduce a representative velocity 
f/oo and a representative length I, a local Nusselt number 

*=-2r=5r 

and a Reynolds number 



Then the value of K = 

is given by solution of the equation 



or KdW = ~-jL, (185) 



where A = aa 

and B = 1 



. 

dx 
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The constants a, b are determined by reference to the similarity 
solutions. The values 

a = 0-2226, b = 0-1046 

are such that agreement is within about i 1 per cent, over the whole 
range of pressure gradients from stagnation-point to separation, and 
with values of the Prandtl number of order unity (typical of gases) and 
order ten (typical of liquids). The numerical accuracy is, accordingly, 
an improvement on that given by Lighthill's method, although the 
presence of two terms on the right-hand side of (184) or (185) precludes 
analytic integration. A simple method of numerical integration, due 
to Thwaites, is to replace the integral in (185) by its Simpson's rule 
value, so that when K(x) and K(x-\-Ji) are known, K(x+2h) follows 
by solution of a quartic algebraic equation. 

It has been remarked earlier that the above approaches are 
essentially high Prandtl number approximations. Curie (196 la) has 
also examined the possibility of deriving acceptable solutions at Prandtl 
numbers of order unity by empirical correction of a low Prandtl number 
solution. At low values of the Prandtl number the temperature 
boundary layer is thicker than the velocity boundary layer, and Curie 
replaces the local velocity u in (176) by f$u l9 where j8 will tend to unity 
as the Prandtl number tends to zero. He also writes 



with a Q constant, in the special case of a single discontinuous jump in 
wall temperature. Then substitution into (176) yields 

^J^L\ = V Q 

dx\Qj pa,Wa(T w -T^ 

which may be integrated and transformed to yield 
Q w (x)= _ 



The generalization to the case of arbitrary distribution of wall tem- 
perature is exactly as described by Liepmann, and yields the Steiltjes 
integral x . x 

Q w (x) = -fc(^)x(*) J J i 

where T is defined as T W T 19 as in (174). The form of this equation 
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is in agreement with that obtained by Morgan et al. (1958), who showed 
that the asymptotically exact value of (i/3a )* as <r - is 

(Jj8 )* - 77-* = 0-5642. 

Curie remarks that when this constant is changed to 0-313, the formula 
(186) yields values of the heat transfer agreeing with the similar solu- 
tions to 25 per cent. Provided that the velocity boundary layer is 
not too thin in comparison with the thermal boundary layer Curie shows 
that a considerable improvement may be effected by writing 

(J]8a )* = 0-5642-0-265o-(l 5A), 

where A is the pressure gradient parameter, defined by (108), which 
may be calculated by means of equations (107) and (108). 

Curie has applied these two methods to the case of flow past a heated 
cylinder, and compared his results with the experiments of Schmidt 
and Wenner (1941). The agreement is good in each case, the solution 
(186) (derived from the low cr approximation) giving particularly good 
agreement. 

7. Spalding's method 

An alternative method of improving Lighthill's method was given 
by Spalding (1958). His method begins with the equation (180) for u, 

I dp 9 



and he remarks that the correction to Lighthill's method will depend 
upon the relative importance of the second term on the right-hand side. 
Now at a position which is representative of the thickness of the thermal 
boundary layer, say 



the ratio of the second term to the first is 



Accordingly he considers the case of a single step in wall temperature, 
writes (178) as 

apfc 3 (2L--5P,) 8 dfrl\ x A 

-- r v * u -f = constant, 
P* T*, dx\Q 3 J 

and suggests that the right-hand side of this equation be replaced by 
a function of the parameter % defined in (187), so that it becomes 
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where x * s the expression (A 4 S 4 /v)(<Z%/(fa;) of Spalding's paper. The 
function F is determined by reference to the available exact similar 
solutions. It is shown as a graph in Spalding's paper, and is tabulated 
here in Table 13. Spalding suggests the solution of (188) by iteration, 

TABLE 13 

The function F(x) 



x 

-4 3-5 

-3 3-8 

-2 4-3 

-1 5-1 

6-4 

+ 1 8-5 

2 l]-6 

3 15-8 

so that a first approximation to Q w , obtained by setting F = 6-4, is 
used to determine a better approximation to F, and so on. For the case 
in which T w = T : when x ^ f , and T w = 2 7 1 +A!T c (f ) when x > f , this 
yields the result 

Q w (x) = - 

where xi is the value of x obtained by using the previous approximation 
to Q w . The solution for arbitrary wall temperature is obtained by 
adding the solutions for the relevant distribution of the above ele- 
mentary increments in wall temperature. 

It should be remarked at this stage that, although the method will 
work effectively (and yield an improvement in accuracy as compared 
with Lighthill's method) provided that F does not differ too much 
from 6-4, it will still break down as separation approaches and the 
parameter ^ -> oo. 

Spalding goes on to show that if only the total heat-transfer rate is 
required, it is not necessary to calculate (and then integrate) the local 
values. A quicker method is to set up an equivalent relationship to 
(188), which will allow this to be done directly. Spalding takes the 
thermal energy integral equation (55), and integrates it with respect 
to x to obtain 



j Q w (x)dx = ~pc p j u^- 
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where A 2 is a measure of boundary -layer thickness. Considering, as 
usual, the single-step distribution of wall temperature, Spalding shows 
that (173), which is Lighthill's approximation, may be written as 

A* = 0-725, 



and that the appropriate generalization may be expressed as 

?. {K A 2 )H == 0-7-0 

(t JLT w) dx 
The first step in an iterative solution of this equation is 

f r / \ i 09? ' U 

%A 2 = 0-7- [Tv>\*dx KA 2 )^^^ , (190) 
( "J If*/ * J r w J 

where again a breakdown occurs as separation is approached. 

8. The method of Davies and Bourne 

We consider now briefly one further method which has been sug- 
gested for improving the accuracy of Lighthill's method, due to Davies 
and Bourne (1956). Their method is to assume that the velocity profile 
u(x, i/j) is known, and may be approximated by the expression 

so that u(x, y] oc yftv -#. (192) 

The values of 6 X and /3 are determined by the conditions of the parti- 
cular problem under consideration, so that the exact and approximate 
velocity profiles are in good agreement over the greater part of the 
thermal boundary layer. As o-~>oo we have j8 = f, as in Lighthill's 
approximation. When considering boundary-layer flow with a of order 
unity, Davies and Bourne chose b 1 and j8 to give good agreement over 
most of the velocity boundary layer, when applying their method to 
the case of the similar solutions. The temperature equation in the von 
Mises form (169), after substituting for u from (191), yields 

8T _ W h 8 (jj a d/1 
dx cr 1 1 dijj\ 

which is a form amenable to analytic treatment. Davies and Bourne's 
results indicate that the percentage error in the predicted heat-transfer 
rate is roughly equal to a typical percentage error in the assumed 
velocity profile. 

The method is neat, but does not appear to be of general use, since 
it relies upon a knowledge of the detailed velocity profile, which cannot 
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easily be estimated with any accuracy. Further difficulties arise if in- 
sufficient evidence is available concerning the relative extents of viscous 
and thermal boundary layers. 

9. The Meksyn-Merk method 

We now turn our attention to a method which was devised by Meksyn, 
applied originally to a number of viscous boundary -layer problems, and 
subsequently adapted by Merk (1959) for calculating thermal boundary 
layers. The applications to viscous boundary layers will not be discussed 
here, and the reader is referred to the book by Meksyn (1961) for 
details and for further references to the method. 

The method as used to calculate thermal boundary layers may be 
illustrated by reference to the problem of section 1 of this chapter, 
namely forced convection from a flat plate. Consider the function a (cr), 
given by (153) as 



U \ / Q \ 

We introduce the function F(r]) defined as 

T? 



1 

so that ,-r = f exp(-erJF)<fy. (193) 

<*o(<r) I 

Now by (65) we have 

f(rj) = 1 a ry*_I V + oPij*..., (194) 



so that F(VI) = --, aT? 3 - (XV + - c*V > ( 19r >) 

3! 6! 9! 

where a = 1-32824 and the series both converge for a finite range of ??. 
The series (195) may be inverted to give 



Then, upon substitution into (193) we have 

00 

-L = Q* J ^(-^ftf -i+^t-^ft ...}dF. 
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This series may be integrated formally term by term, and when the 
reciprocal of the resulting series is calculated it is found that 



-}- < 196 ) 

The limited range of convergence of the series (194) means that (196) is 
in the nature of an asymptotic expansion. However, even the three 
terms shown here are sufficient to give excellent results when a ^ i- 
When a = |-, these three terms yield a (0-5) = 0-520, compared with 
the exact value 0-518, and when a = 0-8 the three terms yield three 
correct decimals in O; O (<T). 

The above indicates in general terms the principles on which the 
method is based. For more general cases Meksyn and Merk make first 
a transformation of the variables which is virtually identical to that 
used by Gortler (1957 a), so that the equations for the velocity and 
temperature functions take a form in which the above type of asymp- 
totic expansion can be applied. 

By developments along these lines, it has been possible to consider 
the boundary layer for Schubauer's experimental pressure distribution 
(Meksyn, 1956), and for retarded flow along a semi-infinite plane (Mek- 
syn, 1950), the exact solutions to these problems being given by Hartree 
(19396) and Howarth (1938) respectively. For the thermal boundary- 
layer problem Merk (1959) has derived solutions corresponding to the 
similar solutions of E. Pohlhausen and Squire, considered in sections 1 
and 3 of this chapter. He has also considered more general cases of 
flow past circular and elliptic cylinders. For the case of a circular 
cylinder, the predicted heat transfer over the front 70 of the cylinder 
agrees exceedingly well with experiment, but the method breaks down 
as separation is approached at about 80. Meksyn (1960) has also shown 
how accurate solutions for the plate -thermometer problem (section 2) 
may be obtained for large values of <r(> 100) where straightforward 
calculations of the integral (159) are inaccurate. 

10. Curie's analysis by Stratford's method 

It has been remarked earlier that one of the simplest and most 
accurate methods devised for calculating the skin -friction in an incom- 
pressible laminar boundary layer is that due to Stratford (1954), the 
method being based on the idea of dividing the boundary layer into 
inner and outer portions. Curie (19616) has shown Jti$^'a ^i^njjar 
technique may be used for calculating the correspon|^n)g 'property for 
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the thermal boundary layer, namely the heat-transfer rate, when the 
wall is at uniform temperature. 

Following the procedure adopted in the case of the viscous boundary 
layer, the thermal boundary layer is divided into an inner and an outer 
region. In the outer region the thermal energy equation is written as 



where (d/ds) denotes differentiation at constant i/f, defined by (110). We 
also expand T(x, \jj) in the form 

x^. (198) 

We consider the special case in which there is zero pressure gradient 
upstream of a position x = X Q . The downstream velocity profile is 
distorted by the pressure gradient, but only near to the wall initially. 
Accordingly we deduce from (197) that the value of (dT/ds)^ is un- 
affected by the pressure gradient, except in the region near to the wall, 
so (198) shows that when > i/fy, say, T(x,\fi) is exactly as it would 
have been in the absence of pressure gradient, provided terms of order 
(x # ) 2 are neglected. Now it may be readily deduced from Pohl- 
hausen's solution (154) that in the absence of pressure gradient 

T B = T^+^-T^r,), 
and (anticipating the approximation below) this may be written 

T B (x, I) = T w +ai(T 1 -Tjpl^)* (199) 

\ P U I 

for small values of $, where a* arises from the approximation (155). 
There is no a priori reason why fa (determining the join in the inner 
and outer velocity profiles) and ^ (determining the join of the tem- 
perature profiles) should be equal. It is assumed, however, that both 
are sufficiently small that the Blasius velocity profile remains sensibly 
linear throughout ijj ^ ^ and ^ ^. It follows from (199) that the 
value of T(x, fy) given from the outer profile is 



T(x, #) = 3 T w +a*(r i -2'j (200) 

and upon differentiating with respect to t 
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Turning now to the inner region, ifj ^ i/^-, we assume a form 



which should be a good approximation, since both 8 2 T/8y 2 and 

are identically zero at y = 0. The values of T and dT/dy at the join, 

as given by (202), may be equated to the values (200) and (201), yielding 



(203) 
and gg = gi^TJUjl TB , ) ? . (204) 



These two equations connecting Q W9 $, Uj 9 y] must be supplemented 
by relationships giving $ and u'j in terms ofy'j. To do this it is assumed 
for simplicity that the expression (126) is adequate throughout ifj ^ fy. 
This means, in reality, that $ <J i/^ but since the inner and outer 
solutions for u have continuous values of u, Su/dy, and d^ujdy 2 , it is 
possible that (126) holds for values of ifj rather greater than ^. Then 

(205) 



and ^; =Ttr ^ + ^ + __^ + i, (206) 

with T M , and a known from the solution for the velocity profile. The 
value of Q w may be determined from (203) to (206), and is 

L//TT _ rp \ i \i 

Q = a i 5l_i^ T jQ-6642+0-3358 ^ *. (207) 

^o I T B) 

Accordingly, once the skin-friction r w has been determined, the heat- 
transfer rate is given by this explicit formula. This, of course, is by 
far the simplest method of calculating heat transfer. The limitations 
of the method are that it can only be used when the wall temperature 
is uniform and the pressure gradients are adverse. 

To estimate the accuracy, we remark that the formula (207) is exact 
when there is zero pressure gradient. For small pressure gradients the 
accuracy remains good. Thus, if the mainstream velocity is 
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the heat-transfer rate is predicted from (207) and (130) to be 



whereas the exact solution replaces the coefficient of x/c by 1-86. We 
should expect the accuracy to decrease as we proceed farther from 
the accurately known starting conditions, that is as separation is ap- 
proached. At separation (207) yields 



Q w = 0-815** i~ TB , (208) 



the numerical coefficient being equal to 0-72 when a = 0-7. For the 
similarity solution which corresponds to a separating boundary layer 
the numerical coefficient is 0*600(%/w )J, which might typically equal 
about 0-57. The agreement is still fairly good. 

11. Squire's method 

We shall now give a brief discussion of a method due to Squire (1942), 
modelled on the integral methods first used by K. Pohlhausen (1921) 
for the viscous layer. Squire's method rests initially upon the idea that 
the 'shape' of the temperature profile in a laminar boundary layer does 
not vary much over a wide range of pressure gradients. Since the tem- 
perature profile is the same as the Blasius velocity profile in the case 
of zero pressure gradient, we may suppose therefore that 



, T-/8604JA /oim 

where = ~|/' s - (210) 



u T-//0- 
= ~|/' 
% \ 



is the Blasius velocity profile, %f'(rj) being tabulated in Table 1. We 
then assume that the velocity profile is always given by (210), with 
the displacement thickness 8 X suitably chosen. Squire chooses 

M 2-960^ f 5 , /011X 

8 i == -^6~- "I**' ( 211 ) 

% J 

which is the value given by an approximate method due to Young and 
Winterbottom (1942). Then, substituting from (209) and (210) into the 
thermal energy integral equation (55), we have 



J 



dy = 
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which may be written as 

1-143* 



d A! ,,/A! \ n lff 

Sp o^ol J ' feT-*' 



M r ~ \ i ~ /ri c\\ 

-{2 Wl A^} = -^-, (212) 

where ^(AJSj), which was calculated by Squire, is defined as 





Upon integrating (212) we find that 

ul&l<j>*= 1-143/c 

o 
and when divided by (211) this yields 



A; ;AA 

~' 



which is an equation for Aj/c^. 

Squire suggests that a good first approximation is to take < as a 
constant on the right-hand side of (213), which may then easily be 
calculated as a function of x. Then the values of AJSj. may be deter- 
mined from a table of the function (Af/8f ^(AJSj), which is given here 
as Table 14. A second approximation, if required, is obtained by using 
the resulting values of ^ to recalculate the right-hand side of (213). 

TABLE 14 

The function ?? 2 <(7?) 

-n 4>W WU) 

i 0-208 0-052 

I 0-257 0-100 

f 0-272 0-121 

f 0-332 0-230 

1 0-386 0-386 
f 0-456 0-713 
*f> 0-499 1-018 
f 0-548 1-522 
|f 0-575 1-901 

2 0-599 2-398 

Squire compares the results of his method with some of the exact 
'similar' solutions over a range of values of a. He finds agreement 
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better than 1 per cent, as regards heat transfer in the absence of pressure 
gradient, but the predicted heat transfer is about 5 per cent, too low 
near a stagnation-point. A similar comparison for the case of the 
separation profile indicates that the predicted heat transfer is 30 per 
cent, too high. Squire has also used his method to calculate the heat 
transfer in flow past a circular cylinder, obtaining reasonable agree- 
ment with experiment. 

12. Free convection from a heated vertical plate 

To conclude this chapter we consider the problem of free convection 
when a heated plate is placed in still air. It is assumed that the tem- 
perature difference between the plate and the air is small, so that the 
fluid properties may be taken as constant, but that the (small) motion 
of the air is caused by a buoyancy force due to density variations. 
Then, with x measured upwards from the lower edge of the plate, y 
normal to it, u, v, the associated velocity components, T w , T Q , the 
temperature of the plate and the air respectively, the equations become 

^ + ^ = 0, (214) 

3u . 3u 3 2 u , 



, dT , dT v JL /c .~ _. 

and u \-v = K -. (216) 

3x ^ dy 3y* v ' 

np /TT 

Upon introducing 6 = -, 

(215) and (216) become 

u ~-\-v = v-^ + g w ~~ (217) 

, 80 , 36 3*6 /010X 

and u \-v = AC -. (218) 

3x dy dy* 

These partial differential equations may be reduced to ordinary differen- 
tial equations, as was first shown by E. Pohlhausen (1921). We write 

/ ^ 



8J> 
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so that u = vC*x*f'(i)), v = v<7ar-i(7]/'--3/). (219) 

Then (217) and (218) become 

f+2ff'-2f'*+g = 



(220) 

with boundary conditions 

0. a = L >n = Q \ 

(221) 



/'= = 0, ,,-*oo 

These equations were integrated numerically by Pohlhausen, with the 
Prandtl number a = 0-733. From his results we deduce that 

*- (222) 



giving the heat-transfer rate at the plate. 



THE COMPRESSIBLE LAMINAR BOUNDARY 
LAYER WITH ZERO PRESSURE GRADIENT 

THE essential difficulty in dealing with compressible laminar boundary 
layers is that the momentum and thermal energy equations are coupled, 
since the temperature field influences the density, which in turn in- 
fluences the velocity field. Mathematically we are faced with the im- 
mense task of solving two simultaneous non -linear partial differential 
equations. For this reason the only flow which has been considered in 
detail over a wide range of conditions is that when there is zero pressure 
gradient. For such a flow the velocity u l9 the temperature T^ and the 
density p x outside the boundary layer are constant. Each of these 
quantities is thus equal to its value at the leading edge, x = 0, and we 
shall accordingly denote them by suffix 0. The boundary-layer equa- 
tions, (1), (12), and (19), then become, upon neglecting body forces 

~(pu) + ^( P v) = 0, (223) 

BU . 8U 3 I dU\ /nnA\ 

pu -- \~pv- = ILL , (224) 

r dx r dy dy\ dy) v 

and 

(225) 

Before proceeding farther it is necessary to discuss what assumptions 
are to be made concerning the values of p, and a, and their dependence 
upon the various other properties of the fluid. 

1. Values for viscosity and Prandtl number 

An accurate representation of the variation of viscosity with absolute 
temperature is given by Sutherland's formula 

(226) 



where c is a constant, taken to equal 114 K for air. The complexity 
of the compressible laminar boundary -layer equations is such that it is 
often necessary to use a simpler formula such as 

H oc T>. (227) 
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By suitable choice of o> and of the constant of proportionality it is 
possible to make (227) agree well with the accurate value (226) over 
a fairly wide range of temperatures. For example, o> = 0-76 has often 
been used, and Cope and Hartree (1948) pointed out that good agree- 
ment in the range 90 K < T < 300 K could be obtained by taking 
o = 8/9. Considerable simplifications in the boundary-layer equations 
can be effected by taking o> = 1, and this fact has been used to good 
effect in many theoretical investigations. 

The value of the Prandtl number, cr, does not vary much with 
temperature for air, and a constant value of about 0-72 or 0-73 has 
usually been assumed. It must be pointed out, however, that con- 
siderable simplifications result from the choice a = 1, and theories have 
often been developed on this basis. 

We can easily see from equation (225) that when a = I the thermal 
energy equation always has a solution 

c p T+$u 2 = constant. (228) 

This solution corresponds to the case of a thermally insulated wall, 
whose temperature is constant when cr = 1 . For the case of zero pressure 
gradient we see further, by comparison between (224) and (225), that 
(225) always has a solution 

c p T-\-%u 2 = (constant)^, 
provided the wall is at uniform temperature. The general solution is 

then 



and since T = T W9 u = 0, when y = 0, and T T Q , u = U Q , when 
y = oo, this becomes 



or T + --V.+(T.-T,'l-, (228) 

4C p UQ 

introducing the temperature T z at which there would be zero heat 
transfer at the wall. 

We note that when a = 1, regardless of the viscosity-temperature 
relationship, the solution of (225) is given by (229) for the case of zero 
pressure gradient, and it remains only to solve (224) for u. We now 
consider various attempts to do this, with different assumptions regard- 
ing viscosity, as well as some solutions with values of cr other than unity. 
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2. The solutions of Busemann and Karman 

Busemann (1935) assumed a = 1, and so deduced as above that T 
is a function only of u. He then reduced the equation (224) for u to 
an ordinary differential equation with 77 = \y(u^v w x}^ as independent 
variable, the form being similar to that of the Blasius equation. This 
equation was then integrated numerically for the case of Mach number 
M = 8-8 and zero heat transfer at the wall. 

As a result of these calculations he found that the velocity varies 
almost linearly with distance from the wall at this Mach number. Since 
a linear velocity profile is a fair approximation even in incompressible 
flow von Karman (1935) considered the consequences of assuming a 
linear velocity profile at all Mach numbers, and calculated approximate 
solutions by the Karman-Pohlhausen integral method. Zero heat 
transfer was taken, of course, and the Prandtl number was set equal 
to unity, so that the temperature profile was given by (229), with 
T w = T s , namely, a 



or 



Now when there is no pressure gradient, the momentum integral equa- 
tion (40) may be written as 

d 



Upon substituting for p, 

/>= 

setting 

and 

we deduce that - Po^. (230) 



If we let the viscosity law be /x oc T^ 9 then fi w is related to the mainstream 
value JLA O by 
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so that (230) integrates to give 

8 



from which we deduce that the skin friction is 
= = f 

yiw \ I 

(231) 

This is effectively the form in which Karman presented his results. 
We note, however, that it may alternatively be written as 

(232) 



We note that (i/(A)}* increases from 0-290 when A == 0, i.e. M = 0, 
to 0-393 when A = 1 ? i.e. Jf ->oo. The accurate values are 0-332 and 
0-400 respectively, so the accuracy is considerable at high Mach numbers, 
decreasing from 13 per cent, at zero Mach number to less than 2 per 
cent, as the Mach number tends to infinity. 

3. The solutions of Karman and Tsien 

In view of the relatively large error at smaller Mach numbers, 
Karman and Tsien (1938) later attempted a more accurate solution of 
the equations of motion. They began with the boundary -layer equations 
in the von Mises form (59), (60), and in the absence of pressure gradient 

these become , 

/ooo\ 

(233) 



and ^T+l*) = +l*. (234) 

Upon taking or = 1, the solution of (234) is given by (229) as 

T + - = T v +(T.-T v )* (235) 

Abp MQ 

so that, once the wall temperature T w has been prescribed, the tempera- 
ture distribution is given explicitly in terms of u, which may then be 
determined from (233). To solve this equation, Karman and Tsien 
introduce non-dimensional quantities, 



m 

T Q 



and 
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where suffix zero refers to (uniform) values in the mainstream, and c is 
a representative length. Then (233) becomes 



* 

which equation may be solved by iteration. As a first approximation 
the function u*(\}f*) appropriate to the incompressible Blasius solution 
is taken, with T* deduced from (235) and //*,*, from the relationships 

ju* = (T*) ' 76 , (237) 

and p* = 1/T* 9 (238) 

this last result following immediately from the equation of state (7), 
since p is constant. Hence the product 



is known. Equation (236) then yields 

u* = C I j d$*, (239) 

J J 
o 

(/ 

J 

and C- 1 = 



This improved value of z(-*(i/f*) may then be used as a starting-point 
for the next approximation. When u*($*) has been determined to 
sufficient accuracy, the value of the true coordinate y normal to the 
wall is calculated as ^ ^* 

- 



Karman and Tsien calculated solutions for flow with zero heat transfer 
over a range of Mach numbers from to 10. The values of the skin- 
friction coefficient are shown in Fig. 1, together with the results of 
other workers. They also carried out some calculations for cases in 
which the wall temperature was kept constant at one-quarter of the 
mainstream temperature. For details of these results reference may be 
made to the original paper. 
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4. The solutions of Emmons and Brainerd 

Some solutions for zero heat transfer, with various viscosity -tempera- 
ture laws and at several values of the Prandtl number, were given by 
Emmons and Brainerd (1941, 1942). These authors introduced a num- 
ber of non-dimensional functions, defined as 

T 



u ^ I x U 

= , v* = v r, 
u \u Q vJ 



JM* 
= y \^} ' 



and showed that the equations of motion (223) to (225) could be reduced 
to three ordinary differential equations for %*, T* } and 



as functions of 77. These equations take the form 

(fit \ AT* 
|?_ w *| = ?_, (241) 

dr] } dr] 

f^!+2T* L* 1 = 0, (242) 

dr] dr]\ drj } 



and ^- w *+?AL*^\+2(y-l)JfjT*|~^l = 0. (243) 

These equations were then integrated on a differential analyser under 
the following assumptions. 

Firstly, a range of solutions was obtained for //,*=!, that is for 
constant viscosity and conductivity, equivalent to taking a> = in the 
viscosity-temperature relationship (227). These solutions were for 
values of the Prandtl number <T = 0,0-25, 0-733, 1-00, 1-20, respectively, 
and at Mach numbers Jf = 0, 0-5, 1, 1-5, 2, and VlO. Some of the 
results for the variation of skin-friction with Mach number are shown 
in Fig. 1. An interesting conclusion which may be drawn from the 
results of these calculations is that when there is zero heat transfer the 
relationship between wall temperature T w and mainstream temperature 
T , is given with considerable accuracy by 

T w = T {l+J(y IJJfefgor*}, (244) 

which is in agreement with the result obtained by E. Pohlhausen for 
low-speed flow. 

In their second paper, Emmons'and Brainerd initially took a = 0-733 
throughout, and examined various possible viscosity-temperature laws. 
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They concluded that the choice co = 0-768, close to the value (237) 
used by Karman and Tsien, was the most satisfactory, and subsequently 
used only this relationship. They then calculated velocity and tem- 
perature distributions at various Mach numbers, and with the same 
Prandtl numbers as they used in their first paper. Some of these further 
results have also been included in the skin-friction coefficients in Fig. 1. 




FIG. 1 . Skin friction on a flat plate with no heat transfer. 

Calculated by various workers; - - - according to equation (272). C & H, 

Cope and Hartree; K & T, von Karman and Tsien; B & E, Brainerd and 

Emmons. 

5. The calculations of Grocco 

An exhaustive investigation of the laminar boundary layer with zero 
pressure gradient has been carried out by Crocco (1946). Upon setting 
the pressure gradient equal to zero and applying the Crocco transforma- 
tion (Chapter 1, section 11), the equations (57), (58) become 



(245) 



and 



n _ \^^ r JL.(" * 

^ cry T -f-l - 



(24c) 



the dependent variables being r and /, the independent variables x 
and u, Crocco looks for a solution in which the velocity and tempera- 
ture profiles are each similar for different values of x (though not similar 
to each other). Thus / is a function of u alone, and r may be expressed 

as T = <f>(x)g(u). 
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After some straightforward analysis, Crocco deduces that (f> oc #-*, and 
accordingly writes / x 

r = gi, (247) 

whence substitution into (245) and (246) leads to 

gg'+ppu = 0, (248) 

and (I"+a)g+(lo)rg f = 0, (249) 

where dashes denote derivatives with respect to u. Non-dimensional 
quantities, similar to those used in previously discussed methods, are 
introduced, so that 

u j^ __ I # __ p 

T ~' Po ' ^ '^ . (250) 

/ *> \ i 
and 



Equations (248) and (249) then take the form 

*u* =.0, (251) 



and (/**+ cr ^2 \G+(l~~-a)I^G r = 0, (252) 



primes now denoting derivatives with respect to u*. The boundary 
conditions are 

G' ^0, /* = ^ at ^* = 0, 
*o 

and (2 = 0, /* = 1 at u* = 1. 

We note that p* = (I*)- 1 

and /i* (/*) w 
are functions of /* alone. 

Crocco considered first solutions with a = 1. When this holds 

(252) becomes simply 2 



so that the relationship between J* and u* is quadratic, as has pre- 
viously been indicated by other methods. It then follows that (251) 
takes the form GG+h( u *) = (253) 

where h is a known function of u*. This equation was solved by Crocco 
by a similar iterative method to that used by Karman and Tsien when 
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they solved (236). Thus we have 

1 ( 
G(u*) == f du* du*, (254) 



( \ 

== f f ^ du*\ 

^ ' 



and a first approximation 19 substituted into the right-hand side of 
this equation, leads to a second approximation G z . Crocco remarked 
that this procedure would not, in general, converge but that it would 
do so rapidly if either %(Gi~\-G%) or (G l 6? 2 )* were taken as the starting 
value for the next iteration. 

Having thus considered the case of a = 1, Crocco turned to the case 
o} = 1, a ^ I. When this holds equation (251) takes the form 

GG"+2u* = 0, (255) 

so that G is independent of /* and a. It is fairly straightforward to 
transform this equation so that it becomes identical with the Blasius 
equation for incompressible flow. The required transformation is 



whence (255) may be shown to become 

r+ff* = o. ( 257 ) 

The boundary conditions are easily shown to be such that the relevant 
solution of (257) is that given in Chapter 2 for the Blasius profile. It 
follows from (255) that G as a function of u* is simply the known 
relationship for \j" as a function of /'. 

In particular, the skin -friction is given by (247) and (250) as 



and by (256), substituting the value of /"(O) from (67), this yields 

( 258) 



x 
which may alternatively be written as 



(259) 

w w 

since pec T. 

Once this relationship is known, it remains to solve (252) for /*. 
Since this is a linear first-order equation for /*', it is a straightforward 
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procedure to obtain /* by integration, and amongst the various forms 
in which the result can be expressed is 

I*(u*) = l-I^OJA^u^+a^Aafott*), (260) 

J o 

where A 1? A 2 are two functions, calculated by Crocco, and shown in 
Table 15 for various values of cr. They are defined as 





In view of (256) these may alternatively be expressed in terms of the 
Blasius functions. Crocco showed, in fact, that 8crA 2 is the function 
0(77) denned by equation (158) and calculated by E. Pohlhausen (1921). 
An alternative method of expressing (260), which is more immediately 
useful when /*'(0) is unknown, is 

where A 3 = 1 -ifo^ ' (262) 

Ai (cr, 0) 

and A 4 = cr{A 2 (cr,^*)- A 2 a? ( A^o-,?/*) . (263) 



For the more general case, when a> ^=- 1, a ^ 1, it is necessary to 
integrate (251) and (252) simultaneously, and this was done by Crocco 
for a large number of cases, over a range of Mach numbers from to 5, 
with oj = 0-5, 0-75, 1, and 1-25, and also in a number of cases in which 
the viscosity -temperature law was given by a parametric form of Suther- 
land's law (226). An important conclusion from these calculations is 
that the variation of enthalpy with velocity, that is the function /*(w*), 
is virtually independent of which viscosity-temperature law is assumed. 
For this reason Crocco suggested that the relationship (260) or (261) 
could be assumed to hold in general, even when to 7^ 1, so that the 
solution of (251) alone remains to be determined. This does depend, 
though, upon the viscosity-temperature law, and the accurate value 
should be used. The solution of the equation is formally exactly as in 
the case cr = 1, since /**p* is known as a function of /*, and hence of 
w* and cr, and (251) takes the form 

GG"+k((j,u*) = 0. 

853502 II 
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This is similar to (253), and may be solved by an analogous procedure, 
so that we write * 



and substitute an approximation 6^ in the right-hand side to obtain 
a second approximation G 2 . Then either |($i+$ 2 ) or (G 1 G^ is the 
starting-point for the next iteration. 

This considerable simplification of the simultaneous equations, (251) 
and (252), was used by van Driest (1952), who carried out extensive 
calculations with a = 0-75, Mach numbers from to 20 and wall-to- 
free-stream temperature ratios, /*(0), from 0*25 to 6. 

We may deduce from Crocco 's work, and in particular from (260), 
that the relationship between heat transfer and wall temperature, that 
is between /*(0) and /*'(0)> is 

J*(0) = l-/* 



(0 = 
Now Q _l k 8T\ _kj8l\ I8u\ __J^ 

y -" ~ ~ ~ - ( ' 



where suffix zero on A x or A 2 denotes values at u* = 0. Accordingly, 
when there is zero heat transfer the wall temperature T w = T z is given 



or 



since it may be deduced from the results of Crocco, van Driest, and 
Brainerd and Emmons, and from the results of E. Pohlhausen for low 
speeds, that A a (<r, 0) =* Jo*. 

A similar method to that of Crocco was used by Hantzsche and 
Wendt (1940, 1942) to consider cases with and without heat transfer. 
Their results are mainly for either a> = 1 or a = 1, and, as they are less 
comprehensive than those of Crocco, they will not be considered further 
here. 
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We will, however, refer briefly to an approximate analysis by 
Monaghan (1949) for the relationship between enthalpy and velocity 
at arbitrary a close to unity. Monaghan makes the major assumption 
that enthalpy and velocity are related exactly as if each depended only 
upon local conditions, and deduces that 

(267) 



Now it is not difficult to express (260) or (261) in this form. We write 

(261) as m m 



^r 

J-Q 



and add and subtract the quantity 



whereupon we have 



J.Q J.Q J.Q 

By making use of (262) and (263) this may be simplified to 



A 20 ). (268) 

-*o J o J o 

Comparison with (267) shows that Monaghan 's solution replaces 
A 3 by ^U/UQ and A 2 A 20 by |(^/w ) 2 . 

Monaghan compares these values for two cases, a = O725 and a = 1-25. 
In each case he finds agreement within 1 per cent, for values of u/u 
less than 0-5, the error increasing to about 3 or 4 per cent, at ufu Q = 0-8, 
and the order of 10 per cent, at u/u =1-0. The formula given by 
Monaghan is exact when cr = 1, and in general we can remark that for 
values of a close to 1 the approximation is extremely good except per- 
haps at the very edge of the boundary layer. 

6. Summary of results for uniform wall temperature 

All the solutions considered so far in this chapter are for flows in 
which the wall is maintained at a uniform temperature. Taken together 
they provide a comprehensive set of results for values of or between 
0-5 and 2 and values of o> between and 1-25, both with and without 
heat transfer. Young (1949) has carried out a detailed analysis of these 
results, and has succeeded in developing simple formulae for skin- 
friction and heat transfer which agree with the accurately calculated 
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values to within about 1 per cent, over the physically important range 
of the parameters considered. 
Young expresses his results as a formula for the product c f *JR x , where 



and 



are the skin-friction coefficient and Reynolds number, respectively. 
He deduces from (258) that 

c f *lR x = 0-664 when o> = 1, (269) 

and shows from the results of Hantzsche and Wendt that 

c f <JR x ^ 0-664{l + 0-3(y~l)Jfgc7*p t >~ 1 >, (270) 

where the coefficient 0*3 and the factor <r* are approximate. These two 
formulae apply when there is zero heat transfer. For the case of heat 
transfer, Young obtains a rough solution of (251), which indicates that 
when or=l 

c f 4 x == 0-664{0-5 + 0-5(r fC /3 7 )+0-l(y-l)JfS}- (271) 

where the numerical coefficients within the brackets are approximate. 
By examination of (269) to (271) we conclude that for arbitrary o>, cr, 
and T w , the value of c f ^R x may be expressed as 



where a, b, c are roughly equal to 

a = 0-5, 6 = 0-5, c = 0-1. 

Young adjusts the values of these constants to get the best overall 
agreement with all of the available results, and thus finds that 

c f <IR x = 0-664{0-45+0-55(TJT )+0-09(y-l)Jfga^^- 1 >. (272) 

The accuracy of this equation is considerable, since it may be shown 
to agree with the exact values to within about 1 per cent, over the 
range of the parameters likely to arise in practice with air as the working 
fluid, as is illustrated in Fig. 1. 

It is now easy to derive a similar formula for the heat transfer rate. 
From (265) and (266) we may deduce that 

& T 1 

/i _ ^w l w im m \ * 
Vw; -- \ J -z~~ J -w) T~> 
PwUo A 10 
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and, since A 10 may be shown to be approximately o-*, this yields 

^ g"~ w a i f (273) 



To conclude the work on the compressible laminar boundary layer 
on a flat plate with uniform wall temperature we should make mention 
of a paper by Mack (1958) who has provided a considerable catalogue 
of results for the case of zero heat transfer over a range of Mach num- 
bers from 0-4 to 5, with viscosity and thermal conductivity prescribed 
accurately as functions of temperature from experimental data, and the 
Prandtl number therefore not absolutely constant. 

7. The solutions of Chapman and Rubesin 

All the preceding solutions apply only when the wall is at a uniform 
temperature. The problem becomes considerably more difficult when 
T w is allowed to vary with # in a prescribed manner. The first serious 
attempt to consider the high-speed boundary layer on a flat plate with 
non-uniform wall temperature was made by Chapman and Rubesin 
(1949), though some solutions for low-speed flow were given earlier by 
Fage and Falkner (1930). 

Chapman and Rubesin begin by writing the boundary -layer equations 
in the von Mises form, (59) and (60). Upon approximating 



which is equivalent to taking p oc T with the constant of proportionality 
left arbitrary, and then setting the pressure gradient equal to zero, 
these equations become 

du 



tf>r . A . 

(274) 

and opo M+Mopo0w ' (275) 



They then introduce non-dimensional quantities 

w * _ u r* - T - J LL* - ^ o* - p 

U* -- , J. - - = - , /i - , p -- , 

^o Zo J o ^o Po 

x* = ~, and 0* == 
L r 



and upon substituting into (274) and (275) this leads to 

(276) 
dx* 
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We note that the momentum equation (276) is now uncoupled from 
the thermal-energy equation, (277), so the solution of the former is 
simply the Blasius solution, 



where 77 = ifj*(x*)~*. 

This suggests a further change of independent variable from (#*,t/r*) to 
(x*, 77), whence (277) becomes 

'. (278) 



Now the general solution of this equation is the sum of a particular 
integral and the complementary function. For a particular integral we 
may choose a function of 77 alone, 

T* = N(r,), 

where N"+afN' = -\a( 7 -\)Ml(f"}*, 

which is very similar to Pohlhausen's equation. A solution is 



where #(77) is given by (159) and 

( r -l)Jf|0(0) 

. (279) 



The complementary function is a solution of 



and Chapman and Rubesin write this solution as 

T* = 7, 

where 7 = f a n a?* n ^(77), 

o 

and T w is the solution of the equation 

T' n +cfT n -2af'nY n - (280) 

with boundary conditions 

r n (o) = i, ^ r B (oo) = o. 

It follows that F(0) = f a n x* n . (281) 
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Hence, when the wall temperature is given by the sum of (279) and 
(281), that is T T 

37r = 3f+2>***. ( 282 ) 

-/o -*o o 

the temperature profile is given by 

Jr= l+*(y-lW0to)+|o** w r07). (283) 

In particular, the heat-transfer rate is given as 



, "T;(0). (284) 

^/ o 

It should be remarked that the above work holds for non-integral values 
of n. Taking a = 0-72, Chapman and Rubesin integrated equations 
(280) for several values of n, giving curves of the functions Y n (y]). In 
particular they gave the values of Y' n (0) which are shown in Table 16. 
These values may be interpolated for intermediate values of n, to yield 
a solution for any wall temperature distributions T*, which can be ade- 
quately approximated as a sum of arbitrary powers of x*. 

TABLE 16| 
Exact and approximate values o 



n (exact) (approx.) 

0-296 0-304 

1 0-489 0-490 

2 0-597 0-594 

3 0-678 0-671 

4 0-744 0-734 

5 0-799 0-787 
10 1-006 0-987 

f Reprinted from the Journal of the Aeronautical Sciences, vol. 16, p. 547, by per- 
mission of the Institute of the Aero/Space Sciences. 

Amongst the conclusions which Chapman and Rubesin draw from 
their results is that when the wall temperature is not uniform, heat- 
transfer rates cannot be estimated from uniform-temperature results. 
In this context, reference should be made to a range of solutions ob- 
tained by Baxter and Fliigge-Lotz (1957, 1958) including flows with 
zero pressure gradient and non -uniform wall temperature. These were 
obtained by means of a step-by-step solution on a digital computer. 
Amongst the cases considered are a number in which the wall tempera- 
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ture is uniform, T w = \T Z or %T Z upstream of a position x = I, with 
dTJdx constant (but not zero) downstream of this position. These 
authors have also considered cases in which the wall temperature is 
uniform except in the region 1 ^ xjl ^ 1-0256, where it increases con- 
tinuously by a maximum of 2 per cent, and then decreases symmetrically 
to its upstream value. 

8. Lighthill's analysis 

A method for calculating heat transfer in low-speed flow, due to 
Lighthill (1950), was described in Chapter 6, section 5. Now Lighthill 
has indicated that in equation (169), or equation (172), ju, and p only 
appear in the product pp. Accordingly, if this product is constant, the 
low-speed solution holds at all Mach numbers provided the appropriate 
value of the constant is taken. Thus we have, rewriting (173) 



Q w (x) - -0'5Ss( 

*i (285) 

For the case of a flat plate we may write 



so that (285) becomes x 

Q (x) = 0-339-^/^0)^-* f (x*x\)-* dT^xJ. (286) 

/* J 

o 

This is, strictly speaking, relevant to the case when T = in the main- 
stream and T = T at the wall. Lighthill shows that at high speeds 
it is also the value of the heat transfer for the case T = T t in the 
mainstream and T = T w at the wall provided T is defined as 

fp _ rn //y\ fp 
J. J-w^) 2 z> 

which, of course, reduces to the form (174) when the Mach number 
tends to zero. Accordingly, upon interpreting the Steiltjes integral in 
(286) we have x 

Q w (x) = -0-339* W^o)*k(0)-.r f +^ f ^^1 (287) 

P \ x I [ J ( x *!/ ) 

Lighthill has compared his results with those given by Chapman 
and Rubesin when the wall temperature is (282). Upon setting 
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(287) yields 



which may readily be shown to be equal to (284) with the expression 
in square brackets replacing -|F^(0). A comparison of results is given 
in Table 16 for the values of n considered by Chapman and Kubesin. 
We note that the error is nowhere as great as 3 per cent., indicating 
the considerable accuracy of Lighthill's method for flow with zero 
pressure gradient. 

Having thus established the accuracy of his method, Lighthill uses 
it to calculate the distribution of temperature which might arise in the 
high-speed flight of a projectile, when heat transfer to the body is 
balanced mainly by radiation from it. Now the rate at which heat is 
radiated at the surface of the body is given by 

Q v (x) = c^*)}*, (289) 

where a is the product of the Stefan-Boltzmann constant and the 
emissivity of the wall. Upon equating (287) and (289) we have 

f x ' \ 

Tf, - -0-339 A a */W^oUL (0) _. T+ ^ f ^A^ (290) 
otp \ x ] [ ' J (x* a?*)* J 

We note at this point that the only possible value for T w (0) is T z . If it 
had any other finite value then it would be necessary for the integral 
to be of order x~* as x -> 0, implying that T' w (x-^ = 0(x^ 1 ), which 
contradicts the assumption that T W (Q) is finite. Physically the boundary 
layer is extremely thin near x = 0, so that heat transfer is extremely 
effective in forcing the wall to take up its equilibrium temperature T a 
very quickly. If we neglect the term T W (Q) T Z , the remainder of the 
equation may be simplified by writing 



a 



--, 

and T* 3 = 0-678a* 

<xp.\ I ] 

whence it becomes z 

' < 291) 





with JF(0) 1. 

Lighthill showed that for small values of z a series solution 

F(z) = 1 l*461z+7-252z 2 46-46z 3 +332-9z 4 ... (292) 
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could be calculated, and that for large z there is an asymptotic expan- 
Slon F(z) ~ 0-84092-* 0-15243-* 0-01953-*... . (293) 

He also showed how the convergence of the series (292) could be ex- 
tended reasonably accurately, and so was able to link the expansions 
successfully. 

The principal conclusion to be drawn from Lighthill's analysis is that 
the temperature decreases extremely rapidly with increasing x near to 
the leading edge, with a more moderate rate of decrease further back. 
The actual temperature gradient predicted at x = is infinite, because 
of the assumption of zero boundary thickness. In practice this is 
countered by conduction of heat within the solid and (though to a 
considerably lesser extent) by the fact that the boundary -layer approxi- 
mation is inadequate there and the boundary-layer thickness should 
be regarded as small but non-zero. Even so one expects the wall 
temperature to fall from a maximum at the nose, so that melting may 
occur at high enough speeds with solidification further back. There is 
some evidence that this happens in the case of some stony meteorites. 

Before concluding this chapter, mention should perhaps be made of 
a much less exact but considerably simpler solution of the above 
problem by the author (Curie, 19586), based upon a method of calculat- 
ing compressible laminar boundary layers which will be considered 
later (Chapter 9, section 5). By this method the heat-transfer rate 
through a laminar boundary layer is given as 



Q w (x) = 0-332{!T.+r w (0)-2T w ( a! )}, (294) 



the Prandtl number being taken as unity. This replaces the more 
complex expression (290), involving an integral, and when equated to 
(289) leads after the appropriate algebra to the following simple alge- 
braic equation for F, namely 

ZF = 0-980(1 JF). (295) 

The solution of this equation is, of course, simple, and it follows fairly 
straightforwardly that the error in the predicted wall temperature in- 
creases from zero at the leading edge to a maximum of 18 per cent, far 
downstream, where the constant 0-8409 in (293) is approximated as 
0-9950. 



8 



THE COMPRESSIBLE LAMINAR BOUNDARY 
LAYER WITH ZERO HEAT TRANSFER 

IT has already been seen that the prob]em of the compressible laminar 
boundary layer with zero pressure gradient becomes extremely difficult 
unless suitable approximations are made concerning the value of the 
Prandtl number and the variation of viscosity with temperature. It is 
to be expected, therefore, that the same will be true when there is zero 
heat transfer but an arbitrary pressure gradient. Accordingly we con- 
sider first an approach to this problem in which the Prandtl number is 
taken to be unity, and the viscosity is assumed proportional to the 
absolute temperature. 

1. Howarth's method 

The equations of motion and continuity are given by (1) and (12), 
and with body forces neglected these become 

du . 8u dp , dl du\ .^^ 

pu -- \-pv- = -- ^-\ -- ha (296) 

^ dx^^ dy dx dy\dy) v ' 

and JL( pU )+JI-( P v) = 0. (297) 

Equation (297) may be satisfied by introducing a stream -function 
such that , , 



where suffix s denotes some standard datum position. Howarth (1948 a) 
assumes that the viscosity and temperature are related by 



and then attempts partially to reduce (296) to an incompressible form 
by making the transformations 







and <f,(x,y) = x (x,Y). (299) 
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dy 

He shows that u = - 

dY 

and that the momentum equation (296) takes the form 

8 X _8* x 8 X _ _ 

( ' 



with boundary conditions 

|=|= * r-0; -., . T~. 
This equation differs from the incompressible form only in the factor 

--** 

which multiplies the pressure gradient term, and in the altered scale 
normal to the wall. 

To proceed farther requires knowledge of the temperature field T(x, y), 
given by the thermal energy equation. When the Prandtl number is 
taken to be unity it follows from (19) that 

c p T-}-^u 2 = constant, 
and Howarth writes this equation in the form 



(302) 

\ u l/ 
It follows that (301) becomes 

> <308) 



so that (300) becomes an equation for x( x > Y) alone. 

Howarth then develops a momentum integral equation in the trans- 
formed variables, by integrating (300) from Y to Y = oc. This 
yields 

where oo 

8i= f /i-J^dF 

(305) 

and 8i= f^fl-^W 1 

J u i\ u il 

o * 

are boundary -layer thicknesses measured in the transformed coordinate, 
F, normal to the wall. He then seeks to develop a general method of 
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solution along the lines followed by Pohlhausen for incompressible flow. 
Thus w/% is expressed as a quartic polynomial in F/8', where S' is the 
boundary-layer thickness in terms of the ^-coordinate, and the coeffi- 
cients of the quartic are chosen to satisfy the corresponding boundary 
conditions to those satisfied in incompressible flow (Chapter 4, section 1). 
Howarth then finds that the relevant parameter is 



which may be compared with the incompressible form (97). Corre- 
sponding to (98), the momentum integral equation (304) yields an 
ordinary differential equation for A, which takes the form 



-% -t 

(306) 

where #(A), ^(A) are the same universal functions as in incompressible 
flow. The inaccuracies inherent in the incompressible form of Pohl- 
hausen's method are, of course, perpetuated in the compressible form. 
Howarth used his method to consider flow with a linearly retarded 
external velocity / \ 

WI = W Q 1- - (307) 

\ c l 

for various values of the Mach number M Q w / a o- He finds that the 
position of separation moves towards the leading edge as the Mach 
number is increased, a result later obtained by the more acceptable 
method of Stewartson (1949). 

2. Young's method 

An alternative method of generalizing Pohlhausen 's method was 
suggested by Young (1949). This method will not be discussed here, 
as it has since been improved and further generalized to flow with heat 
transfer at the wall (Luxton and Young, 1960). It is sufficient at this 
stage to point out that the method is valid for arbitrary cr and o>, where 
a) is the temperature -viscosity relationship index, and it is discussed 
in Chapter 9, section 6. 

3. The Ste war tson-Illing worth transformation 

The crux of Howarth 's method, which partially reduces the com- 
pressible boundary-layer equations to an incompressible form, is that 
the external velocity remains unchanged. By relaxing this condition 
Stewartson (1949) and Illingworth (1949) have independently shown 
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that it is possible to effect a complete transformation between a com- 
pressible boundary layer and an incompressible boundary layer with 
a different mainstream. The analysis given here follows Stewartson's 
approach. 

Stewartson assumes that a 1, so that for the case of zero heat 
transfer the temperature field is given by (302). He then transforms 
the ^-coordinate by writing 



Y = ' ~ f . dy, 

a Q V l J pQ 



(308) 





where suffix refers to any standard reference position, say at the 
leading edge in the mainstream. He also introduces a stream function 
ijj such that , , 



P u = p 4--, pv = -~Po4-~> (309) 

which satisfies (297) automatically. Then upon substituting into (296), 
transforming from coordinates (x, y) to (X, Y), making use of the assump- 
tion that CD = 1, i.e. /z oc T, and the temperature profile (307), Stewart- 
son finds that 



__ __ 

8Y 3 8Y8XdY^~8XdY* l dX' 



where X = | (~ ] * ? dx, (311) 

J \ a o/ 
o 

and U = a M^ = **%. (312) 

Since the relevant boundary conditions were shown by Stewartson to 
be 

\L = -?L- = when Y = ; ~- -> U l as F -> oo, 
0J: Si 

it follows that an exact correlation has been established between the 
compressible boundary layer in the (x,y) plane with external velocity 
%(#), and an associated incompressible boundary layer in the (X, Y) 
plane with external velocity U^X). It should be remarked that in 
deriving this same correlation Illingworth (1949) worked with the von 
Mises form of the boundary -layer equations. 

Since many methods are available for calculating the development 
of the incompressible laminar boundary layer, the use of this trans- 
formation enables compressible laminar boundary layers to be calcu- 
lated, subject to the restrictions made in the transformation, namely 
zero heat transfer, cr = 1, \L oc T. Stewartson himself used the trans- 
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formation simultaneously with Howarth's method (Chapter 4, section 3) 
to calculate the effect of compressibility on the separation of the laminar 
boundary layer with linearly retarded mainstream. His results are in 
qualitative agreement with the less accurate results obtained by Howarth, 
using the method of Chapter 8, section 1, and with the results given 
by Young's method in Chapter 8, section 2. Later, Rott and Crabtree 
(1952) suggested that for general calculations of compressible laminar 
boundary layers the transformation should be used with Thwaites's 
method, and this is done by many people. In view of the simplicity 
and overall accuracy of the latter method, this combination probably 
leads to good results in general, though it is possible that errors in the 
incompressible solution may be exaggerated at high Mach number, for 
example. 

This possibility is borne out by an accurate solution obtained by the 
Mathematics Division, National Physical Laboratory, for a particular 
external velocity distribution, 

% = U (IX/C), 

the velocity U Q at the leading edge being such that the Mach number 
is M = ^o/ a o = 4. It is known that when M = 0, boundary-layer 
separation occurs at x/c 0-120 (Howarth, 1938), and according to the 
precise numerical solution of N.P.L., separation occurs at x/c == 0-04(5) 
when M Q = 4. The corresponding solutions by Thwaites's method (as 
amended by Curie and Skan) and the Stewartson-Illingworth trans- 
formation are x/c = 0-123 when M Q = and x\c 0-06(7) when Jf = 4. 
We see that the error increases somewhat with Mach number in this 
particular case. 

Stewartson has also shown that, by applying the transformation, a 
considerable simplification can be effected even when there is heat 
transfer at the wall. In such a case, he writes the temperature as 

S, (313) 



where 8 -> as Y -> oo, but the boundary condition on S at Y = 
depends upon the temperature of the wall, and is S(X, 0) = only 
when there is zero heat transfer. It may then be shown fairly easily 
that 



_ 

d7 3 3Y , 

(314) 

, ty 38 di/j dS 3*8 ( 

and - ---- _ = 

3Y 3X 3X 3Y 3Y* 
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and these equations must be solved simultaneously for S and ijj. Some 
solutions of these equations were given by Cohen and Reshotko (1956 a), 
and were made the basis of a general method for calculating com- 
pressible laminar boundary layers with pressure gradient and heat 
transfer, which will be discussed in Chapter 9, sections 1 and 3. 

4. Rott's method 

Stewartson's transformation holds rigorously only when a == 1 and 
/x oc T. Rott (1953) showed how it may be generalized approximately 
to the case of arbitrary a and general temperature-viscosity law. The 
three approximations upon which Rott's method is based are as follows. 
Firstly the assumption is made that the wall temperature for which 
there will be zero heat transfer is given by 

where the recovery factor, r(cr), is given to a good accuracy as 

r(o-) ^= or*. 

This formula is exact when cr = 1, and has been found to be of good 
numerical accuracy for other values of a close to unity. Secondly, the 
expression m 



is accepted as an adequate approximate solution for the temperature 
profile. This satisfies the boundary conditions at the wall, (315), and 
at the edge of the boundary layer (where jP- 2\), and is exact when 
a = I. Thirdly it is assumed that the viscosity -temperature relation- 



a tr\ T (Y\ 

/*"i0v*7 L w\ J[ ') 

so that jit is proportional to T at constant x 9 but the variation with x 
is chosen so that the viscosity at the wall is given accurately, by 
Sutherland's law or by a carefully chosen power law. With these three 
approximations, the Stewartson transformation can be generalized, and 
the equations (308), (311), and (312) become respectively 

r=Kv n - 
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and U 



This transformation may then be used, in conjunction with say 
Thwaites's method, to provide a convenient approximate method for 
calculating the development of compressible laminar boundary layers 
with zero heat transfer, arbitrary a, and general viscosity-temperature 
law. 

5. The method of Oswatitsch and Weighardt 

One of the earliest methods developed for calculating compressible 
laminar boundary layers was that due to Oswatitsch and Weighardt 
(1943), their method being a plausible generalization of the method of 
Walz (1941), using the 'similar' solutions of the boundary-layer equa- 
tions. They write the momentum integral equation (40) in the alterna- 
tive form 

(316) 

v 



and further assume that the velocity profiles appropriate to specified 
values of the function 



are those given by the low-speed solutions of Falkner and Skan (1930), 
and used by Walz. It is then fairly straightforward to obtain BJSg and 
$z/ u i(d u ldy)w as functions of A* and M^ the dependence upon M l arising 
because of the variation of density with Mach number. Then (316) 
becomes 



and this is easily integrated by a step-by-step procedure. 

As no numerical examples were calculated by Oswatitsch and Weig- 
hardt it is difficult to estimate what accuracy is to be expected from 
the method. However, the assumptions made are somewhat arbitrary, 
so it is unlikely that the method is of particularly wide application. 

6. The work of Cope and Hartree 

A considerable attack on the problem of the compressible laminar 
boundary layer with pressure gradient and zero heat transfer was made 
by Cope and Hartree (1948), involving the use of a high-speed com- 
puting machine. These authors make first the following transformation 
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of variables * 

v- x ap 



l (^ (tody, 
2\v x) J p, J " 



.. TTr 

and W = --- = 



P 

1 dp! x \l 

-- 2 



with ?/r defined as in (168). They then derive the equations which u* t 
r, W must satisfy, and seek solutions of the form 



with two further expansions, 



and = ( 

Mo o 

introduced for analytic convenience. 

It is found that the zero-order terms, h Q , r , / , xo> ^o satisfy non- 
linear equations, and these have been integrated by Cope and Hartree 
for various values of Jf , assuming a = 0-715 and p, oc 3P 8/9 . Not sur- 
prisingly the solutions are barely dependent upon M , since they would 
be identical under the similar assumption that JJL oc T. Some of the 
first-order functions were also calculated approximately, but these 
exhibit a considerable dependence upon Jf . The wth-order equations 
(n ^ 1) are all linear, and depend upon the functions of lower order. 
In principle, therefore, they may be calculated successively in order to 
obtain a complete solution. 

Cope and Hartree have also given some consideration to the possi- 
bility of devising a solution along the lines of Pohlhausen's method 
for incompressible flow. They assume that the Prandtl number is close 
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to unity, so that an adequate solution of the thermal energy equation 
is given by (302) when there is zero heat transfer. The momentum 
integral equation (40) is first written in the alternative form 



and the velocity profile approximated as 



where 8 is a measure of the boundary -layer thickness, and the function 
/ is chosen to satisfy a number of the boundary conditions on the 
velocity u and its derivatives. In the absence of a pressure gradient 
Cope and Hartree take 

f(rj) = sin^, 

and obtain results for the skin-friction which agree well with the calcula- 
tions of Brainerd and Emmons. 

When there is a pressure gradient the algebra becomes considerably 
more complicated, and Cope and Hartree concluded that the Pohl- 
hausen method did not offer much promise of progress, a conclusion 
which was perhaps a little pessimistic. 

7. The work of Illingworth, Frankl, and Gruschwitz 

We might perhaps mention at this point some of the other investiga- 
tions of the Pohlhausen-type which have been made by various authors. 
Illingworth (1946) has considered the case of zero heat transfer and 
arbitrary pressure gradient, with the assumption that viscosity ^ and 
thermal conductivity k are constant. The velocity was defined by 

u 



and the temperature (or density) by 



Pl Pl \ Pl) 

where x = S'/S is the ratio of thermal boundary-layer thickness to 
velocity boundary-layer thickness. The functions / and g are chosen 
to satisfy the usual boundary conditions, derived from the differential 
forms of the equations. The momentum integral equation and the 
thermal energy integral equation then yield two relationships for the 
three unknowns p w /pi, x> an( i & 2 I u i(d* u ly 2 )w> a third relationship being 
given by the boundary condition which follows from setting y = in 
the thermal energy equation (18). 
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Illingworth has carried out detailed calculations for flow with zero 
pressure gradient, using a number of possible forms of the functions / 
and g. In particular he considered cases when / and g were of the same 
order in 77, namely quadratic, cubic, and quartic, and one case in which 
/ was quintic and g sextic. Further, he also considered trigonometric 
forms, in which /, g are 



f - i / ^Y 

f = sinfTTT], g = lsin~-M 



Upon comparing his results with the accurate calculations of Emmons 
and Brainerd (1942), Illingworth concludes that the trigonometric and 
cubic forms appear to give the most accurate values. He notes, how- 
ever, that the particular form of the function g is critical and that the 
predictions of wall temperature leave much to be desired in all the cases 
considered. 

Illingworth 's work may be regarded as a generalization of the earlier 
work of Frankl (1934), who made the further approximation that the 
viscous and thermal boundary layers have equal thickness (x = 1), 
although the Prandtl number was not assumed to be unity. Various 
possible forms of the functions /( 77), g(rj) were considered, and in spite 
of the apparent additional assumption (as compared with the later 
work of Illingworth) Frankl's results for zero pressure gradient, with 
/ and g both quartics, were better than any of those obtained by 
Illingworth. This is perhaps indicative of the uncertainty which besets 
attempts to make direct extensions of Pohlhausen's method without 
first making a suitable transformation of coordinates and/or variables. 

A slightly different approach was adopted by Gruschwitz (1950) who 
introduced a new boundary -layer thickness 



CxY 

= ( dy, 
J ft 



A 

and expressed u\u as a quartic in 



PI 
o 

Further, instead of expressing p/p l in terms of 77, he chose to write the 
product {I (U/UI)}(PI/P) as a quintic in 77, satisfying appropriate 
boundary conditions. Substitution into the momentum integral equa- 
tion and the kinetic-energy integral equation leads, after some con- 
siderable calculation, to a first-order non-linear ordinary differential 
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equation for the momentum thickness 8 2 in terms of the known % and 
M 1 and the unknown T w /T l9 and an explicit expression for T W /T 1 in 
terms of M l and S 2 . These two equations involve universal functions, 
analogous to those arising in Pohlhausen's method, which are dependent 
only upon the Prandtl number a- and the shape parameter 

$ 2 T 
\ _ d 2 f J-w 

7 17 T" 
v i 2 i 

It should be remarked that the method does not make use of the 
thermal -energy integral equation, and accordingly is suspect unless 
a = 1, when the kinetic -energy integral equation is identical with it. 

8. The investigations of Gadd 

We conclude this chapter with a brief discussion of some work of 
Gadd (1953 a), which investigates approximately the effects on the 
separation position of variations in Prandtl number a and the index 
to in the viscosity-temperature law (227). 

He begins by considering an artificial case in which a = 0. In such 
circumstances it is clear, by inspection of ( 1 9) , that temperature gradients 
can remain finite only if 

%r) = 0> 

3T 

so that Lt = 

dy 

when the heat transfer at the wall is zero, and thus 

T(x,y) = T&). 

With this solution of the thermal-energy equation, Gadd shows how 
the momentum equation may be transformed into an incompressible 
form by means of a transformation analogous to that of Stewartson 
and Illingworth. This equation may then be solved approximately by 
any convenient method, and Gadd gives some solutions (at various 
Mach numbers) by Thwaites's method for the case when a> 1 and 
the external velocity is of the form (307). Comparisons with solutions 
for the case to a = I suggest that the distance to separation increases 
with decreasing a, and the solutions for a = provide an upper bound 
to this variation. 

Gadd then considers approximate solutions for cases in which either 
co or a is equal to unity and the other does not differ much therefrom. 
In each case, he writes the equations in the Crocco form, with inde- 
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pendent variables x, u, and then further transforms to variables x and 
/ = ufa. 

When a = 1, a> = 1, the thermal energy equation has the simple 
solution (228) for temperature in terms of velocity. The solution of the 
momentum equation was then obtained by an approximate method. 
When a) = 1, a ^ 1, the process is slightly more difficult, as it is 
necessary to obtain first an approximate solution for the temperature 
distribution, and Gadd did this by means of a small perturbation about 
the solution (228) valid for small values of jl <r|. 

The conclusions to be drawn from Gadd's calculations may be illus- 
trated as follows. When a) 1, the proportionate increase in the 
distance to separation when a = 0-72 (as compared with a = 1) is 
about 8 per cent, when the Mach number at the leading edge is M = 4, 
and about 25 per cent, when M =10, these results being for the case 
of a linearly retarded external velocity. Further, when <r = 1, the 
proportionate increase in the distance to separation when w = | (as 
compared with o> = 1) is about 2 per cent, when M Q = 4 and about 
6 per cent, when M Q = 10. Gadd has also derived results by direct 
application of the approximate method of Young. At the lower Mach 
number, M = 4, the results are almost identical with those quoted 
above, but at Jf = 10 Young's method predicts proportionate changes 
in separation distance which are in each case approximately half those 
derived by Gadd's method. 

We deduce from these results that, when there is zero heat transfer 
at any rate, the errors induced by assuming o> = or = I will be no more 
than about 10 per cent, when M Q = 4 and 30 per cent, when M = 10, 
15 per cent, if Young's values are to be preferred, in practical applica- 
tions with air as the constituent gas. These conclusions are broadly 
similar to those which would be reached if the method of Rott were 
used in the calculations. 



THE COMPRESSIBLE LAMINAR BOUNDARY 

LAYER WITH PRESSURE GRADIENT AND 

HEAT TRANSFER 

As far as the formal boundary -layer problem is concerned, the case 
with both pressure gradient and heat transfer is by far the most difficult 
to deal with theoretically. The momentum and energy equations are 
simultaneous non-linear partial differential equations for which no 
simple approach is known. The problem becomes even more difficult 
if the wall temperature is not uniform (or very nearly so), a situation 
which could arise in practice at very high speeds. The mathematical 
and numerical difficulties encountered in solving these equations are 
such that, even with the assistance of modern high-speed computing 
machines, very few precise numerical solutions for special cases have 
been attempted, and some of these are of doubtful accuracy. 

We shall begin this chapter by discussing briefly these numerical 
solutions, continuing with some of the approximate methods developed 
for this very general problem. In incompressible flow, where the pressure 
gradient alone is sufficient to prescribe the problem, there are now 
about eight really accurate solutions for particular cases, half of which 
are for flows in which the pressure gradient is everywhere adverse. 
These eight solutions are only just sufficient to enable us to make a 
reliable assessment of the accuracy and value of various approximate 
methods. In compressible flow things are further complicated by the 
presence of additional parameters such as Mach number and wall 
temperature. Clearly at least six or eight accurate solutions are re- 
quired for several values of these parameters before we can assess the 
usefulness of the approximate methods with any certainty, and this 
point will accordingly not be considered as thoroughly as in the incom- 
pressible case. 

1. Accurate numerical solutions for special cases 

A range of accurate solutions are, of course, available for the case of 
zero pressure gradient, covering various values of the Prandtl number 
a, several possible viscosity-temperature laws, and a few values of the 
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wall-temperature. These were discussed and summarized in Chapter 7, 
and will not be further considered here. 

For flows with pressure gradient, the earliest accurate solutions (as 
far as the author is aware) were obtained by the Mathematics Division, 
National Physical Laboratory, and are first referred to in a paper by 
Gadd (1952). Three particular cases were considered, the assumptions 
that a = 1 and // oc T being made in each case. All were such that the 
wall temperature was uniform and the pressure gradient everywhere 
adverse, so that flow near to the leading edge was like that on a flat 
plate with zero pressure gradient. Further details of the three cases 
are as follows: (i) In the first case a uniform adverse velocity gradient 
was assumed, so that the external velocity was given by (307), the 
Mach number at the leading edge was M ( = 4 and zero heat transfer 
was considered, (ii) The data for this case differed from case (i) only 
in that the wall was taken at a temperature TQ, which is the external 
stream temperature at the leading edge. Since TJT is equal, by (244) 
to l + |(y l)Jfg, upon putting y = 1-4 we see that the ratio of wall 
temperature T w (= T ) to its zero-heat-transfer value T z is ^. In other 
words the wall has been considerably cooled, (iii) In the third case, 
a uniform adverse pressure gradient was assumed, so that 

p = , 

The Mach number at the leading edge was M = 2, and the wall tem- 
perature was T Q (defined as above), equal to five-ninths of its zero- 
heat-transfer value; the cooling at this lower Mach number is less 
severe. 

For each of these cases detailed calculations were made of both the 
viscous and thermal boundary layers in the region between the leading 
edge and the position of zero skin-friction, and some of these results 
have been reported in a paper by Curie (1959). For further details, 
reference may be made to the Mathematics Division, National Physical 
Laboratory. 

Some time later a series of 'similar' solutions were derived by Cohen 
and Reshotko (1956 a), analogous to those found for incompressible 
flow by Falkner and Skan and by Squire. Following Li and Nagamatsu 
(1953), Cohen and Reshotko look for solutions of (314) in which the 
external velocity in the transformed plane is equal to 

UJiX) = CX m , 
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and introduce new variables /, 77, where 



and 



Then equations (314) reduce to the form 

j jj __ ^qpjw I ^ 317 j 

S"+ofS' = J 

Here primes denote derivatives with respect to 77, and the boundary 
conditions are 



= /'(O) = o, 

/'->!, ->0 as 

Equations (317) and (318) were integrated numerically for a series of 
values of the parameters m and S w , corresponding respectively to 
different pressure gradients and wall temperatures, and the results are 
tabulated and plotted in their paper. Although the same reservations 
must be observed in applying these results to non-similar flows as were 
relevant to the Falkner-Skan solutions, there is no doubt that the 
solutions of Cohen and Reshotko are extremely valuable, none the less. 
Subsequent to the above solutions Baxter and Flugge-Lotz (1957, 
1958) have reported a series of some sixty examples which were calcu- 
lated on a digital computer by means of a step-by-step solution, starting 
from Crocco's form of the equations, and have given detailed results 
for the skin-friction and heat transfer. As an illustration of the type 
of example considered by these authors, we mention four cases in which 
the external velocity is constant, equal to U Q , when x ^ I, and varies 
when x ^ I according to the law 



\AJJJ rQ 1 7 

dx I \l 

The four cases differ only in that T w = \T Z in two of them, with T w = 3T S 
in the other two, and M Q = | in two cases, with Jf = 3 in the other 
two. In all the examples considered the Prandtl number was a = 0-72, 
and Sutherland's law was used for the viscosity- temperature relation- 
ship. For details of the method used and for the numerical results, 
reference may be made to the two papers by these authors. 
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2. Kalikhman's method 

Turning now to approximate methods for calculating the compressible 
laminar boundary layer with heat transfer and pressure gradient, we 
begin with a generalization of the Pohlhausen method, due to Kalikhman 
(1946). The method is valid for arbitrary constant Prandtl number, and 
assumes that the viscosity-temperature relationship is 



Kalikhman introduces a new independent variable 77 in place of?/, where 



Ps 

o 



and suffix s denotes a standard reference condition, which he takes to 
be the free-stream stagnation value. He also introduces new dependent 
variables v, E*, defined by 



where T H is the total temperature, as defined in equation (48). The 
continuity, momentum, and thermal energy equations (1), (12), and 
(19) then take the form 

^ + ?=0, (319) 

fix dr) 

8u , .du p, du, , d(p(T\<du) /OOAX 

u -- \-v = ^Ui -+v n { 7^ 1, (320) 

rJz drj p 1 dx^ dr)\p 8 \Tj drjl V } 

, dT* , _aT* dTf , 8 {pIT\> d \-. , /I \-'\\ /001 . 
and u -- \-v - = ^ - -1 + 1/_ (H=H M + -- 1 (321) 
dx ^ dri l dx ^ S drj\p h \Tj drj[ ^\<r / J) 

It will be observed that the various transformations have gone some 
way towards reducing these equations to an incompressible form, in 
much the same way as the related transformation of Howarth (1948 a). 
Kalikhman then integrates these equations with respect to 17, and 
obtains the results 



and 



* = -, (322) 

\ J S / Ps U l 
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where m 



o* 

8 



r u L u 

= 1 -- 
J Uj\ U 



C u I T* 

I "* I i -* 

= d 1 "? 5 

J %\ ^1 



and 8*, 8*, and 8| are analogous to the displacement thickness, momen- 
tum thickness, and enthalpy thickness of the boundary layer as pre- 
viously defined. Following the idea of Pohlhausen, he then makes the 
approximations 



A 

where the ten coefficients A n , B n are chosen to satisfy the usual boundary 
conditions, namely /QX _ 

and 

u _ T* du_d 2 ^__^^. __ - - __ __ 



together with the two relationships obtained by setting r] = in the 
momentum and energy equations, (320) and (321). Then (322) and 
(323) yield two simultaneous ordinary differential equations relating 
the boundary -layer thicknesses, 8 and A, which appear in the para- 
meters 

A = 



* 12+A 

and A, = - ?- -- 



where a = (1 < 

No numerical results were given by Kalikhman, so the accuracy of the 
method is uncertain. It may be remarked, however, that the partial 
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transformation of the equations to incompressible form is a point in 
favour of the method, although the arbitrary choice of the order of the 
polynomials in (324) and (325) is unfortunately somewhat restrictive; 
an examination of which order would yield best results (by comparison 
with accurate solutions) would be extremely lengthy. 

3. The method of Cohen and Reshotko 

Mention has already been made of some 'similar' solutions for com- 
pressible laminar boundary layers with pressure gradient and heat 
transfer. These solutions were calculated by Cohen and Reshotko 
(1956 a) and were subsequently used (1956 6) as the basis of an approxi- 
mate method for general use. 

These authors begin by applying the Stewartson-Illingworth trans- 
formation, in a form which neglects the factor v appearing in the stream 
function as defined by Stewartson. Thus they write 

dib ddt 

pu = Po -, f=-p*te, 

and introduce the new independent variables, 

fViY 3y ~ 1)/(y - 1) , v i f p * 
= M dx, F = ~<%, 

J W a oJ Po 



which differ from (308) and (311) only as regards the factor v. Then 
with fl / o / 

TJ y v d r 

W dX' 

the continuity, momentum, and thermal-energy equations may be 
written as ^ TT ^ T7 



(327) 

. TT dS v dS d*S 

and U + V = v n , (328) 



where 8 is defined by equation (313). The boundary conditions are 

U(X, 0) = V(X, 0) = 0, S w = S(X, 0) = ?_ 1 \ 



U -> C/i, S -> 0, as F -> oo / 
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Upon defining boundary-layer thicknesses, measured in the F-coordi- 
nate, as follows, ^ 



and 3* = 

o" 

equations (326) to (328) may be integrated with respect to Y from 
Y = to Y = oo to yield 

^+^(8*+28J) = ^/|^j , (330) 

and ^J+^SJ^ -(H\ (331) 



Cohen and Reshotko now introduce a number of non-dimensional 
quantities, analogous to those introduced by Thwaites (1949) for in- 
compressible flow, namely 

I - B *( 8U \ 
1 rfl"flv) ' 
U l\ dY lw 



and ^* = S*/S*. (335) 

Upon substitution into (330) it is found that the following result may 
be obtained, namely 



= L, (336) 

which is analogous to the low-speed form (103). 

Cohen and Reshotko make the assumption at this stage that L, H*, 
and I are all functions of n and S w only, these functions being derived 
from the 'similar' solutions which they had calculated. This assump- 
tion is analogous to that made by Walz (1941). One consequence of 
the assumption is that equations (330) and (331) cannot be satisfied 



9.3 WITH PRESSURE GRADIENT AND HEAT TRANSFER 127 

simultaneously, and Cohen and Reshotko choose to satisfy (330) whilst 
rejecting (331). One expects intuitively, therefore, that in certain 
circumstances the method may give the temperature field with only 
low accuracy, whilst predicting the velocity field more successfully. 
When L(n,S w ) is thus prescribed, equation (336) may be integrated 
numerically as in the incompressible case. When the wall temperature 
is uniform it is often possible to approximate L as a linear function 
of n, and an analytic solution of (336) then becomes possible. Having 
calculated 8*, 8* follows from (335), (dUldT) w from (332) and (BS/dY) w 
from (334). It is then a simple matter to derive from these results the 
true boundary -layer thicknesses, the skin-friction, and the heat-transfer 
as measured in the physical plane. The relevant formulae are 



'du 



(dT\ T, T r 
and \~a~} ^T^/TT ~~* 



The functions I, H* 9 L, and r are plotted in Cohen and Reshotko 's paper 
as functions of n for various constant values ofS w . 

Comparisons have been made with some of the accurate solutions 
discussed in section 1, and the agreement is quite good. For example, 
for the case of a uniform adverse external velocity gradient at a Mach 
number of 4, with the wall cooled to a temperature T w = jfcT z , separa- 
tion is predicted when the external velocity has decreased by about 
17| per cent., as compared with the accurate value of 22J per cent. 
This comparison indicates fairly well the overall accuracy of the method 
for this particular example. 

4. Monaghan's method 

Important simplifications and improvements to the method of Cohen 
and Reshotko have been suggested by Monaghan (1960). In the course 
of a more general study of the compressible laminar boundary layer, 
Monaghan found that by careful choice of parameters some of Cohen 
and Reshotko 's curves for the function Z, at different values of S w , may 
be collapsed onto a single curve. In particular it appears that in regions 
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of favourable pressure gradient the function l(m,8 w ), where 



is almost independent of S w , whilst in regions of adverse pressure the 
functions Z{(m/w gep ), S w } fall on a single curve, m gep being the value of 
m at which I = 0. A similar collapse onto a single curve is found for 
a function which is related to the r of Cohen and Reshotko's method. 
Monaghan further generalizes the transformation of coordinates used 
by Cohen and Reshotko, and writes 



?T? dx 



His final equations take the form 



8 2 = 0.664^*, 
\ %/ 



^ , 

and Q w = 



, A . , , T w X 

where m= ~ 0-441 -g - , 

T z M dx 

and the functions l(m), h(m) are given to a good approximation as 



and, when m > 0, h = 0-225 

and Z is given as a function of m/m gep by means of a graph in Monaghan 's 
paper. 

It remains to define g lt g 2 , and m gep . To a first approximation g : and 
<7 2 are given by 

7 = 3+2' and = -?. 



The value of m sep is determined by means of a graph relating w gep /m 
to T W /T Z , where m is theoretically equal to 0-0681. However, Monaghan 
suggests an empirical improvement in which m is taken to equal the 
value of m at which the same external velocity distribution would lead 
to separation in incompressible flow. If we may judge by experience 
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in incompressible laminar boundary-layer theory, it should be taken 
to equal about 0-090. 

Comparisons with the exact solution for the case considered earlier 
(M = 4, T w fT z = |p uniform adverse velocity gradient) indicate that, 
with all the empirical improvements noted above, Monaghan's method 
is remarkably accurate. 

5. Curie's method 

The two methods above, namely those of Cohen and Reshotko and 
of Monaghan, may be regarded as possible generalizations of the method 
of Thwaites (1949), or more precisely that of Walz (1941). An alterna- 
tive method of doing this was presented by Curie (1959) in a paper 
giving a generalization of Thwaites 's method starting from the Howarth 
transformation, as defined by equations (298) and (299). Upon formally 
applying this transformation we find (following Howarth's analysis) 



20} 



and this may be formally integrated with respect to F, from F = to 
Y = oo, to yield 

N 

| p + ^tti 8i+(2-iyJff)8i+ f (~-l) rf7 = vj^l , (337) 

ax j \i l ] j \v*/w 

where 8^ and 8 2 are defined by equation (305). Further progress depends 
upon a knowledge of the temperature profile. 

Now it is well known that when a = 1 the total temperature is 
constant if there is zero heat transfer, and is a linear function of the 
velocity u if there is zero pressure gradient. For the purpose of obtain- 
ing a straightforward solution of (337) for the velocity boundary layer 
(but not for calculating the heat transfer, for example) Curie suggests 
writing / 7 2\ 

tji n~i\imi>rF\\T?nnl w \ ^Q^tt\ 

JL ff = J[,,.-\~[j. z -*-w)~ (Jtl.-/il ~~o|> \OOOj 

where the last term should be zero when either the pressure gradient 
or the heat transfer is zero. The value 

T e T w 

XTi = 1 , z ck(t "^ $> 

where 
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is suggested for use in regions of retarded flow. The temperature T , 
velocity u , and Mach number M Q are values at the leading edge if the 
velocity there is non-zero, otherwise they are values at the pressure 
minimum. The analysis cannot be directly applied to regions of ac- 
celerated flow, where Curie indicates that Cohen and Reshotko's method 
should give good accuracy. 

With the approximation (338) it is easy to show that 



J (|- 



dY = r_i8i+(ff+i(y- l)Jff)8i, 
whence (337) becomes 



(339) 

\ GJ -/w 

Curie now introduces non-dimensional parameters 

H' = S(IS' 2 , 
V = **( 



whereupon (339) becomes 

(340) 



where L' = 2l f ~2X'(H f +2). 

By considering the consequences of determining relationships between 
L', H', Z', and A' by means of velocity profiles prescribed in the F-plaiie, 
Curie concludes that the functions L'(\'}, H'(\ f ], Z'(A') are largely inde- 
pendent of wall temperature and Mach number, and are given approxi- 
mately by their incompressible values; the values appropriate to 
Thwaites's method are accordingly recommended. Since we can now 
use the approximation to L'(A') analogous to the Thwaites value (106) 

wehave i'^O-46-eA', 

and (340) is integrable in the form of a simple quadrature. Upon trans- 
forming from the F-plane into the physical plane, the final equations 
of the method are x 

j^ 
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where 



x 

, = expf2 f (J + 2+^-pf?)^ dx\ (342) 

( J w l u i j 



(343) 
(344) 



The method is perfectly straightforward to apply, and comparisons 
with the three accurate solutions calculated by the Mathematics Divi- 
sion, National Physical Laboratory, indicate that the accuracy bears 
reasonable relationship to the work involved. 

The method as discussed above is valid for or = 1 , ft oc T, and 
uniform T ir . In two later papers Curie (19586, 19616) has removed 
these restrictions. For full details reference may be made to these 
papers. Briefly, the removal of these restrictions leaves the general 
procedure unaltered, but amends some of the details, so that (for 
example) the coefficients of TJT^ and Jff in (342) and (345) are pre- 
scribed functions of a, and an additional factor to account approxi- 
mately for the generalized viscosity-temperature relationship appears 
in (344) and in the definition of the function G v 

6. The method of Luxton and Young 

We now consider a method which is essentially a generalization of 
the Pohlhausen procedure, due to Luxton and Young (1960). It was 
originally developed by Young (1949) for the case of zero heat transfer, 
but was later generalized and improved by these two authors for use 
when there is heat transfer, and is applicable for arbitrary values of 
the Prandtl number a and with the viscosity-temperature relationship 
given by (227) for arbitrary a>. 

The method assumes that the velocity profile is of the form 

u Y 



where T == 
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and a, b, c, d are chosen to satisfy the usual boundary conditions. Thus 
6 = A, c = 2 A, d = 1 A, 



SS2 

and A is given by A = u\ . (346) 



It follows that - = . (347) 

6% 8 



The form parameter H = Sj/Sg 

is not chosen from the assumed velocity profile, but is taken to be 

H = 2 .59 + fc 2 A + i( y -l)^i, (348) 



which is approximately correct in the case of zero pressure gradient 
(when cr is close to unity), and agrees fairly well with values derived 
from Cohen and Reshotko's solutions provided & 2 is a suitably chosen 
function of T W /T Z . Further the ratio of total boundary -layer thickness 
to momentum thickness is not derived from the velocity profile, but is 
given by 

_- = / = 9-072[o-45+0-55^-' + 0-09(y l)Jffa*l "(l+J^A), (349) 

S 2 L T l \ 

which again is correct in the case of zero pressure gradient, and gives 
approximately the correct effect of pressure gradient when compared 
with Cohen and Reshotko's values provided lc is correctly chosen as 
a function of TJT Z . 

Upon substituting from (346) to (349) into the momentum integral 
equation, written first in the form 



it is found that 

which may be formally integrated to give S 2 as a function of x. Luxton 
and Young suggest a numerical method of integration, but upon intro- 
ducing a preliminary quadrature to give 

G a = exp 2 



J {H+i- 
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it may be integrated as 



= !_ f 
rfej 



Having obtained 8 2 , 8 X follows from (348) and then r w from (347), making 
use of (349). 

It will be recognized that this method contains elements of both 
Pohlhausen's method and of the similar solutions. At first sight the 
relative contributions might appear to be rather arbitrary. On the 
other hand, Pohlhausen's method, applied formally in incompressible 
flow, usually yields values of the skin -friction which are too high, with 
too late a separation, the value of the pressure gradient parameter at 
separation being A = 0-157. Methods based on the similar solutions, 
such as that due to Walz, usually err in the opposite direction, with 
A 0-068 at separation. It is interesting to note that the method 
of Luxton and Young yields A = 0-090 at separation in incom- 
pressible flow, which was shown by Curie and Skan (1957) to be the 
optimum value. Comparisons with some of the accurate solutions of 
section 1 suggest that at a Mach number of 4 the method still yields 
results of good accuracy. 

7. The method of Foots 

A method has recently been given by Foots (1960) which is theo- 
retically more acceptable than any of the methods given in this chapter, 
but becomes much more complicated and lengthy in practice. The 
method begins with the Stewartson transformation, so that upon further 
introducing transformed velocities 

//__ <W_ v - _^ 
BY' dX' 

equations (314) become 



+ F 

where X, Y are defined by (308) and (311), and the boundary conditions 
are given by (329). By formally integrating (351) and (352) with respect 
to Y from Y to Y = oo, and by multiplying (351) by U and then 
integrating, three equations are obtained, analogous to the momentum, 
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kinetic-energy, and thermal-energy integral equations. These may be 
written as 



and Ul ~ +2U '^ == - 2 ^ (355) 

where ^ 



CO 

and A* - J 8 dY, 



the definition of 8* being different from that used by Cohen and 
Beshotko. 

Foots now expresses U/U-^ and S as quartics in F/8*, where 8* is 
a boundary-layer thickness in the transformed plane. These quartics 
are chosen to satisfy the conditions 

U = 0, 8 = 8 W at F = 0, 
and 

-1 at 7-3* 
' [ZL "~ ' 

but, following the method devised by Tani (1954) for incompressible 
flow, the conditions obtained by setting Y = in (351) and (352) are 



dY ~~ 31* " dY 



Yi y\ 3 
ll I 

8*\ S*/ 
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not satisfied, a free coefficient being left in each of the quartics. Thus, 
for the case when S w = 1, 

U Y 2 / Y Y 2 \ Y 
_L_ f___/6~8 4-S 14-tf 

Di~8* a \ * S**/ + 8 
y2 / y 

and ^ 1 -H 6 - 8 J+ 3 

Equations (353) to (355) must be solved simultaneously for 8*, a, b. 
The problem has thus been reduced to that of solving three simultaneous 
ordinary differential equations. 

Foots has expressed the various factors in (353) to (355) in terms of 
8*, a, and b (for the case when S w = 1), and has gone on to show how 
the simultaneous integration might be done, taking as an example the 
case when the external velocity in the transformed plane is 

U, = U (l-X), 

the wall temperature, of course, being T w = 2T Z . He finds that separ- 
ation of the boundary layer occurs at 

X 8 = 0-075, 

and this value agrees with that obtained by expanding a power series 
solution of the partial differential equations, valid at small values of X, 
and continuing the solution by a numerical step-by-step method. 

This method provides simultaneously solutions for skin -friction, heat 
transfer, and indeed any desired boundary -layer property. It is, there- 
fore, of wider application than the methods previously considered, 
where the thermal energy equation is sacrificed in an attempt to obtain 
a quicker approximate solution for the viscous boundary layer. If the 
computing facilities are available the method should yield good accuracy, 
even in cases of non-uniform wall-temperature. 

8. The methods of Lilley and Illingworth 

Two rather similar methods of calculating the skin-friction and the 
heat-transfer rate in a compressible laminar boundary layer have been 
given by Illingworth (1954) and Lilley (1959). In what follows a brief 
summary of Lilley 's approach will be given, together with some dis- 
cussion of the numerical examples considered by both these authors. 

The method adopted by these two authors is designed to extend the 
analysis of Lighthill (1950 a) to the case of compressible flow with both 
heat transfer and pressure gradient. Unlike Illingworth, Lilley applies 
first the Stewartson-Illingworth transformation and then, following 
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Lighthill, transforms the equations further to the von Mises form, in 
which the independent variables are X and i/j, and the dependent 
variables are 8 and Z, where 



and S, X } $, are defined respectively by (313), (311), and (309). The 
momentum and energy equations then take the form 



,356, 
and - 



and the assumptions or = 1, ^ oc T (not made in Lilley's analysis) are 
made here for simplicity. Lilley rewrites (356) in the form 

] 8(i,,t)dU*(r,)) = vU~(z+ ] 8( 

/ r V 

and notes that at both if; = and ijj = ao this may be written as 

X 

where G(X, i/j) = Z+ J 8 fa i/j) dUKrj). (359) 

o 

Lilley assumes that (358) is an adequate approximation right across 
the boundary layer, and solves it subject to the boundary conditions 

0(X 9 oo) = 0, #(0, 0) = 0, (360) 

and for small iff 

(361) 





and since, by (313), S u , = |f - 1, 

I z 

this further simplifies to 

x 

G(X, 0) = E7J(+0)- f/ 2 (Z, t)+ j ^M dUKr,). (362) 



Lilley indicates that equation (358) may be integrated, subject to the 
boundary conditions (360) and (362), by a method analogous to that 
used by Lighthill, the term U 2 (X, ^) in (362) being first replaced by 
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where ^c is the viscosity calculated at the reference conditions. His 
solution is 



(363) 

Lilley further transforms back into the physical plane (for which see 
his paper), and also shows how (363) may be approximately inverted 
to yield r w as an explicit function of U^X) and T W (K) ; again no details 
will be given here. It is easily seen that (363) reduces identically to 
the corresponding expression in Lighthill's paper when incompressible 
flow is considered. 

Lilley deals with equation (357) in a similar manner, the boundary 
conditions being OQ) = g(0> ^ = ^ 



and B=-l+ 

J z c p- L z\P'0 T 

where Q* = JkJ_J -. 

The solution follows Lighthill 's solution for the heat transfer, and is 




(364) 

Lilley indicates how the solutions given above may be improved by 
taking & ^ I, by introducing a factor C(X), following Rott (1953), to 
account approximately for the viscosity-temperature law, and by adapt- 
ing the technique of Spalding (1958) so as to extend the range of con- 
ditions over which the theory is accurate. 

Both Lilley and Illingworth (by his basically similar method) have 
considered in detail the numerical results for a number of special cases, 
and obtain fair agreement with the accurate solutions. In particular, 
Lilley has considered the three cases considered by the Mathematics Divi- 
sion, National Physical Laboratory, and obtains an accuracy which, 
for the velocity field at any rate, is comparable with that given by 
the other methods of this chapter. No comparisons were made for heat- 
transfer rates. Illingworth has not made any calculations for flows 
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with both heat transfer and pressure gradient at high Mach numbers. 
He has, however, made some comparisons with the low-speed 'similar' 
solutions, and finds that the skin-friction is satisfactory in accelerated 
flow but not in retarded flow, especially in the vicinity of the true 
separation position. On the other hand, the heat-transfer rate as 
calculated using the erroneous skin-friction is exceedingly accurate 
owing to a fortuitous cancellation of errors. It is unlikely that this 
will happen in general. Illing worth has also calculated the separation 
positions for a uniform adverse external velocity gradient with zero heat 
transfer for several values of the Mach number. His solution at zero 
Mach number differs from the exact value by less than 8 per cent, and 
at MQ = 4 appears to agree with the solution calculated by the Mathe- 
matics Division, National Physical Laboratory, to within about 1 per 
cent. It is difficult to estimate what accuracy might be expected in 
a general case, in view of the respective failure and success of the 
method in predicting separation in the two examples cited. 

9. Curie's method for calculating heat transfer 

In view of the predominance of methods which aim at calculating 
skin-friction as accurately as possible, whilst sacrificing accuracy of 
prediction of heat transfer, it is valuable to conclude this chapter with 
a method which aims at calculating the heat-transfer rate when the 
skin-friction is assumed known. The method, due to Curie (1962 a) is 
an extension to high-speed flow of the method described in section 6 
of Chapter 6. 

The method begins with the Stewartson transformation, and for 
convenience the factor v^ is neglected in the definition of the stream 
function, so that the transformed velocity and temperature fields are 
given by (326) to (329). If it is assumed that the velocity field is 
adequately known by some other method, it remains to calculate the 
temperature field, which is given by the equation 

T 88_ F 0/Sf ___ 2 
+ ~ 



Following the technique used at low speeds (Curie, 1961 a), based on 
the ideas of Lighthill and Liepmann, this equation is integrated to 
yield 

(365) 



dX 





Now Lighthill indicated that for Prandtl numbers of order unity and 
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above the velocity could be taken in an approximate form valid near 
to the wall, since the thermal boundary layer is usually thinner than 
the viscous layer. Thus we write 

..., (see) 

where the first term only was retained by Lighthill, and the third term 
is negligible in the low-speed case provided temperature differences are 
small. In the case under consideration we assume that (8U/dY) w is 
known, (d*UldY*) w is obtained by setting 7 = in (327), and is 



whilst (d 3 f//c)F 3 ) M ,, obtained by differentiating (327) with respect to Y 
before setting Y -- 0, is equal to 

ldS .,, 



Thus, writing (dS/dY) w as S' w , (366) becomes 



and hence (365) may be written 



1>- 





"0 L o 



. (367) 



Following Liepmann, we write 

' n+2 

(368) 



for the case of a wall-temperature distribution which is either uniform 
or possibly uniform with zero heat transfer upstream of say X = X 19 
with a step to a new uniform value downstream of this position. Then 
(367), after formal integration with respect to X, and substitution from 
(368), takes the form 



<389) 
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where a, 6, c are positive constants, which will be determined at a later 
stage. Finally, if we write 

m LS' L M 



which is a Nusselt number in the transformed plane, (369) becomes 



- HI) - ^(S 



The three terms whose coefficients include a, 6, c are of decreasing 
orders of importance. Thus Lighthill found that upon retaining only 
the first term on the right-hand side of (370) (thus setting b and c equal 
to zero) the equation could be integrated analytically, and a formula 
obtained for F which agreed well with accurate solutions except near 
to separation, where (dUjdY) w is zero. Similarly Curie (196 la) found 
that for the case of small temperature differences (i.e. S w very nearly 
zero) the retention of the second term, involving the coefficient 6, 
enabled the accuracy to be maintained right up to separation. The 
values suggested by Curie for the two coefficients are 

a = 0-2226, b = 0-1046. (371) 

It might be expected that these two terms alone would be equally 
satisfactory in more general cases, with S w not almost zero, unless both 
almost vanish simultaneously. This could happen near to separation 
when 8 W is close to 1, that is for very low wall temperatures. Now 
examination of the Cohen and Beshotko similar solutions bears out 
this point, and Curie suggests retention of the third term to cover this 
relatively unlikely occurrence. It will be noted that the simplicity of 
equation (369) is not at all reduced by the retention of this term. The 
value of c is chosen to obtain good agreement with the accurately known 
similar solutions, and the term in square brackets in (369) becomes 



where b is given by (371) and 

b-lc = 0-0833. (372) 

The application of the method is limited strictly to the case a = 1 
(although a factor a* might well adequately convert the values of F 
into those appropriate to alternative values of a close to unity), and 
to the case of uniform wall temperature, unless the assumption be made 
that the skin-friction is little dependent upon wall temperature, when 
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the heat-transfer rates appropriate to each of a number of elementary 
steps in wall temperature may be added. Within these limits the 
accuracy is considerable, and for the similar solutions heat-transfer 
rates are given to within about i 2 per cent, for the considerable range 
of wall temperatures and pressure gradients considered. 
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INTERACTIONS BETWEEN SHOCK WAVES 
AND BOUNDARY LAYERS 

IT has already been remarked, in the first chapter of this book, that 
interactions between the boundary layer and the external stream are 
much more important in supersonic flow than in subsonic flow. Thus, 
although it is possible to treat the boundary layer and the external 
flow separately in many practical problems at low speeds, this is so 
only rarely in the supersonic case, since the presence of the boundary 
layer can make a tremendous difference to the external flow, particularly 
if the layer is a laminar one. Thus, although the inviscid theory of 
supersonic flow does not allow a disturbance at a point to have an up- 
stream influence, the boundary layer provides the means whereby an 
upstream influence can exist. 

One of the mechanisms whereby this can occur was suggested by 
Oswatitsch and Weighardt (1943), who remark that if a disturbance 
causes a positive pressure gradient the boundary layer will thicken, so 
that the external streamlines must begin to curve slightly upstream. 
Now this curvature will itself produce a pressure gradient in the same 
sense, and the boundary layer grows in equilibrium with the pressure 
gradient caused by its own growth. When the disturbance is strong 
enough to cause the boundary layer to separate the presence of a region 
of reversed flow (possibly extensive) makes a considerable difference 
to the mechanism, and, indeed, upstream spreading of a position of 
separation can occur in subsonic flow equally well. An important 
practical example of the interaction of a boundary layer with the ex- 
ternal flow is the problem of the interaction between a shock wave and 
a boundary layer. In practice the nature of the interaction will depend 
vitally upon the state of the boundary layer. Three main possibilities 
exist, namely (i) when the boundary layer is laminar throughout the 
region of interaction; (ii) when the boundary layer goes turbulent down- 
stream of the separation position, but is laminar upstream; (iii) when 
the boundary layer is turbulent throughout the interaction region. 
In this chapter it will only be possible to consider the case of boundary 
layers which are laminar throughout the region of interaction (and these 
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only very briefly). For details of the other cases, reference may be 
made to reviews by Holder and Gadd (1955) and by Gadd (1957 a). 

1. Principal results of experimental investigations 

In practical cases of shock- wave and boundary-layer interaction, such 
as those occurring on an aerofoil surface at transonic speeds, it is some- 
times difficult to separate the fundamental effects of the interaction 
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FIG. 2. Boundary layer and externally generated shock. 



Flow 
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FIG. 3. Shock generated from within boundary layer by wedge on surface. 

from those which are due to the limited extent of the supersonic region 
or to the pressure gradients associated with surface curvatures. For 
this reason the majority of the basic experimental investigations have 
been carried out with a boundary layer on a flat surface, so that the 
pressure gradients would be entirely absent if there were no shock 
waves. The shock is produced either by placing a wedge in the super- 
sonic stream above the plate (Fig. 2) or by attaching a wedge (Fig. 3) 
or a step to the surface of the plate. If 9 is the angle of deflexion through 
the external shock, then the streamlines at the edge of the boundary 
layer are turned through an angle 26 in passing through the shock and 
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the reflected expansion. A wedge of angle 26 on the plate will give 
roughly the same pressure rise, and in this case the external stream- 
lines change direction continuously but are deflected discontinuously 
through an angle 26 relative to the wall. 

It is found experimentally that the pressure distribution is roughly 
symmetrical about either the apex of the wedge on the surface or the 
position where the shock strikes the boundary layer. When separation 
occurs the pressure distribution has three points of inflexion, which 
arise for the following reasons. The pressure gradient rises from zero 
well upstream of the interaction to a maximum value at separation. 
Downstream of separation the pressure gradient decreases with distance 
due to the thickening of the dead-air region, and then increases again 
when the dead-air region becomes thinner downstream of where the 
shock strikes the layer. It increases to a maximum value at the position 
of reattachment and then gradually decreases to zero as the pressure 
tends to its limiting downstream value. 

In the case of zero heat transfer a considerable body of experimental 
evidence has been accumulated concerning the dependence of this 
pressure distribution upon such parameters as Mach number, Reynolds 
number, and shock strength. Chapman, Kuehn, and Larson (1958) 
have correlated some of these results. In particular they find that for 
a range of upstream Mach numbers 1-1 < M < 3-6, and with a Rey- 
nolds number (based on external conditions in the interaction region) 
10 4 < R < 10 6 , the pressure coefficient C p at separation takes the value 

, (373) 



the pressure coefficient being defined as 



The scatter in the experimental pressure coefficients is of order 10 per 
cent.f Various theoretical investigations give formulae like (373) with 
the numerical coefficient taking values ranging from 0-83 to 1-15. 

More recently attention has been turned towards the effects of heat 
transfer upon the pressure distribution. Early experiments (Gadd 
19576) suggested that the pressure distribution was barely affected by 
heating or cooling the wall. This conclusion did not agree with the 
theoretical treatments then available, and further experiments carried 

f A recent paper by Sterrett and Emery (1960) indicates that the numerical 
coefficient appropriate to their experimental observations decreases as M Q increases 
from 4 to 6. 
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out by Gadd and Attridge (1961) with considerably greater rates of 
heat transfer have indicated that, provided the boundary layer remains 
laminar throughout, the length of the interaction region increases with 
wall temperature, with a corresponding decrease in the pressure 
gradients. Agreement between these experimental results and the latest 
theoretical investigations is promising. 

2. Summary of early theoretical investigations 

One of the first theoretical attempts to throw some light on the 
interaction of a shock wave with a laminar boundary layer was made 
by Howarth (19486). Howarth's model consisted of two semi-infinite 
streams moving parallel to one another, one at a supersonic speed and 
the other at a subsonic speed. A small steady disturbance originates 
in the supersonic region and is propagated along the appropriate Mach 
line. At the boundary with the subsonic region it is partly reflected 
and partly transmitted. The portion transmitted into the subsonic 
region affects the whole of that part of the flow, including the interface, 
and causes further disturbances in the supersonic region which are 
propagated along the other family of Mach lines. For certain particular 
forms of incident wave, Howarth was able to calculate the details of 
the complete flow. 

Howarth's results were somewhat inconclusive, because of the as- 
sumption of an infinite depth of subsonic flow, as compared with the 
thin subsonic region in a practical boundary layer. Accordingly Tsien 
and Finston (1949) improved on Howarth's model by making the sub- 
sonic region of only finite thickness and bounded on the other side by 
a plane wall, but again no conclusive results could be obtained as no 
indication was available as to what Mach number or thickness would 
be appropriate for the subsonic region. Lighthill (19506) therefore 
further refined the analysis by replacing the subsonic layer by a region 
in which the Mach number falls continuously from M 1 in the main- 
stream to zero at the wall. Unfortunately these assumptions led to 
predictions of the upstream influence of the disturbance which were an 
order of magnitude smaller than experimental values. 

Now the source of the error in LighthilTs work, as he himself pointed 
out in later papers (Lighthill 1953 a, 6), was one of inconsistent approxi- 
mation. All three contributions discussed above assume that viscosity 
may be neglected except in the actual setting up of the basic steady 
flow, and that squares of disturbances may be neglected. Now in a 
boundary layer the pressure gradient at the wall is exactly equal to 
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the lateral derivative of the viscous shear stress, so the neglect of 
viscosity will force the pressure at the wall to remain constant unless 
a non-zero Mach number at the wall, Jf 2 , is postulated, allowing the 
pressure gradient to be balanced by inertial forces. Thus the continuous 
variation in Mach number from M l > 1 to M 2 < 1 is an improvement 
on the Tsien and Finston model, but the value of M 2 still remains to 
be found. 

Lighthill then makes the crucial point that viscous forces have not 
been entirely neglected, but only disturbances to the viscous forces. 
Now, following the ideas of Stratford (1954), we know that the viscous 
forces in the outer part of the boundary layer are almost exactly what 
they would be in the absence of pressure gradient; in other words, the 
disturbances to the viscous forces are non-negligible only in an inner 
portion of the boundary layer. Lighthill assumes (and verifies after- 
wards) that this inner layer is sufficiently thin for the approximation 
of incompressibility to be made (because velocities very near to the 
wall are low). He assumes also that the velocity may be assumed to 
be linear within this region, which is probably adequate provided that 
the disturbance is not great enough to cause the boundary layer to be 
near to separation, and finds that the value of the effective Mach number 
at the wall for use in the inviscid theory is 



where M' (0) is the lateral Mach number gradient at the wall in the 
undisturbed flow, which is proportional to the skin-friction, and L is 
a length defined in terms of the kinematic viscosity and skin-friction 
at the wall, and the wave number of the disturbance. 

Comparisons of the magnitude of the upstream influence with those 
obtained experimentally show reasonable agreement. As an example, 
the upstream disturbance falls by a factor e~ l in about 5 or 6 boundary- 
layer thicknesses when the Reynolds number is of order 10 5 to 10 6 and 
the free-stream Mach number is about 2 to 3. Accordingly it will have 
fallen to 5 per cent, in about 3 times this distance, say 15 to 20 boundary- 
layer thicknesses. 

3. Gadd's analyses for interactions causing separation 

Provided that the incident shock is strong enough the boundary layer 
will separate and the assumptions of the preceding analysis will break 
down. Alternative methods have therefore been developed for con- 
sidering shock-wave and boundary-layer interaction in cases where 
separation occurs well upstream of the agency which causes it. This 
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situation is of some considerable practical importance and can arise in 
a variety of ways, such as when a boundary layer encounters a step or 
a wedge on the surface, or when separation occurs upstream of a suffi- 
ciently strong externally generated shock. In such circumstances the 
pressure at the wall and the development of the boundary layer are 
governed by the equilibrium between the pressure gradient, the resulting 
boundary-layer thickening and the deflexion of the external stream- 
lines. The relationship between pressure and the deflexion of the 
external streamlines is derived from the theory of supersonic flow and, 
provided the deflexion angle is not too great, that given by linear theory 
will be adequate, namely 



where suffix zero refers to conditions in the undisturbed stream, about 
which the linear perturbation is taken, and a is the angle through which 
the external streamlines are deflected. In the shock- wave and boundary- 
layer interaction problem it is usual to equate the deflexion of the 
external streamlines (caused by boundary -layer thickening) to the rate 
of growth of the displacement thickness, that is 



where d{B w (x)}/dx represents the rate of growth of displacement thick- 
ness in the absence of pressure gradient, which is often neglected if the 
interaction takes place far enough from the leading edge of the boundary 
layer. Combining (374) and (375) we have 

:~{8i(*) 8io(s)}- (376) 

A second relationship connecting p and 8 X is given from the basic theory 
of laminar boundary -layer development in the presence of a pressure 
gradient, and three papers by Gadd differ mainly in the method used 
to treat this part of the problem. 

In the earliest paper (Gadd, 19536) it is assumed that in the region 
upstream of the interaction, where there is no pressure gradient, the 
velocity within the boundary layer is given by 

- = sm|^ 



where the boundary -layer thickness 8 is chosen so that the displace- 
ment thickness is in agreement with that given by Crocco's calculations. 

853502 L 2 
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If we assume that there is no heat transfer at the wall, and that a = 1 
and ^ oc T, this yields 

fi l-721(l+0-277Jf8) /v, 



{l-(1422/lf )tan~ 1 (0-4472Jf )}\ U Q ) 

In the region where the boundary layer thickens under the influence 
of the shock it is assumed that 

fi_/^\ y y [y__ 

l 



when 

\PI 

u 'u 

and = ri+^o, when < | < nX, 





( y \ 3 fi M (y+1)/2y O 
where r, = (l- ^) smj^J AJ. 

The form of r^ is determined by the condition that in the outer part of 
the boundary layer changes in the profile caused by the pressure 
gradient may be determined by essentially inviscid considerations, 
following the ideas used by Stratford (1954). The form of r 2 ensures 
continuity of u, duldy, and d^ujdy^ at the join, y nX, with the inner 
profile, and one of the two parameters n, A, is determined by the balance 
between pressure gradient and viscous stresses at the wall. The con- 
dition on mass flow (corresponding to Stratford's condition that if/ is 
continuous at the join) is also approximately satisfied provided that 
(plp Q ) l and A both remain small. The other relationship between n 
ond A, which is required to complete the solution, is given by (376). 

Gadd made a number of comparisons between his theory and experi- 
ment. In particular the predicted value of the pressure coefficient at 
separation is of order 10 per cent, greater than the mean experimental 
value (373), which is extremely satisfactory. 

In a second paper Gadd (1956) considers again the case when there 
is zero heat transfer and the Prandtl number is equal to unity, and 
assumes that the velocity profiles are given either by 

pu __ tanhm+tanh{(2//) m} 
1+tanhm 



or by a similar form which will not be considered here. The virtues of 
this family of profiles are firstly that it is simple, and secondly that for 
large values of m there is an extensive region of low- velocity flow near 
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to the wall, rather like that occurring in a well -separated laminar 
boundary layer. The main drawback is that no region of negative 
velocity occurs, so that the position of separation cannot be predicted. 
The value of I is determined from the displacement thickness in the 
region upstream of the interaction, and is assumed to be constant. 



1-f-tanhm 
so that just upstream of the interaction, where m = 0, we have 



The parameter m is determined from the condition that the pressure 
gradient at the surface is given by 

dp 



__ 

v 

which follows from the momentum equation (12) upon setting 
equal to zero for the case of zero heat transfer, and a second relation- 
ship between m and p is given by (376) with S l given by (377). Simul- 
taneous solution of these equations yields p as a function of x. Gadd 
remarks that in the case of an externally generated shock (strong 
enough to deflect the external streamlines through an angle 6) the 
streamlines are deflected instantaneously through an angle 26 at the 
position where the shock strikes the boundary layer, and this corre- 
sponds to the discontinuous change in the slope of the wall when the 
interaction is due to a wedge on the surface. 
Gadd's results may be expressed as 

C UM*-1)BV - 0-90 1 + 0-20*3 
L(M !)* _ 90 



? _ j )t . 

in which dQ/ds is a numerically defined universal function of s s, s is 
a non-dimensional value of a?, and s is the value of s corresponding to the 
upstream end of the interaction region. A comparison of the predic- 
tions of the upstream influence with the experimentally observed values 
shows fair agreement. 

In a third paper Gadd (1957c) presents a more complete treatment 
of the problem (including the effects of heat transfer) by a similar 
approach to those discussed above, but involving fewer approximations. 
For the outer part of the boundary layer, following more closely the 
analysis used by Stratford (1954) in low-speed flow, it is assumed that 
the distortion of the upstream zero-pressure-gradient profile may again 
be determined by largely inviscid considerations, and in 
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of the boundary layer (close to the wall) the profile is assumed given 
by a form which is essentially equivalent to 

pu T ( y Jy\ 2 I 

w*o ~~ T u \yj W \y>l 

The three coefficients a, 6, c, and the position of the join y^ are then 
determined by the conditions that i/r (== j pudy\, pu, d(pu)ldy, and 

d z (pu)/8y* are all continuous at the join. Thus the complete profile is 
determined in principle at any rate. In practice two further approxi- 
mations are made, namely that the inner profile is thin (so that the 
usual simplifications may be made to the joining conditions) and that 
the changes in displacement thickness under the action of the pressure 
gradient are mainly due to changes in the shape of the inner profile. 
Subject to these approximations it is indicated that the pressure at 
separation is given by an expression similar to the empirical form (373), 
but with the coefficient 0-93 replaced by 1-13. It is also predicted that 
the wall temperature T w has no influence on the pressure coefficient at 
separation, but that the stream wise extent of the region of interaction 
varies with T w as T^. As has already been remarked, the earliest 
experiments on the effects of wall temperature suggested that the 
pressure distribution is barely dependent upon wall temperature. More 
recent experiments, however, with greater rates of heat transfer, have 
shown that the stream wise extent of the interaction region does indeed 
increase with T w , though probably not as quickly as TJ,. 

4. The analysis of Hakkinen, Greber, Trilling, and Abarbanel 

An interesting analysis, based on similar ideas to the above, has been 
given by Hakkinen et al. (1959). The approximations made, though 
numerous, are all physically sound, and the analysis deserves mention 
because of its essential simplicity. 

It is assumed first that the boundary layer consists of an outer 
portion, differing little in 'shape' from the upstream undisturbed layer, 
and an inner sub -layer of small momentum. The join between the layers 
occurs at y = Ay say, and the apparent displacement of the undisturbed 
layer is A, with A ^ Ay. Two consequences of this assumption are that 

/ZW(Ay) == T (Ay A) 

and r(Ay) == T O , 

where r is the local viscous stress and T O the skin-friction immediately 
upstream of the interaction. Since the momentum therein is small, it 
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may be assumed that pressure and viscous forces balance throughout 
the inner layer. Thus 

dr dp , . A 

=2: .. when y < A,- 
By dx y 3 

and hence T (&) - 



where r w is the local skin-friction and, more generally, 

^ u==T ^+l-^ 2 when y<^- 

The conditions of continuity in u and r at the join, y = Ay, may then 
be shown to yield / _ x 2 



If it is assumed that the pressure rise due to the thickening of the 
boundary layer is main]y due to the growth in A, then 

2 ^ A 



and elimination of A between (378) and (379) yields 

X 

To (JfJ-l)tJ J C p (fa = (T O -T W ), 

cx) 
C 

or i/>o!T (JlfJ-l)^ J C, fa = (T,-T W ). (380) 

co 

Since it is found experimentally that dG p /dx is very nearly constant, 
except in the region where C p is very small, we may write approxi- 
mately x 

^2 j C p dx^lC* 9 (381) 

00 

so that (380) becomes 



or C p (Ml-l)*(lc f tf ^ c fQ -c fw , (382) 

where c /0 = ^- and c fw = -^H 



are skin-friction coefficients. Thus, at separation, (382) yields 



and, if p. oc y, it follows from (270) that this becomes 

C ps {(M*-l)R}i = (1-328)* = 1-15, (383) 
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where R is the local Reynolds number based on conditions at the edge 
of the boundary layer. This value of the pressure coefficient at separa- 
tion is about 20 per cent, higher than the mean experimental value, 
so the agreement is fair. This is probably partly fortuitous, since the 
approximation (381) alone could introduce an error of order 20 per 
cent, in C p , so some cancellation of the various sources of error has 
presumably been present. 

These same authors have also carried out a similar analysis, making 
direct use of the momentum integral equation, and Greber (1960) has 
shown how to take account of convex wall curvature. This analysis 
has been somewhat improved upon by Gadd (Bray, Gadd, and Woodger, 
1960), and will now be discussed. 

5. The work of Gadd and Greber 

The description given here follows that of Gadd, and will be restricted 
to the case of solid walls, although Gadd in fact developed the method 
in a form which makes it applicable when suction or injection takes 
place. Gadd writes the displacement thickness 3 1 as 



(384) 

-here , 



m 

and T H is the total temperature. If, for simplicity, it is assumed that 
to a good approximation we may write 

(385) 

which holds when there is zero heat transfer and when there is zero 
pressure gradient (provided a 1), then (384) becomes 



2 , (386) 

where r 2 = i f [ 1 - - ) dy. 

S 2 J PI\ %/ 
o 

This form will be accepted in what follows. 

Now in an interaction which takes place over a distance which is 
small compared with the distance from the leading edge, changes in 
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displacement thickness are likely to be large compared with changes 
in momentum thickness. This is borne out by the experimental results 
and was surmised by Gadd by inspection of the momentum integral 
equation. Accordingly it may be assumed that 8 2 is approximately 
constant, as is M . Then it follows from (386) that 



so that Si-S^ ^ s (r 2 -r 20 ) S 20 , (387) 

J o 

where suffix zero refers to (constant) values in the region upstream of 
the interaction. 

Gadd further deduces from the momentum equation (12) that 



r^ ( 388 ) 

UlA, U% J. w 

where g= ~( 1 

\_dz ^I/JM, 

and 
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He then assumes that g is a function only of r 2 , being independent of 
such parameters as the ratio of wall temperature to mainstream tem- 
perature, and deduces from (376) and (387) that 



Tx (Ml-l)*d^ lw 
By equating this value of dp/dx to that given by (388) it follows that 



dx* 

Multiplication by 2(dr 2 /dx) and formal integration with respect to x 

yield fl 

*3 = ^".(JQ-l)* f (f } ^U (39Q) 

& Tj y^oSto J ?U I l ' 



so that ^r 2 /^ is known as a function of r 2 . This relationship is universal 
except for the scaling factor, which is such that the extent of the region 
of interaction is proportional to T w . This is in better agreement with 
experiment than the earlier predictions, in which the extent of the 
interaction was found to be proportional to T^. 



154 INTERACTIONS BETWEEN 10.5 

Gadd assumes that the universal relationship between g and r 2 is 
precisely that appropriate to the similar solutions with zero heat trans- 
fer, and deduces from (390) the relationship between r z and dr^dx, 
after which the value of r 2 (or p) as a function of x easily follows. He 
finds, in particular, that the value of the pressure coefficient at separa- 
tion is equal to ^ = () . 94{(Jf 2_ l)B} - lf (39] ) 

which differs from the mean experimental value by only 1 per cent., 
in the case of zero heat transfer at least, over a range of Mach numbers 
up to about 4. 

Gadd has also given a second method which does not involve the use 
of the approximate temperature relationship (385), but rather assumes 
that the function g is dependent upon r 1 and T. u ., the dependence being 
given by the similar solutions with heat transfer. The results are 
broadly the same as those given by the method discussed above but 
indicate slightly less dependence upon T H , in the case of the extent of 
the interaction region, and a slight dependence of the pressure coefficient 
at separation upon T w . 

In the same paper Bray, Gadd, and Woodger have also given a simpli- 
fied outline of a method originally due to Crocco and Lees (1952), and 
have discussed solutions of the shock- wave and boundary-layer inter- 
action problem obtained by that method, making use of the National 
Physical Laboratory DEUCE computer. This work will not be discussed 
here, partly for reasons of space, and partly because it is the author's 
view that, for investigations of this type of interaction, the method is 
less convenient than others discussed in this book. 

6. The method of Curie 

We conclude this book with a brief description of a method due to 
Curie (19616), which bears considerable similarity to that discussed in 
the previous section and, indeed, was developed almost simultaneously. 
The method begins with the approximate equation (341 ) for the momen- 
tum thickness of a compressible laminar boundary layer, from which 
it is deduced that if the pressure rise is a short sharp one (as has already 
been seen to be often true in interactions between shock waves and 
laminar boundary layers) then it automatically follows that approxi- 
mately a _ , , 
J o 2 constant. 

Under the same conditions it follows from (345) that 

81-810 = 8 20 J'W)-#'(0)}, (392) 
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where A' is defined by (343) and H f (X') is a universal function. From 
this definition of A' it follows that when the pressure (and velocity) 
changes are small we have 

V _~_ A > _. ^AT W dp 



Then elimination of 8 X and A' between equations (376), (392), and (393) 
leads to a single equation for p or u v This equation may be expressed 
non-dim ensionally by writing 

/ ?/ 2 \ 

F = 0-4096(Jf g- 1)1^1 1--4 , (394) 

\ U Q! 



X - l-820(Jf*-l)U-l, (395) 

*W\ X 8 I 

and A -- l-6380jR*(Jfg-l)-l, (396) 

where X R represents the position of separation (or any other convenient 
position within the region of interaction), the Reynolds number R is 



and 20 is the angle of deflexion of the streamlines relative to the wall 
caused by either the shock (Fig. 2) or the wedge (Fig. 3). We then find 
that 

[ F dX = H f (F f )-H f (0)+A(X-X 1 ), (397) 

oo 

where e = upstream of the position X X where the shock strikes 
the boundary layer, and e = 1 downstream of this position. Solution 
of (397) yields F and F' as functions of X. This is done by multiplying 
the equation by F" and formally integrating with respect to X from 
oo to X, using integration by parts where necessary. This yields 

(F-~Ae) 2 = 2F'H'(F')-2 j H'(z)dz = P(F'), (398) 

o 

where P(F') may easily be determined once H' is known. Thus F is 
known as a function of F', from which we easily deduce X as 

z=f^. 



The condition that F is continuous at X ~ X l implies, by virtue of 
(398), that 
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It is then easy to see that in taking the square root of equation (398) 
we must write p = + / P( j/m x < X 1 } 

^ l V n (399) 

and F = A-{P(F')}l, X > Xj 

It will be noted that the pressure distribution is symmetrical about 
X = X v The results are shown in Table 17, where the relationship 
between F, F', and X is given for X < X lf 

TABLE 17 

Universal functions for shock-wave and boundary -layer 
interactions 



X" 

0-0103 0-02 -7-03 

0-0237 0-05 -5-12 

0-0351 0-08 -4-09 

0-0479 0-12 -3-14 

0-0612 0-17 -2-21 

0-0736 0-23 -1-32 

0-0832 0-29 -0-55 

0-0900 0-338 

0-0885 0-40 -fO-70 

0-0828 0-50 1-86 

0-0645 0-60 3-21 

0-0465 0-70 5-03 

0-0323 0-80 7-61 

0-0174 0-90 11-75 

0-0101 0-95 15-52 

0-0042 1-00 23-33 

1-03 oo 

At separation F' = A' = 0-090, and from Table 17 we note that 

F 8 = 0-338, 
and we may then deduce from (394) that 

C ps {(Ml~I)R}* = 0-825. (400) 

This value is in good agreement with the mean experimental value (373) 
for zero heat transfer. That C ps is independent of wall temperature is 
also predicted by the method described in the preceding section, due to 
Gadd. 

It follows from (395) that the length of the region of interaction is 
directly proportional to T w , which result is also predicted by Gadd's 
method, and is in at least reasonable accord with the most recent ex- 
perimental results. 

A similar approach has also been used (Curie, 19626) to consider the 
interaction when a laminar boundary layer meets an expansive corner 
followed by a compressive agency, such as a step. 
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